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and  Ethylene  Dibromide 

C-21 

C-9 . 

Potential  Atmospheric  Sources  of  Bromine  in 

1975  (in  106  kg  Br./yr) 

C-21 

CA-1 . 

Estimated  World  Chlorocarbon  Emissions,  1973 

0-24 

D-l . 

Major  Volcanic  Eruptions,  1955  to  1976 

D-4 

D-2. 

Stratospheric  Source  Function  of  the  Agung 
Eruption  Cloud 

D-5 

D-3. 

A Recommended  List  of  Stratospheric  Injections 

Due  to  Major  Volcanic  Eruptions  From  1 9 55— 1 9 7 5 a 
Scaled  in  Terms  of  the  Agung  Injection  of  Table 
D-2 

D_  5 

E-l . 

Effective  Area  of  the  Polar  Caps  C,  and  C of 

Fig.  E-2 

E-ll 

E-2 . 

Significant  Solar  Proton  Events,  1955-1973 

E-  2 

E-3. 

Some  Large  Solar1  Proton  Events  and  Their 
Ionisation 

E-26 

EA-1. 

Photochemical  Lifetimes  of  NO  and  HO 

X X 

E-3  5 

EA-2 . 

Representative  Diffusion  Time  Steps  Through 

5-km  Height  Intervals 

E-3  7 

EB-i. 

Some  Major  Solar  Cosmic  Ray  Observations 

During  1956-1961 

E—  4 2 

EB-2 . 

Relation  Between  Rigidity  and  Kinetic  Energy 
for  Protons 

E—  4 3 

EB-  3 . 

Integral  Exponential  Spectra  for  Solar 

Frotons — EQ  (D-3) 

E-4  4 

EB- A . 

) _ ' 

Solar  Protons  (cm  -sec)-  of  Energy  Greater 
than  20  MeV  at  Various  Times  Following  17 

Solar  Flares 

E-49 

EB-5. 

Jhar&e't eristic  Rigidity,  P0(t'i,  and  Ratios  of 
Proton  Fluxes  > 3 MeV,  >50  MeV,  and  > 100  MeV, 
to  that  > 20  MeV  for  Various  Times  in  an 
"Average"  Solar  Proton  Event 

E-50 

x iv 


Sources  of  Fixed  Nitrogen,  ca  1975 


World  Annual  industrial  Production  of  Fixed 
Nitrogen,  1955  to  1975 

Global  Production  and  Consumption  of  Fertil- 
izer Nitrogen,  1953  to  197^ 

Estimate  of  Anthropogenic  Nitrogen  Fixation 
Rate,  1955  to  1975 


SUMMARY 


Here  we  compile  a catalog  of  time-dependent  natural  and 
anthropogenic  sources  of  NO  , HO  , and  C1X  during  the  time  frame 

X X 

1955-1975  as  well  as  other  time-dependent  perturbations  that  would 
be  expected  to  give  rise  to  time  variations  of  total  ozone  other 
than  the  diurnal  and  seasonal  fluctuations.  Factors  considered 
include  atmospheric  nuclear  explosions,  halocarbons,  volcanic 
eruptions,  solar  and  galactic  cosmic  rays,  industrially  fixed 
nitrogen,  solar  ultraviolet  (UV),  and  atmospheric  motions.  The 
catalog  presents  atmospheric  or  stratospheric  inputs  as  a function 
of  time  and  latitude,  to  provide  an  appropriate  input  for  two- 
dimensional  computer  models  so  that  one  can  establish  how  well 
the  existing  models  explain  the  actual  ozone  record. 

Figure  S-l  lists,  in  schematic  form,  the  major  sources  con- 
sidered, and  Table  S-l  scales  the  injections  due  to  each  source 
in  terms  of  the  ambient  stratospheric  loading  of  NO  , HO  , and 

X X 

C1X  (see  Appendix  A). 

In  terms  of  current  (1977-1978)  ideas  of  stratospheric 
photochemistry,  which  suggest  that  the  relative  importance  of 
C1X  as  against  N0x  injections  on  ozone  is  greater  than  was 
thought  in  1974-1975,  the  effects  of  halocarbon  emissions,  ex- 
plosive volcanic  eruptions,  and  atmospheric  nuclear  explosions 
will  be  most  important,  with  smaller  effects  due  to  cosmic  rays 
and  to  nitrogen  fixation.  This  refers  to  the  injection  of 
chemical  species  into  the  stratosphere,  but  the  biggest  single 
effect  on  stratospheric  ozone  may  well  arise  from  the  variation 
of  solar  UV  irradiance  with  the  11-year  sunspot  cycle.  However, 
the  interaction  between  NO  , HO  , and  C1X  species  in  the  atmosphere, 

X X 
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FIGURE  S - 1 . Important  Happenings  During  the  1955-1975  Time  Frame 


Estimate  of  ground-level  Injectants  that  are  sufficiently  long-lived  to  enter  the  stratosphere. 
6?  of  total  nitrogen  fixation  representing  N„0  which  can  enter  the  stratosphere  and  produce  NO 
by  reaction  with  0(  D) . 


For  the  convenience  of  the  reader,  the  work  is  discussed 
and  the  results  are  presented  in  two  levels  of  detail.  The 
first  level,  which  includes  recommendations  for  inputs  due  tc 
specific  sources,  is  given  in  Section  2,  followed  by  a discus- 
sion in  Section  3 of  some  other  factors  that  may  be  important, 
in  particular,  of  changes  in  solar  UV  irradiance  with  sunspot 
number  (see  Section  3-1)  > but  which  have  not  been  considered  ir. 
as  much  detail  as  the  other  items.  Detailed  review's  of  the 
appropriate  stratospheric  or  atmospheric  injections  between  1955 
and  1975  are  given  as  Appendices  B (atmospheric  nuclear  explo- 
sions), C (halocarbons ) , D (volcanic  eruptions),  E (cosmic  rays), 
and  F (nitrogen  fixation).  Appendix  A provides  estimates  of 
stratospheric  burdens  of  the  various  trace  species  of  interest, 
to  provide  a crude  indication  of  the  relative  significance  of 
the  various  perturbing  sources. 


1.  INTRODUCTION 


NO  + NO  + ...),  oxides  of  hydrogen  ("HO  " = 

c.  X 


It  is  known  that  stratospheric  ozone  is  destroyed  by  a 
variety  of  catalytic  reaction  cycles  involving  oxides  of  nitro- 
gen ("MO  " 
x 

OH  + HO^  + •••),  and  chlorine  species  ("C1X"  = Cl  + CIO  + HC1 
+ CIO  N02)>  also  bromine,  "BrX,"  and  possibly  also  other  species 
The  stratospheric  concentrations  of  ozone,  NO^,  HOx , and  C1X, 
are  affected  by  a variety  of  natural  and  anthropogenic  sources 
such  as  solar  UV,  solar  and  galactic  cosmic  rays,  volcanic 
eruptions,  atmospheric  nuclear  explosions,  halocarbons,  and 
(industrially  and  naturally)  fixed  nitrogen.  Other  geophysical 
factors,  such  as  possible  variations  in  stratospheric  temper- 
ature and  in  stratospheric  motions  may  also  affect  stratospheric 
ozone.  Because  of  the  interaction  of  the  NO  , HO  , and  C1X 

X X 

species  with  one  another,  as  well  as  with  ozone,  with  radiation 
and  dynamics  it  is  necessary  to  incorporate  all  these  factors 
simultaneously  in  a photochemical/dynamic  computer  model  to  eval 
uate  the  total  effect. 


There  exist  data  on  ozone  as  a function  of  altitude,  lati- 
tude, and  time  with  moderately  adequate  global  coverage  since 
the  mid-1950's,  by  which  time  a global  network  of  some  45  to  90 
Dobson  spectrophotometers  and  M-83  filter  photometers  had  been 
established;  since  the  1970  launch  of  Nimbus  IV,  there  have 
been  extensive  satellite  measurements  which  will  provide  a 
more  detailed  and  fine-grained  record  once  the  data  are  cali- 
brated and  reduced  so  that  they  fit  into  the  long-term  record. 
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Long-term  global  trends  of  total  ozone  have  been  evaluated 
by  a number  of  workers,  in  particular  by  Komhyr  et  al . (1971), 
Johnston  et  al.  (1973) , London  and  Kelley  (197*0;  see  in  partic- 
ular Angell  and  Korshover  (1973,  1976,  1978)  and  work  referenced 
therein,  also  the  detailed  atlas  of  London  et  al.  (1976).  Angell 
and  Korshover  (1978)  find  a significant  5 percent  increase  in 
total  ozone  in  the  Northern  Hemisphere  between  i960  and  1970, 
a marginally  significant  change  thereafter.  In  the  Southern 
Hemisphere,  there  appeared  to  be  a 2 percent  increase  in  total 
ozone  between  i960  and  1968,  and  a 1 percent  decrease  thereafter; 
see  Pig.  1.1,  which  demonstrates  both  the  temporal  variations  and 
also  the  error  bounds  in  total  ozone  data.  The  other  workers 
find  comparable  effects. 

All  these  measurements  use  the  strong  Hartley-Huggins  ab- 
sorption of  ozone  below  300  nm.  Angione  et  al.  (1976)  have 
analyzed  astronomical  records  in  the  Smithsonian  astrophysical 
observation  program  between  1912  and  1950  which  cover  the  weak 
Chappuis  absorption  of  ozone  in  the  500-  to  700-nm  band.  There 
appear  to  be  some  long-term  variations,  but  it  is  not  apparent 
that  any  inference  on  global  trends  can  be  made  prior  to  1955. 
Reference  should  be  made  here  to  the  very  careful  analysis  of 
the  ozone  record  at  Arosa  from  1926  to  1971  by  Birrer  (1974), 
who  finds  at  this  single  station  that  trend  values  decrease 
with  increasing  periods  of  observation  and  raises  warnings 
against  inferring  trends  from  an  inadequate  period  of  observations 

The  basic  reason  for  compiling  this  catalog  is  to  provide  a 
check-case  for  two-dimensional  photochemical  models  of  the 
stratosphere,  to  see  how  well  a given  model  can  reproduce  the 
time-dependent  ozone  data.  Because  of  this  emphasis  on  time- 
variations  in  ozone  (on  a scale  greater  than  several  months), 
steady-state  sources,  such  as  the  continuous  natural  injections 
of  N»0,  chloro-  and  bromocarbons  from  the  biosphere,  as  well  as 
the  steady-state  production  of  HO  arising  from  natural  water 
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Source:  Angel!  and  korshover,  1978 

FIGURE  1.1.  Temporal  Variation  in  Total  °zone.  Expressed  as  a percentage 
Deviation  from  the  Mean  for  Northern  and  Southern  Hemispheres 
and  che  World.  The  annual  oscillation  has  been  eliminated  and 
the  data  smoothed  by  application  of  a 1-2-1  smoothing  twice  to 
successive  seasonal  values  (last  data  are  the  northern  summer 
of  1976).  The  vertical  bars  represent  confidence  limits  such 
that  there  is  about  a 5 percent  chance  that  the  true  value  of 
the  mean  lies  outside  the  vertical  extent  of  these  bars.  The 
vertical  arrows  indicate  the  time  of  quasi-biennial  west  wind 
maximum  at  50  mbar  (21  km)  in  the  tropics. 
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transport  into  the  stratosphere  are  not  considered  here.  For  an 
overall  review  of  these  problems,  see,  for  instance.  National 
Academy  of  Sciences  (1976,  1977),  Hudson  (1977)}  and  Oliver  et  al. 
(1977) . 
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2.  RECOMMENDED  INJECTIONS 

2.1  INTRODUCTION 

This  section  consists  of  a list  of  inputs  for  two-dimen- 
sional photochemical/dynamical  computer  models,  extracted  from 
the  detailed  discussions  listed  in  Appendices  B to  F.  Thus,  an 
explanation  for  the  choice  of  numerical  values  recommended  here 
is  given  in  the  appropriate  appendices, 

2.2  ATMOSPHERIC  NUCLEAR  EXPLOSIONS 

Table  2.1  (Table  B-12  of  Appendix  B)  lists  estimates  of 

stratospheric  injections  of  NO  due  to  past  atmospheric  nuclear 

x 

explosions.  The  N0x  injection  rate  used  here  is  1.0  x 10J 
molecules/Mt  yield,  as  against  a best  estimate  range  of  (0.4  - 
1.5)  or  (0.6  - 1,1),  x 10^  molecules/Mt  (see  Table  B-l  of 
Appendix  B).  The  latitude  of  given  injections  may  be  found  in 
Table  2.2  (Table  B-6  of  Appendix  B) . The  distribution  of  in- 
jection heights  as  a function  of  yield  and  latitude  comes  from 
Fig.  B-5  of  Appendix  B;  within  each  given  altitude  range  the  dis- 
tribution of  Table  2.3  (Table  B-2  of  Appendix  B)  is  suggested. 

Note  that  submegaton  yield  explosions  are  not  considered 
(see  Section  B.5  of  Appendix  B)  , and  the  injection  of  H.^O  and  of 
0^  is  unlikely  to  be  of  critical  importance  in  the  present  con- 
text (see  Section  B.4  of  Appendix  B) . Also,  detrainment  is  not 
specifically  taken  into  account  (see  Section  B.3  of  Appendix  B)  ; 
its  effect  would  be  to  reduce  the  effective  injection  rate  of 
N0x  per  Mt , but  the  quantitative  change  is  unlikely  to  be  sig- 
nificant in  the  context  of  other  uncertainties. 
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TAL5LE  2-1 


STRATOSPHERIC  INJECTION  OF  NOx  BY  NUCLEAR  EXPLOSIONS 


TABLE  2-1.  STRATOSPHERIC  INJECTION  OF  NOx  BY  NUCLEAR  EXPLOSIONS  (cont) 


Date 

Country 

Location 

4 Nov  61 

USSR 

N? 

2 May  62 

USA 

XI 

10  Jun  62 

USA 

XI 

27  Jun  62 

USA 

X I 

30  Jun  62 

USA 

XI 

11  Jul  62 

USA 

XI 

5 Aug  62 

USSR 

N7 

20  Auq  6? 

USSR 

N7 

22  Aug  62 

USSR 

NZ 

25  Aug  62 

USSP 

N7 

27  Aug  62 

USSR 

NZ 

8 Sep  62 

USSR 

NZ 

15  Sep  62 

USSR 

NZ 

16  Sep  62 

USSP 

NZ 

18  Sep  62 

USSR 

NZ 

19  Sep  62 

USSR 

NZ 

21  Sep  62 

USSP 

NZ 

25  Sep  62 

USSR 

NZ 

27  Sep  62 

USSR 

N7 

18  Oct  62 

USA 

J1 

22  Oct  62 

USSR 

NZ 

30  Oct  62 

USA 

JI 

23  Dec  6? 

USSR 

N 7 

25  Dec  62 

USSR 

NZ 

17  Jun  67 

CHINA 

IN 

24  Aug  68 

FRANCE 

CEP 

8 Sep  68 

FRANCE 

CEP 

27  Dec  68 

CHINA 

IN 

29  Sep  69 

CHINA 

IN 

30  May  70 

FRANCE 

CEP 

30  Jun  70 

FPANCE 

CEP 

14  Oct  70 

china 

i n 

14  Aug  71 

FRANCE 

CEP 

26  Jul  73 

CHINA 

LN 

17  Nov  76 

CHINA 

IN 

NOTES:  For 

1 :•  cations 

see  Table 

?’*  (USSR 

, 1953-195: 

Arc  (USSR 

* 1958;  tal 

Tor 

effective 

yields,  S ( 

actu 

a 1 .’eld. 

NO  inject 

For 

*»"  . '-njtct 

ion,  j e e T i 

p Idu 

• ’Me  Hoi 

ce . 

For 

i ’ t i ? j d e r a * ’ e , see  1 

r ete 

» s<  • • ^ 1 9 

’ C ' Curve-. 

' - 

n tnnre. 

Table  p - 2 

\0 

' r J e c 1 1 o n 

a • * k t ro 

Effective 
Me  1-1, 


NO*  Injection 
1Q^  mo  is 


Altitude  r'l H'H- 
10  - 20.7 


■ - r L’Olar  bursts.  pitz  ' I t ,•  j- 
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TABLE  2-2.  TEST  SITES 


U.S.A. 


Eniwetok  (EN) 

11° 

N, 

162° 

E 

Bikini  ( B K ) 

11° 

N, 

165° 

E 

Chris tmas  Island  (XI) 

2° 

N , 

157° 

W 

Johnston  Island  ( J I ) 

17° 

N , 

169° 

w 

U.S.S.R. 

Novaya  Zeml ya/Arc ti c (NZ) 

75° 

N, 

55° 

E 

Si beria 

52° 

N, 

78° 

E 

Semi palatinsk  (SP) 

50° 

N , 

80° 

E 

U.K. 

Christmas  Island  (XI) 

2° 

N, 

157° 

W 

FRANCE 

Centre  d ' Exper i men ta ti on  du 
Pacifique  (CEP),  Tuamotu 
Archipelago 

21° 

s, 

1 37° 

W 

CHINA 

Lop  Nor  ( LN ) 

O 

o 

N, 

90° 

E 

TABLE  2-3.  INITIAL  VERTICAL  DISTRIBUTION  OF  RADIOACTIVITY 
OR  NOx  ASSUMED  WITHIN  "MUSHROOM  CAP" 


Layer,  Fraction  of 
"Mushroom  Cap" 
(from  base  to  top) 

n 1 


Percentage  of  Activity 
Within  Layer 


1 

7 

2 

7 

3 
7 

4 
7 

5 
7 

6 
7 

Source : 


2 

7 

3 
7 

4 
7 

5 
7 

6 
7 

7 

7 

Peterson,  1970 
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1 


14 


25 


25 


1 5 


1 5 


2.3  HALOCARBONS 

Atmospheric  injection  rates  per  year  for  fluorocarbons  11 
(CFCl.^,  12  (CF2C12),  and  22  (CHF,C1)  are  listed  in  Table  2.4 
(Table  C-2  of  Appendix  C),  and  release  rates  of  various  chloro- 
carbons  are  given  in  Table  2.5  (from  Tables  2-3,  2 —A , and  2-5  of 
Appendix  C).  These  are  all  ground-level  injections;  for  the 
latitude  range  of  injection,  we  suggest  Option  I,  with  Option 
II  as  a secondary  alternative: 


Option 

I 

Option  II 

50-60° 

N 

5 % 

30-55°  N 90% 

30-50° 

N 

90% 

O 

o 

on 

1 

s 

2% 

30-40°  S 4% 

30-40° 

s 

3% 

where  the  emission  should  be  distributed  uniformly  within  each 
relevant  latitude  band. 

During  the  time  period  of  interest  here,  the  anthropogenic 
release  of  bromocarbors  is  so  small  relative  to  the  natural  pro- 
duction that  it  should  not  be  considered  (see  Section  C.4  and 
especially  Table  C-9  of  Appendix  C). 


2.4  VOLCANIC  ERUPTIONS 

The  volcanic  eruptions  to  be  considered  are  listed  in  Table 
2-6  (Table  D-3  of  Appendix  D) , and  the  suggested  source  function 
of  Agung  is  given  in  Table  2-7  (Table  D-2  of  Appendix  D). 

Note  there  is  a large  uncertainty  in  both  relative  and  ab- 
solute amounts  of  injectants,  and  also  in  the  effective  injection 
height.  This  simply  reflects  our  lack  of  knowledge  in  the  area. 


TABLE  2-4. 

GLOBAL 

PRODUCTION 

AND  RELEASE 

RATES  FOR 

FC-11  , 

-12,  AND 

-22  THROUGH 

1976  (in  10 

6 kg) 

FC-11 

FC- 

12 

FC 

-22 

fear 

Production 

Rel ease 

Production 

Rel ease 

Producti on 

Release 

1931 

0.5 

1932 

0.1 

1933 

0.3 

0.1 

1934 

0.7 

0.2 

1 935 

1 .0 

0.2 

1936 

0.1 

1 .7 

0.4 

1937 

0.1 

3.1 

0.6 

1 938 

0.1 

2.8 

0.9 

1 939 

0.1 

0.1 

3 . 9 

1.8 

1940 

0.2 

0. 1 

4.5 

1 . 7 

1 941 

0.  3 

0.  1 

6.3 

2.3 

1 942 

0.3 

0.1 

5.9 

2.9 

1 943 

0.4 

0.2 

8.2 

3.5 

1 944 

0.4 

0.2 

16.7 

4.7 

1945 

0.4 

0.3 

20.1 

6.  1 

0.1 

1946 

0.7 

0.6 

16.6 

11.9 

1 947 

1 . 3 

1.2 

20.1 

19.0 

0.1 

1948 

3.0 

2.3 

24.8 

22.2 

0.2 

1 949 

4.5 

3.7 

26.1 

24.2 

0.3 

0.1 

1 950 

6.6 

5.4 

34.9 

27  . 4 

0.8 

0.1 

1951 

9.1 

7.5 

37.0 

30.8 

1.0 

0.2 

1 952 

13.7 

10.9 

38.5 

32.6 

1 .6 

0.3 

1953 

17.6 

14.9 

48.2 

37.0 

2.2 

0.4 

1954 

21.4 

18.7 

51  . 3 

42.2 

2.9 

0.6 

i 955 

27.0 

23.2 

60.3 

47.5 

3.5 

0.8 

1956 

33.5 

29.1 

71.9 

55.4 

6.1 

1 .2 

1957 

35.3 

32.8 

77.9 

63.1 

6.4 

1 . 1 

1 958 

31  . 3 

31  .2 

77.8 

66.4 

7.7 

1 . 4 

1959 

37.9 

32.3 

92.5 

73.9 

11.4 

2.7 

1 960 

52.5 

42.2 

105.0 

88.1 

12.3 

3.3 

1 961 

63.9 

54  . 1 

114.7 

98.6 

12.3 

3.6 

1962 

82.  1 

67.6 

135.0 

113.1 

15.6 

4.7 

1 963 

98.0 

82.6 

154.0 

132.5 

17.8 

5.7 

1964 

116.6 

97.9 

178.4 

153.8 

22.6 

6.9 

1 965 

129.1 

111.7 

199.2 

173.6 

25.4 

8.2 

1 966 

148.  1 

125.1 

226.1 

192.9 

31  . 9 

9.9 

1967 

167.8 

141.9 

2 53.5 

217.6 

37.8 

12.5 

1 968 

192.1 

161.6 

279.3 

248.9 

46.1 

16.5 

1969 

227.3 

187.0 

310.0 

271.1 

56.6 

20.6 

1970 

249.2 

212.3 

334.8 

296.1 

59.3 

23.5 

1 971 

275.5 

233.2 

356.3 

317.7 

65.2 

26.5 

1972 

315.8 

260.3 

391.1 

34  3.5 

70.8 

31  . 3 

1973 

359.5 

295.6 

434.  7 

378.0 

77.4 

36.1 

1 974 

379.7 

324.1 

455.0 

408.4 

88.7 

42.3 

1975 

326.5 

314.9 

396.1 

396.3 

74.4 

47.6 

1976 

355.2 

308.2 

427.9 

380.1 

94.1 

60.2 

For  FC-11  and  -1?  the  data  from  Table  5 of  DuPont  (1977).  For  FC-22,  the  production 
figures  are  fron  '.rant  (1977)  corrected  for  global  production  as  indicated  in  the 
text  (Se.tion  2).  The  FC-22  release  figures  are  obtained  from  the  production  figures 
by  multiplying  by  the  rel ease/production  ratio  from  Table  7 of  McCarthy  et  al . (1977) 
on  a year-by-year  basis;  but  for  1976  the  1975  ratio,  47.1/73.6  = 0.640,  was  used. 
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Data  for  carbon  tetrachloride  come  from  Table  C-4,  and  for  methyl  chloroform  from  Table  C-6.  In  each  case, 
we  list  here  the  arithmetic  mean  of  all  the  release  figures  quoted  in  the  approximate  table.  All  the 
other  figures  coine  from  Little  (1975)  using  the  analysis  of  Annex  A to  Appendix  C. 


TABLE  2-6.  A RECOMMENDED  LIST  OF  STRATOSPHERIC  INJECTIONS  DUE  TO 


TABLE  2-7  STRATOSPHERIC  SOURCE  FUNCTION  OF  THE 
AGUNG  ERUPTION  CLOUD 


Mass,  kg 


Estimate  of 
Cadle  et  al . 
(1977) 

Estimate  of 

S.  Self  (private 
communication, 

May  1978) 

Suggested  h 
Value 

c 

Total  erupted  material 

S x ,o» 

1.4  x 1012 

12 

2.6  x 10 

Fine  ash  d 

o 

P~"o 

X 

9 

9 x 10 

9.5  x 109 

Total  gas 

- 

2.9  x 1010 

- 

Water  vapor  e 

1.2  x 1011 

2.2  x 1010 

5.1  x 1010 

Total  S as  H2S04  f 

1.8  x 1010 

2.1  x 109 

6.1  x 109 

Hrl  9 

1.2  x 109 

7.3  x 108 

9.4  x 108 

d Based  on  17  March  and  16  May  eruptions. 
b Geometric  mean  of  the  Cadle  et  al . and  Self  estimates. 
c Most  of  this  material  does  not  enter  the  stratosphere. 

d Cadle  et  al (1977)  do  not  specify  the  particle  size;  Self  refers  to  particles 
< 2 um  diameter. 

e 12 

The  mean  stratospheric  water  vapor  loading  is  ^ 2 x 10  kg,  see  Appendix  A. 

f 9 .34 

5.1  x 10  kg  HoS0.  contain  3.3  x 10  S atoms  so  that  presumably  (3-6)  x 

34  * ^ 

10  OH  radicals  are  consumed  in  the  formation  of  sulfuric  acid. 

n 8 33 

9 9.4  x 10  kg  HC1  contain  9 x 10  molecules. 


2.5  CARBONYL  SULFIDE 

Crutzen  (1976)  suggested  that  the  gaseous  species  carbonyl 

sulfide  (COS)  may  be  a significant  source  for  the  stratospheric 

■erosol  layer  in  volcanically  quiet  times.  This  species,  which 

is  produced  in  a variety  of  ways — biologically,  in  volcanic 

eruptions,  and  in  coal  combustion — is  relatively  inert  in  the 

troposphere  in  comparison  with  the  common  atmospheric  sulfur- 

■ontaining  gases,  H S and  SO  , but  will  presumably  be  oxidized 

2 ^ I 

o sulfuric  acid  in  the  stratosphere,  and  thus  will  contribute 

to  the  Junge  layer,  and  conceivably  may  affect  stratospheric 

ozone,  cf.  Cadle  et  al.  (1975). 

Observations  of  Maroulis  et  al.  (1977)  and  of  Sandalls 
and  Penkett  (1977)  indicate  that  the  species  is  fairly  widely 
and  uniformly  distributed,  which  supports  Crutzen's  suggestion 
and  implies  either  a long  atmospheric  life  time  (compared  with 
H^S  and  SO.,)  or  uniform  terrestrial  sources  (i.e.  not  concen- 
trated in  urban  or  other  limited  areas),  or  both. 

At  present  there  exist  insufficient  data  to  make  a reason- 
able estimate  of  atmospheric  COS  input  rates,  but  its  production 
h.  coal  combustion  makes  it  at  least  a potential  future  source 
f"or  the  stratospheric  aerosol  layer  (see  Turco  et  al.  (1978) 
for  a discussion  of  this  topic). 

2.6  COSMIC-RAY-INDUCED  IONIZATION 

Cosmic  rays  produce  ionization  in  the  atmosphere;  explicitly 
Jn  the  stratosphere  one  ion  pair  gives  rise  to  1.3  NO  molecules 
end  2.0  OH  radicals  on  the  average  (see  Section  E.3,  Appendix  E). 

The  OH  radicals  are  short-lived,  so  that  vertical  transport  for 
distances  greater  than  about  1 km  is  never  important;  below  about 
■3  km,  NO  molecules  are  sufficiently  long-lived  so  that  vertical 
transport  into  regions  where  they  can  affect  the  bulk  of  strato- 
spheric ozone  may  be  significant. 

* 

Lrs.  A.  C.  Aikin  and  J.  P.  Friend  suggested  that  this  species 
hould  be  discussed  separately. 
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For  ionization  in  the  stratosphere,  cosmic  ray  energies  in 
the  range  30  to  500  MeV  are  required,  and  thus  we  are  interested 
in  galactic  cosmic  rays  (OCRs)  and  in  solar  proton  events  (SPEs). 
Galactic  cosmic  rays  are  modulated  by  the  solar  wind  and  peak 
during  sunspot  minimum;  these  are  comparatively  energetic  par- 
ticles whose  ionization  peaks  near  15  km  [see  Fig.  2.1  (Fig.  E— 6 ) ] . 
By  contrast,  solar  protons  are  produced  by  active  solar  regions 
which  give  rise  to  intense  solar  flares  near  solar  maximum;  the 
kind  of  large  event  which  affects  the  stratosphere  occurs  perhaps 
1 to  5 times  per  11-yr  solar  sunspot  cycle  and  has  a duration  of 
1 to  3 days  on  the  average.  Figure  2.1  also  shows  the  ioniza- 
tion profiles  of  some  large  SPEs.  Note  that  the  mean  energy  of 
solar  protons  is  significantly  lower  than  that  of  galactic  cos- 
mic rays.  Table  2.8  lists  those  large  SPEs  during  the  1955- 
1975  time  period  that  should  be  considered  in  the  present  con- 
text . 

Cosmic  rays  are  affected  by  the  earth's  magnetic  field,  and 
thus  their  ionization  is  particularly  important  in  the  polar  caps 
or  regions  of  high  geomagnetic  latitude  A.  Figure  2.2  (Fig.  E-2) 
shows  both  the  "outer"  polar  caps  CQ  which  correspond  roughly  to 
A % 60°,  and  also  the  "inner"  polar  caps  corresponding  A K 6 5°. 
The  relatively  high  energy  GCRs  are  significant  in  the  outer  polar 
caps  and  also  at  low  latitudes  (see  Table  2-9),  while  the  less 
energetic  solar  protons  are  only  important  within  the  inner  polar 
caps.  (See  Table  E-l  and  Fig.  E-3  for  a representation  of  both 
sets  of  polar  caps  in  terms  of  geographic  latitude  6 for  use  in 
two-dimensional  photochemical  models.) 
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TABLE  2-9.  TOTAL  IONIZATION  DUE  TO  GCRs 


32 

(in  10  ion  pairs/yr) 

Solar  Sunspot 

Solar  Sunspot 

Minimum  (i  .e. , 

Maximum  (i  .e. , 

maximum  ionization) 

minimum  ionization) 

High  Lati tudes--in  curves  CQ  13 

7 

Low  Lati  tules- -outside  curves  CQ  26 

20 

'onl  ition  profile  due  to  OCRs  may  be  expressed  as  a 
■ Lme  (relative  to  ;he  year  of  GCR  maximum, 

. . , t ■ minimum)  and  of  latitude  (where  we  distir,- 

;ulst  : ■ . th  • Lar  cap  region,  defined  in  terms  of  the 

of  ■"  : (or  A > 1 3 and  the  rest  of  the  earth, 

■ curve  ’’1-  ' " (or  A < 60°)  as  follows  in  terms  of  equations 
of  dvtio:  1.4  in  Appendix  E. 

High  geomagnetic  latitude,  : Low  geomagnetic  latitude 

° (1"Co) 

a > 20  km  Eq  (E4.1)  2 > 31  km  Eq  (E4.4) 

1 ■■  km  < z < 2 0 km  Eq  (E4.2) 

L0  km  < z < 15  km  Eq  (E4.3)  18  km  < x < 31  km  Eq 

(E4 . 5) 

The  ionization  profile  due  to  SPEs  may  be  taken  as  uniform 

with  latitude  over  the  polar  caps  defined  in  terms  of  the  inner 

:urves  C.  of  Pig,  2.2,  or  equivalently,  for  geomagnetic  latitudes 
> 

A o,  f . The  altitude  profile  of  ionization  due  to  SPEs  is 
shown  in  Fig.  2.1  (Pig.  E-6,  Section  E. 5) . 

A rude  appreciation  of  the  overall  significance  of  cosmic 
ation  on  stratospheric  ozone  can  be  seen  from  Tables  ?.  • 
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and  2.9,  which  show  a source  of  order  20  x lO3  ion  pairs  per 
event  or  per  year.  Since  one  ion  pair  rives  rise  to  1.3  NO 
molecules  and  2.0  OH  radicals  on  the  average,  we  have  a source 
of  about  25  x 10^2  NO  molecules  and  40  x 10  OH  radicals  per 
event  or  per  year.  By  comparison: 

• The  Soviet  atmospheric  nuclear  test  series  of  Fall  1961 
(the  largest  ever,  corresponding  to  l80  Mt  total  yield) 
injected  (7C-270)  x 1 02  2 NO  molecules  into  the  strato- 
sphere . 

• The  global  stratospheric  loading  of  NO,,  = NO  + N0o  + HN0o  is 
of  order  (4-15)  x lO-'  molecules  (see  Appendix  A). 

• The  global  stratospheric  loading  of  HO  radicals  (=  OH  + 

52  x 

HO  ) is  of  order  (7-16)  x 10  molecules,  but  with  large 
uncertainties  and  with  a very  rapid  turnover  rate  (see 
Appendix  A). 

Thus  we  conclude  that  while  cosmic  ray  ionization  effects 
are  not  negligible,  they  are  not  of  primary  importance  in  the 
present  context,  as  compared  either  with  nuclear  test  series 
which  put  more  NO  into  the  stratosphere  or  in  particular  with 
sources  of  C1X  such  as  halocarbons  and  volcanic  eruptions. 

2.7  ANTHROPOGENIC  NITROGEN  FIXATION 

Recommended  time-varying  anthropogenic  nitrogen  fixation 
rates  are  listed  in  Table  2-10.  The  numbers  come  from  the  dis- 
cussion in  Appendix  F,  and  the  following  points  should  be  noted: 

a.  Commercial  nitrogen  fixation,  which  is  principally 

(>,  75  percent)  fertilizer  production,  has  been  increas- 
ing at  a rate  of  10  percent  per  annum,  much  faster  than 
the  combustion  source  (5  percent). 


During  the  1955-1975  time  frame,  there  have  been  just  three 
major  series  of  SPEs  that  should  be  considered  in  the  present 
context,  namely,  those  of  May-July  1959,  November  I960,  and 
August  1972. 
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b.  The  numbers  listed  are  production  rates  of  fixed 

nit  r -■  . 'h  r is  si  gn ifi  ;ant  iela y b<  f re  th 

fixed  nitrogen  is  denitrified,  a process  ' r:  w.u ! 
approximately  6 percent  of  total  nitrogen  goes  ':-.v 
::  ( see  ST,  1971  ) , whi  1 ■ the  rem;  inder  . *ts 

• Then  there  are  further  delays  before  this  N 
mixed  throughout  the  troposphere,  and  : f th< 

is  decomposed  in  the  stratosphere  to  product  ;J0, 
which,  in  turn,  reacts  with  ozone.  The  total  tire- 
delay  involved  here  is  probably  in  excess  of  1 .arc, 
and  thus  because  of  the  rapid  increase  in  anthropo- 
genic nitrogen  fixation  it  appears  unlikely  that  a 
significant  effect  on  ozone  would  be  apparent  as  early 
as  1970-1975.  This  point  is  discussed  in  Secticn  F-i 
of  Appendix  F. 

c.  Because  of  this  long  delay,  the  latitude  variate  n jf 
production  or  injection  should  not  be  critical,  but  see 
Table  F-3  for  the  latitude  variation. 


TABLE  2-10.  ESTIMATE  OF  ANTHROPOGENIC  NITROGEN 
FIXATION  RANGE,  1955-1957,  MT(N)/yr 

Year 


1955 

1960 

1965 

1970 

1975 

Commercial  nitrogen  fixation 

7.4 

11.9 

19.1 

30.7 

49.5 

Combustion 

10.9 

14.0 

17.8 

22.7 

29.0 

Total* 

18.3 

25.9 

36.9 

53.4 

78.5 

★ 

Effect  of  soybeans  and  other  leguminous  crops  (see  Section  F.2.3)  is 
not  included  here. 
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3.  OTHER  FACTORS 

3.1  SOLAR  ULTRAVIOLET  RADIATION  AS  A FUNCTION  OF  SOLAR  ACTIVITY 

The  solar  ultraviolet  irradiance  determines  the  overall 
level  of  stratospheric  ozone,  principally  by  the  photodissocia- 
tion processes  which  lead  to  the  formation  of  ozone  in  the  firs 
place,  and  also  by  the  local  heating  resulting  from  the  absorp- 
tion of  solar  UV  and  the  changes  in  temperature,  photochemistry 
and  dynamics  resulting  from  this  heating. 

Different  measurements  of  solar  UV  radiation  in  the  120  to 
300/400  nm  wavelength  range  have  been  reviewed  by  Heath  and 
Thekaekara  (1977),  by  Delaboudiniere  et  al.  (1977),  and  by 
Simon  (1978),  all  of  whom  infer  a significant  variation  of  the 
solar  200  to  300  nm  irradiance  (20  to  50  percent  variation) 
with  the  11-year  solar  sunspot  cycle,  with  more  radiation  at 
solar  maximum  than  at  solar  minimum.  This  is  not  inconsistent 
with  comparable  or  larger  effects — factor  two  at  170  to  200  nm 
and  at  Lyman-  a [see  also  Vidal-Mad jar  (1977)]  at  shorter 
wavelengths . 

The  results  of  Heath  and  Thekaekara  (1977)  are  shown  in 
Fig.  3*1,  were  redrawn  by  the  1977  NASA  Committee  [see  Hudson 
(1977)]  as  shown  in  Fig.  3.2,  and  have  been  incorporated  in 
several  model  calculations  of  stratospheric  ozone,  see  Brasseur 
and  Simon  (1977);  Nealy,  Callis,  and  Natarajan  (1977);  and  in 

i 

These  are  not  treated  at  the  level  of  detail  of  the  items  in 
Chapter  2,  but  they  are  potentially  of  considerable  importance 
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Source:  Heath  and  Thekaekara,  1977,  Fig.  8 

FIGURE  3-1.  Variations  in  Solar  Spectral  Irradiance  Aoparently  Related  to 
the  11-year  Sunspot  Cycle,  Based  on  Observations  from  1964  to 
1975. 


Source:  Hudson,  1977,  Fig.  17 

FIGURE  3-2.  Ratio  of  Solar  Flux  Measured  Near  Solar-Cycle  Minimum  to  that 
Measured  Near  Solar-Cycle  Maximum  Versus  Wavelength. 
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particular  Callis  and  Nealy  (1978).  These  modelers  typically 
find  an  increase  in  total  ozone  of  8 percent  as  one  goes  from 
solar  minimum  to  solar  maximum.  This  is  a very  large  effect, 
which  could  overwhelm  all  others  during  the  1955  to  1975  time 
period . 

The  conclusions  of  Heath  and  Thekaekara  (1977)  have  been 
criticized  by  a number  of  workers,  see  White  (1977),  pp . 21° 
and  325;  also  Delaboudiniere  et  al.  (1977);  while  Simon  (1978) 
estimates  a somewhat  smaller  effect;  indeed  Heath  (private 
communication,  February  1978)  considers  that  while  the  effect 
is  real,  its  magnitude  is  not  yet  established. 

I suggest  that  the  correlation  of  Heath  and  Thekaekara 
(1977)  should  be  used  but  with  the  caveat  that  the  existing 
length  of  record  (from  1964  on)  is  inadequate  to  establish  such 
a variation:  two  cycles  (22  years)  would  be  needed  as  a mini- 
mum. The  Heath-Thekaekara  estimate  of  Fig.  3-1  is  the  largest 
variation  currently  proposed;  if  it  withstands  criticism,  UV 
variability  may  well  be  the  single  most  important  contributor 
to  changes  in  stratospheric  ozone  during  the  1955  to  1975  time 
period,  but  at  present  (June  1978)  there  is  not  sufficient 
evidence  to  provide  a quantitative  estimate  of  the  effect. 

3.2  STRATOSPHERIC  WATER  VAPOR 

The  stratospheric  odd-oxygen  concentration  is  reduced  to  a 


significant  extent  by  the  HO  cycle: 

03  + OH  H02  + 02 

(a) 

H0o  + 0 -+  OH  + 02 

(b) 

(3.1) 

Net  : 0 + 0,  ->  2 02 

(c) 
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Also,  the  HO  species  react  with  NO  and  C1X  species  In  a variety 
A.  X 

of  ways.  The  HO  comes  from  stratospheric  water,  and  0(  D); 

1 x 

the  0 ( D)  is  produced  by  photodissociation  of  0^  (and  of  0^ ) 
while  the  water  comes  mainly  by  transport  from  the  troposphere, 
but  partly  by  oxidation  in  the  stratosphere  of  methane  ( CH^ ) 
transported  from  the  ground. 

Our  knowledge  of  the  burden  of  stratospheric  water  vapor  is 
very  poor— see  Kuhn  et  al.  (1975)  and  Harries  (1976).  There  are 
data  in  England  and  in  Washington,  D.  C.,  extending  over  almost 
20  years,  which  give  some  indication  that  the  overall  level  of 
stratospheric  water  vapor  may  have  been  increasing  over  the  1950 
to  1970  time  period--see,  e.g.,  Pigs.  5 and  6 of  Harries  (1976). 
However,  the  stratosphere  is  so  very  dry  compared  with  the  tropo- 
sphere that  measurements  are  very  difficult,  and  there  is  a con- 
siderable spread  of  data  given  by  different  methods.  Indeed, 
Harries  does  not  think  that  an  increasing  trend  has  been  estab- 
lished. 

There  are  other  analyses  which  support  the  suggestion  of  an 
increase  in  stratospheric  water  vapor  during  this  time  period. 
Thus,  Angell  and  Korshover  (197*0  find  some  evidence  for  an  in- 
crease in  the  tropopause  pressure  and  temperature  in  the  i960  to 
1970  decade,  particularly  in  tropical  latitudes,  which  may  per- 
haps be  related  to  an  increase  in  stratospheric  water  vapor  in 
terms  of  the  still-current  model,  due  originally  to  Brewer  (19^9) » 
that  the  concentration  of  water  vapor  In  the  stratosphere  is 
limited  by  the  "cold  trap"  in  the  tropical  torpopuase.  Also, 
systematic  reports  of  noctilucent  clouds  (see  Newell  et  al., 

197*4,  p.  265)  indicate  a relative  maximum  in  cloud  sightings  in 
1967  to  1968,  which  may  imply  a corresponding  maximum  in  strato- 
spheric water  vapor  (or  a minimum  in  stratopause/mesophere 
temperature ! ) . 
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Dverall,  the  current  status  of  stratospheric  water  vapor 
data  bears  some  resemblance  to  the  ozone  data  of  the  early  1950s, 
when  there  had  been  some  decades  of  observations  at  a few  closely 
spaced  locations  in  Western  Europe  (mainly  Oxford,  Arosa,  and 
Tromso),  but  simply  not  enough  global  coverage  to  permit  any 
inference  of  long-term  worldwide  trends. 

Further,  we  do  not  yet  understand  the  budget,  sources,  and 

sinks  of  stratospheric  water  vapor  very  well,  particularly  as  it 
relates  to  general  stratospheric  transport.  The  current  inter- 
pretation is  based  on  transport  both  into  and  out  of  the  strat- 
osphere by  the  Hadley  circulation,  with  considerable  disagreement 
about  what  happens  at  higher  latitudes  (the  Hadley  circulation 
is  restricted  to  latitudes  less  than  30  to  35°,  north  or  south). 
Thus,  Kuhn  et  al.  (1975)  and  Ellsaesser  (197*0  present  quite 
different  models. 

Long-term  global  measurements  of  stratospheric  water  vapor- 
ana  relative  humidity  in  the  lower  stratosphere  and  upper  tropo- 
sphere are  irrently  needed.  The  best  platform  would  presumably 
be  a sat<  Hite,  ‘alibrated  against  balloon-borne  (or  other 
ground-based)  data,  in  a limb  scan  mode.  For  budget  purposes, 
the  concent  rat  ion  lata  clearly  have  to  be  supplemental  with 
velocities  of  upwelling  and  subsidence  as  a function  of  place 
and  time. 

3.3  VARIATIONS  IN  STRATOSPHERIC  TEMPERATURE 

There  is  an  interaction  between  solar  irradiance,  ozone 
concentration,  temperature,  and  dynamics  which  is  now  being 
characterized  by  two-dimensional  models  [see,  e.g.,  Louis, 

London,  and  Danielsen  (197*0].  Further,  the  rate  constants  for 
most  of  the  chemical  and  photochemical  processes  that  determine 
stratospheric  ozone  depend  on  temperature,  and  thus  changes  in 
stratospheric  temperature — by  affecting  the  rates  of  different 
reactions  dif f erently--may  be  expected  to  change  the  total  ozone 
: ;ent ration  [see,  e.g.,  Boughner  (1978)  and  Luther  and  Duewer 
(1978 ) ] . 
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En  addition  to  <1 1 urna  1 vari  tl  ions , the  stratospheric  t< -m- 
perature  shows  also  semi-annual,  mnual,  biennial,  and  1 
term  variations--see  Newell  et  al.  (197*0.  A number  of  workers 

iriv(  looked  at  longer  term  variations,  in  particular  in  its  reln- 

# 

t : on  t the  (tropical)  quasi-biennial  oscillation  (QBQ)  and  the 
solar  cycle.  Tims,  diurnal  and  seasonal  variations  are  well 
established,  as  is  the  QBO  at  low  latitudes  (£  30°)  (cf.  e.g., 
Angell  and  Korshover  (1973 > 1976,  1978).  Figure  3-3  shows  some 
of  the  correlations  between  total  ozone,  tropical  tropopause 
temperature,  and  stratospheric  water  vapor  measured  at  Washington 
D.  C.  Correlations  between  total  ozone  and  temperature  at 
different  temperatures  in  the  stratosphere  are  found  at  a number 
of  mid-latitude  stations  on  a scale  of  1 to  3 years  [cf.  Seki- 
hara  et  al . (1977)].  This  may  just  be  a representation  of  the 
QBO.  There  are  also  longer  term  effects,  such  as  an  apparent 
breakdown  of  the  QBO  associated  with  the  volcanic  eruption  of 
Mt.  Agung  in  1963,  which  is  discussed  further  in  Section  3-5, 
and  a significant  long-term  cooling  in  the  Northern  Hemisphere 
at  21  to  24  km  altitudes  of  almost  1°  K between  1957  and  1972 
(cf.  Zorefos  and  Mantis,  1977)  • 

It  is  not  apparent  that  these  effects  can  yet  be  incor- 
porated in  models  of  stratospheric  photochemistry  in  a signifi- 
cant way . 

3.4  DYNAMICAL  EFFECTS 

Two-dimensional  photochemical  models  incorporate  atmospheric 
transport  in  a semi-phenomenological  way,  typically  separating 
large-scale  mean  motions  from  smaller  scale  (sub-grid)  eddy  trans 
port,  see  Danielsen  and  Louis  (1977).  The  parameterization  gives 
seasonal  averages  and  provides  a reasonably  satisfactory  repre- 
sentation of  ozone,  which  is  created  and  mostly  destroyed  in  the 
_ 

The  QBO  is  apparently  due  to  a tidal  forcing  by  the  planet  Mars 
with  a period  of  25.6  months  (see,  e.g.,  Fairbridge,  1978). 
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Source:  Angel!  and  Korshover,  1978 

FIGURE  3-3.  Comparison  of  Time  Trends  in  Given  Parameters  Where  a 1-2-1 
Smoothing  Has  Been  Applied  to  Successive  Annual  Values.  The 
16-km  temperature  at  the  equator  approximates  the  tropopause 
temperature  there.  The  data  are  expressed  as  either  percent- 
age deviations  or  temperature  deviations  (degrees  Celsius) 
from  the  mean.  The  arrow  at  the  top  indicates  an  intense 
polar  cap  absorption  (PCA)  event. 
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stratosphere,  as  well  as  • " some  radioactive  tracers  assoei 
with  atmospheric  nuclear  explosions.  It  can  also  be  used  to 
represent  water  vapor,  which  is  brought  into  the  si; rat ■ \u  h-r- 
by  the  Hadley  cell  upwelling  through  the  tropical  tropopaust , 
i.e. , by  large-scale  motions , sup?  Lemented  t some  exl  nl  by 
smaller-scale  motions  such  a-,  cumulus  • nvect  ion.  The  princij 
return  flow  is  near  30°  North  and  South  at  the  trot opauce  f : •- 
ing  associated  with  the  subtropical  jets  and  the  downdraft  in 
both  Hadley  and  Ferrell  cells.  In  the  stratosphere  th-'-r  i. 
a significant  horizontal  transport,  supplemented  to  sons,  e>; 
by  upwelling  near  the  subpolar  jet  stream  associated  with  \ be 
updraft  near  60'  North  and  South  associated  with  both  Ferrell 
and  polar  cells  of  circulation.  This  is  strong  in  winter,  : r 
often  nondetectable  in  summer,  with  some  downward  transrn  nt 
near  the  poles.  A good  review  of  both  the  physics  and  para- 
meterization is  given  by  Danielsen  and  Louis  (1977),  who,  ■ •- 

ever,  do  not  say  much  about  the  problems  of  water  vapor  • rain;:  . 

There  is  some  evidence  that  this  circulation  applies  t 
tracers,  e.g. , Vedder  et  al.  (1978),  who  find  higher  so rat  u ri 
concentrations  of  halocarbcns  and  N.,0  near  the  equator  .Fun  ■ 
higher  latitudes  consistent  with  the  picture  of  Hadley  cell  up- 
welling.  There  is  a strong  seasonal  variation:  thus,  mi  i-lal  11  id 
stratosphere-to-troposphere  transfer  is  much  more  visors.  ■ in 
winter  than  in  summer  [see,  e.g.,  Reiter  (1975),  and  work  qs.tei 
therein] . 

This  picture  describes  an  average  of  motions  on  a variety 
time  (and  space)  scales  from  very  short  to  rather  long,  whose 
relative  contribution  is  not  yet  understood.  Thus,  for  Instance, 
stratospheric  sudden  warmings  (SSWs)  [see  Quiroz  et  al . (197;  ), 
Mclnturff  (1978),  and  R.  S.  Quiroz,  private  communication]  are 
very  dramatic  phenomena,  occurring  in  variable  strength  dur: 
many  winters,  and  of  different  strengths  in  the  Northern  and 
Southern  Hemispheres,  leading  to  dramatic  enhancements  in  tr  r 
sphere/stratosphere  mixing.  In  particular,  there  is  some 
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evidence  that  in  the  years  of  major  stratospheric  warmings 
(which  are  characterized  by  reversals  in  zonal  circulation  of 
the  polar  vortex  at  10  mbar)  the  normal  springtime  maximum  in 
stratospheric  ozone  occurs  significantly  earlier  than  usual  — 
cf.  Zullig  (1973),  P-  133  in  Mclnturff  (1978),  Ellsaesser  (1976). 
Table  3-1  gives  a representative  listing  of  Northern  Hemisphere 
SSWs;  there  is  much  less  evidence  on  Southern  Hemisphere  SSWs, 
but  one  (or  more)  minor  SSWs  appear  to  occur  during  each  Northern 
and  Southern  Hemisphere  winter  (i.e.,  January/February  and 
July/August,  respectively).  However,  to  date,  no  major  warmings 
have  been  reported  in  the  Southern  Hemisphere.  In  any  case,  to 
date,  SSWs  have  not  been  incorporated  in  models  of  stratospheric 
ozone . 

It  must  also  be  stressed  that  the  parameterization  of  motions 
as  used  in  the  two-dimensional  models,  with  a mean  wind  u (hori- 
zontal, N-S),  and  u (vertical)  supplemented  by  subscale  eddy 

transport  coefficients  K , K , K is  not  necessarily  adequate; 

yy  yz  zz  J 

certainly,  smaller  scale  transport  is  not  always  downgradient , 
so  that  the  parameterization  in  terms  of  eddy  transports  is  open 
to  serious  question. 

In  conclusion,  while  the  standard  parameterization  can  de- 
scribe the  long-term  averaged  overall  motion  of  ozone  and  water 
in  the  stratosphere,  not  all  sets  of  values  used  for  the  numeri- 
cal coefficients  K , etc.,  are  equally  satisfactory.  Further, 

y y 

we  know  that  there  have  been  changes  in  stratospheric  motions 
in  different  years,  associated  with  such  events  as  major  volcanic 
eruptions  (e.g.,  Agung  in  1963),  the  quasi-biennial  oscillations 
and  SSWs  [cf.  Fig.  2-1  and  also  Newell  et  al.  (197*0,  Chapter  10]. 
The  effects  of  these  phenomena  on  stratospheric  transport  and 
hence  on  the  action  of  the  various  time-dependent  injectants  on 
stratospheric  ozone  cannot  yet  be  modeled  adequately. 
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TABLE  3-1.  A LIST  OF  NORTHERN  HEMISPHERE 
STRATOSPHERIC  SUDDEN  WARMINGS  (SSWs) 


3.5  SOME  EFFECTS  OF  MAJOR  VOLCANIC  ERUPTIONS 


The  large  eruption  of  Mt . Agung  in  March  1963  occurred  at 
a time  when  there  was  available  a significant  global  network 
capable  of  providing  estimates  of  upper-atmosphere  as  well  as 
of  ground-level  effects.  Thus,  we  have  very  much  more  infor- 
mation on  the  Agung  eruption  than  on  the  significantly  larger 
ones  of  Krakatoa  ( 1 8 8 3 ) and  of  Tambora  (I8l6)  for  which,  in 
addition  to  information  on  physical  destruction  and  on  enhanced 
"dust  veils"  (Lamb,  1970),  we  only  have  data  on  decreases  in 
surface  temperature  [e.g.,  Newell  et  al . (1974),  P-  2 h 6 ] . 

Thus,  we  know  that  the  volcanic  eruption  of  Agung  not  only 
injected  large  amounts  of  material  into  the  stratosphere  (se<_- 
Section  2. A and  Appendix  D),  but  also  led  to  an  increase  in 
temperature  in  the  tropical  and  subtropical  lower  stratosphere- 
(presumably  due  to  absorption  of  solar  radiation  by  aerosols 
volcanic  origin)  and  a probably  related  decrease  in  tropospheric 
and  ground-level  temperatures  [see  Newell  et  al.  (197*0,  pt  . 
264-267]. 

Since  water  is  thought  to  enter  the  stratosphere  principally 
through  the  "cold  trap"  at  the  tropical  tropopause,  such  ar. 
crease  in  tropical  tropopause  temperature  would  be  expeotc  i *. 
lead  to  an  increase  in  stratospheric  water  vapor  content,  and 
Newell  (op.  cit.)  suggests  that  the  maximum  in  noctilucent  cl 'Ud 
observations  in  1967  and  1968  may  thus  be  due  to  the  Agung 
eruption,  which  raised  the  temperature  of  the  tropical  tropo- 
pause and  thus  permitted  more  water  to  enter  the  stratosphere. 
Presumably,  the  4-year  delay  between  the  eruption  and  the  effect 
on  noctilucent  clouds  is  related  to  the  time  required  for  water 
to  be  transported  from  the  tropical  tropopause  to  the  polar 
stratopause . 
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3.6  TROPOSPHERIC  POLLUTANTS  AND  CHEMISTRY  EFFECTS 


As  our  understanding  of  the  complex  photochemical  inter- 
actions of  trace  species  in  the  atmosphere  has  progressed, 
additional  processes  by  which  stratospheric  ozone  may  be  affected 
have  been  noted.  Thus,  it  has  been  suggested  by  various  workers 
in  this  field  that  any  tropospher ic  pollutant  which  affects  the 
OH  concentration  can  result  in  changes  to  stratospheric  ozone 
because  reaction  with  OH  is  the  principal  tropospheric  destruc- 
tion mechanism  for  species  such  as  CHj,  and  CH,C1.  If  OH  is 
reduced  (e.g.,  by  reaction  with  CO),  or  enhanced  e.g.,  by  the 
reaction 


HOp  + NO  - NOp  + OH,  (3.2) 

the  tropospheric  lifetimes  of  these  compounds  are  altered  and 
their  fluxes  into  the  stratosphere  change.  Increases  in  CH^ 
lead  to  increased  Ho0  at  altitudes  where  HOx-mediated  ozone 
depletion  is  important,  while  increases  in  CH^Cl  increase  the 
C1X  effects. 

"These  preliminary  suggestions  ought  to  be  considered  in 
detailed  modeling  studies.  In  particular,  if  ground-level 
sources  of  halocarbons,  hydrocarbons,  and  fixed  nitrogen  are 
taken  into  account,  it  is  important  that  tropospheric  as  well  as 
stratospheric  photochemistry  should  be  included  in  the  models. 

Note  that  the  atmospheric  carbon  dioxide  level  has  been 
increasing  during  the  past  100  years  or  so,  principally  as  a 
result  of  fossil  fuel  combustion  resulting  from  the  industrial 
revolution.  Presumably,  this  has  led  to  a small  increase  in 
stratospheric  ozone  through  changes  in  stratospheric  temperature, 
see  Boughner  (1978)  and  Groves  et  al . (1978). 
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AIRCRAFT  SOURCES  (Of  NOy,  H^O,  ETC)  AND  ROCKET  SOURCES 


A listing  of  aircraft  so  j i- ■ • < f <p  1075  ! ;;  ; 

3.2  as  a fun  ifcion  of  alt  itud<  and  lat.  itude . 

A . D . Little  (1976 ) , corrected  a;  indicat  1 f 
( 1977 ) » P • 2-30  ff  for  the  fad  that  v 

X 

engines  are  roughly  equal  t-  those  ,2200  • .•  ; 

rection  amounts  to  approximately  14  percent.  /-  ••  1 

operations  have  been  increasing  very  rapidly  . i:- 
aircraft  came  into  widespread  use  in  the  1930:;. 
figures  on  passenger  travel  per  year  (Source  : Lrr  <•  . # 

1974 , p.  977;  ibid.  October  1976 , p.  927;  ibid.  ;1  l r • ; . 

p . 957 ) show  a mean  annual  increase  of  ),4  percei  • 

pounded.  While  passenger  travel  on  scheduled  a ' 1 : 1 ■ r and  ' 

NO  emissions  are  not  strictly  proportional,  and  indeed  the  fl 
altitude  has  been  increasing  s >mewhat  sines  th<  nid-19  s, 
reasonable  scaling  factors  for  reducing  the  1975  in;e;aio::  ; it 
of  Table  3.2  to  those  for  earlier  y<  ars  a ;lven  in  rat  .3. 

Regarding  the  effect  of  a steady-s  it<  airpl  me  inject!  . , 
note  that  scaling  from  the  calculations  of  Hidalgo  and  0 rut sen 
( 1977  ) and  o f Widh"  : , Jlatt , Kramer  t<  take  ini 

account  changes  in  fuel  flow  ral  uni  i,  and  al:  a revisi 

in  the  rate  constant  for  tl  re  stioi  • • . ! ) by  a factor  40 
[Howard  and  Evenson  (1977)  and  other  workers]  one  would  expect 
■ >m<  -'Zone  r<  list!  ibution  and,  in  pari  ; :ular,  a net  ozone  in- 
crease of  order  l.  percent  1 ) result  fr  m a si  ;ady-stat<  air 
plane  injection  at  the  1975  rate.  Evidently , the  effect  f 
actual  aircraft  operations  during  the  1955  to  1975  time  period 
is  much  smaller  than  this  because  of  the  very  rapid  increase  i 
injections  shown  in  Table  3.3. 

Rocket  sources  were  very  small  during  the  1955  to  1975  to. 
period:  current  estimates  indicate  that  60  space  shuttle  launches 
per  year — whicti  would  represent  a far  greater  activity  than  any 
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TABLE  3-2.  (Continued) 
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TABLt  3-2.  (Continued) 
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dui  Lng  the  19S5  to  19  75  time  frame — would  produce  a steady— state 
ozone  depletion  of  order  0.2  percent  [see  NASA  (1977)  addendum, 
also  Wuebbles  (1977)]. 


TABLE  3-3.  AIRCRAFT  EMISSION  RATES  AS  A FUNCTION  OF  TIME, 
NORMALIZED  TO  THE  1975  EMISSION  RATES  OF  TABLE  3.2. 


,'gar  1955  1960  1965  1970  1975 

Emission  Factor  0.17  0.27  0.41  0.64  1.00 
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4.  DISCUSSION  OF  RELATIVE  EFFECTS  AND  UNCERTAINTIES 

Our  understanding  of  the  relative  importance  of  NO^ , HO^  , 

and  C1X  injections  to  ozone  has  been  changing  rapidly  since 

1974.  Thus,  it  is  now  thought  that  C1X  injections  are  twice 

as  effective  in  reducing  ozone  (relative  to  NO  injections) 

x 

as  was  the  feeling  in  1976 — see,  e.g.  National  Academy  of 
Sciences  (1977).  These  changes  have  resulted  from  revised 
estimates  of  reaction  rates,  the  most  important  of  which  is  for 
HO,,  + NO  going  to  OH  and  NO,,.  The  first  direct  measurement 
of  this  reaction  rate  was  made  by  Howard  and  Evenson  (1977). 

See  Hampson  and  Garvin  (1973),  p.  33,  for  a listing  of  other 
measurements . 

Predictions  of  ozone  reduction  due  to  the  nuclear  test 
series  of  1961  and  1962  have  been  changing;  thus  Chang  et  al. 
(1977),  using  a one-dimensional  model,  predicted  an  8 percent 
ozone  reduction  using  1974  chemistry  as  against  a 4 percent  re- 
duction using  1977  chemistry.  Their  estimates  of  injection 
height  are  higher  than  those  adopted  here;  our  estimate  of  in- 
jection heights  should  reduce  the  expected  ozone  depletion  still 
further.  Explicitly,  in  1974  it  was  thought  that  NO  injections 

X 

reduce  total  ozone  above  13  km,  but  increase  it  (by  "smog"  re- 
actions) below  13  km  (cf.  Johnston,  1975).  By  late  1977,  this 
cross-over  height  had  gone  up  to  17  km  (cf.  H.  S.  Johnston, 
private  communication)  which,  together  with  couplings  between 
the  HO  , NO  , and  C1X  cycles,  contribute  to  the  Chang  et  al. 
(1977)  result. 
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Overall,  It  now  appears  that  tne  predicted  oz'ne  reduction 
arising  from  a given  C1X  injection  is  roughly  twice  as  large  as 
was  expected  In  1976 — set,  e.  g.  , NASA,  1977  Addendum.  Thus, 
the  effects  of  halocarbons  and  of  volcanic  eruptions  [especially 
if  the  latter  inject  as  much  HC1  into  the  stratosphere  as  is 
suggested  by  Cadle  et  al.  (1977)1  may  well  b<  as  imj  rtanl  as 
were  the  1961  to  1962  atmospheric  nuclear  ext losions. 

It  is  probably  impossible  to  make  a reliable  quantitative 
estimate  of  the  various  uncertainties  in  the  catalog,  but  some 
indications  can  be  provided.  As  far  as  nuclear  explosions  go, 
there  are  uncertainties  in  NO  product ion/Mt , from  (0.4  to  1.5) 

"5  0 ^ 

x 10J“  molecules/Mt , in  injection  height  (see  Fig.  B-5  of  Appen- 
dix B),  and  in  detrainment  (see  Section  B.3-3  of  Appendix  B). 
Regarding  halocarbons,  production  and  injection  rates  are  pre- 
sumably quite  accurate  (+  10  percent  for  production,  somewhat 
worse  for  atmospheric  injections),  although  different  estimates 

for  current  carbon  tetrachloride  production  differ  by  2C  to  50 

a * 

percent;  the  main  uncertainty  here  is  provided  by  tropospheric 
transport  and  photochemistry,  in  particular,  the  concentration 
of  OH  and/or  other  reactive  species  which  destroy  these  materials 
in  the  troposphere.  The  injection  estimates  for  volcanic  erup- 
tions listed  in  Tables  2-6  and  2-7  are  evidently  quite  poor.  As 
far  as  cosmic  ray  effects  go,  data  on  the  1959  and  I960  SPEs  are 
quite  limited,  but  since  the  effects  of  SPEs  (or  of  OCRs,  for 
that  matter)  are  not  very  large,  this  uncertainty  is  probably 
not  very  important.  Finally,  regarding  nitrogen  fixation,  there 


'A  very  recent  note  by  Stief  et  al.  (1978)  on  the  reaction  of 
stratospheric  formaldehyde  with  atomic  chlorine  to  produce 
the  relatively  inert  HC1  suggests  approximately  a 10  percent 
reduction  in  sensitivity  of  stratospheric  ozone  to  chlorine 
perturbations . 

< 

This  uncertainty  is  not  very  important  since  during  the  1955 
to  1975  time  period  the  principal  stratospheric  effect  of  CCl^ 
comes  from  the  large  pre-1950  injections  (which  are  relatively 
well  known),  rather  than  from  the  small  post-1950  injections 
where  there  is  a major  discrepancy. 
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are  significant  uncertainties  in  the  nitrogen  fixation  by  com- 
bustion, in  the  rate  and  product  of  denitrification,  and,  in 
particular,  in  the  fate  and  budget  of  Nn0  in  the  troposphere . 

Among  all  these  factors,  I think  that  the  uncertainties  in 
volcanic  injections  may  well  be  of  greatest  significance  in 
the  present  context. 

Compared  with  these  injections,  the  different  factors  dis- 
cussed in  Chapter  3 are  far  more  uncertain  and  ill-defined,  and 
limitations  in  our  overall  understanding  may  well  affect  overall 
predictions  to  a significant  extent.  In  particular,  if  the  es- 
timates of  solar  UV  variation  shown  in  Fig.  3-2  apply,  and  if  the 
corresponding  predictions  of  an  8 percent  ozone  variation  be- 
tween solar  minimum  and  maximum  quoted  in  Section  3.1  are  appro- 
priate, then  this  is  probably  the  biggest  single  time-varying 

1 

effect  during  the  1955  to  1975  time  frame. 
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APPENDIX  A 


STRATOSPHERIC  BURDENS  OF  NO  , HO  , C1X,  ETC. 

x x 

For  orientation,  it  is  desirable  to  have  estimates  of  the 
total  stratospheric  burden  of  the  NO  HO  , and  C1X  groups  of 

X X 

molecules.  The  present  estimates  are  very  poor;  they  are  based 
primarily  on  available  experimental  data,  supplemented  by  cur- 
rent theoretical  computations  of  Crutzen  et  al . (1978)  and  Wofsy 

* 

(1978).  It  should  be  noted  that  the  stratospheric  mean  life- 
time of  HO  is  significantly  shorter  than  those  of  NO  and  C1X, 
so  the  relative  importance  of  HO^,  NO^ , and  C1X  perturbations  is 
not  simply  proportional  to  the  respective  number  of  molecules 
inj ected . 

For  some  of  the  more  reactive  species  (OH,  HOp  , Cl,  CIO) 

we  only  have  a daytime  profile  near  30°  N latitude  which  is 

taken  as  a global  mean.  Thus,  for  NO  and  C1X  species,  the 

global  burden  is  taken  as  this  integrated  profile  multiplied 

1 A 2 

by  the  surface  area  of  the  earth,  = 5.11  x 10  m . However, 

for  HO^  species  which  are  produced  by  sunlight  and  have  a mean 

lifetime  less  than  12  hours,  we  use  only  the  sunlit  fraction  of 

the  earth's  surface  area.  A,  = 1/2  A . 

d s 

A. 1 OXIDES  OF  NITROGEN 

NO  = NO  + N09,  NO  = NO  + UNO  (N_Or,  etc.,  are  sometimes 
x 2 y x 3 2 5 o’i 

included).  Estimates  for  NO  range  from  4.4  x 19 J molecules 

(Johnston,  1975,  p.  5-109),  to  8.5  x 103^  (after  Penndorf,  1977) 

or  15  X 1034  (COMESA,  1975,  p.  26).  We  suggest  (4-15)  x 1034 

as  a plausible  range. 

* 

Some  of  the  numbers  are  written  with  three  digits  to  minimize 
internal  round-off  errors:  this  does  not  imply  significance. 
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A.  2 OXIDES  OF  HYDROGEN 


HO  = OH  + HO ^ (HO  should  be  included).  For  OH,  we 
use  Hudson  (1977),  p.  154,  Fig.  28. 

50  km 

Aj  J"  [OH]  dz  = 1.1  x in^‘  molecules  (Is 

11  km 

70  km 

A^  f [OH]  dz  = 0.4  x 10 ^2  molecules  (It 

r , . (not  included  here) 

5 0 km 

Figure  la  from  Crutzen  et  al.  (1978)  gives  0.95  x 10j2 

82 

molecules,  and  Fig.  6b  from  Wofsy  (1978)  gives  1.5  x 10 
molecules,  both  for  the  stratosphere  alone,  i.e.,  11  to  50  km 
altitude . 


Measurements  of  H0-,  are  currently  being  made  by  at  least 
the  Max  Planck  Institute  and  the  University  of  Michigan  groups, 
but  I have  not  yet  seen  sufficient  data  to  make  an  estimate, 
thus  present  model  estimates: 

I n ..  in-1  / f ' t , r i rr 1 . - ~ u "l 


50  km 

/ 


[HO  ] dz  = 


tropopause 


3 x 10' 


6 x 10' 
15  x 10: 


( C FAI’ , Mi >n< igrap h 1 , 
p.  3-96,  representa- 
tive values) 

(Crutzen  ct  al.,  data 
from  Fig.  1) 

(Wofsy,  data  from 
Figure  6b) 


Those  experimental  numbers  that  I have  seen  support 

Crutzen ’s  rather  than  Wofsy’ s numbers.  Thus,  one  can  say  that 

current  estimates  for  the  HO  = OH  + H0„  burden  in  the  strato- 

x 32  2 

sphere  lie  in  the  range  7 to  16  x 10  molecules,  but  evidently 
the  reliability  of  these  figures  is  not  high.  In  fact,  both  the 
production  and  destruction  rates  of  HO^  radicals  in  the  strato- 
sphere are  very  high,  so  that  while  the  burden  of  HO  is  low, 

the  injection  of  additional  HO  will  affect  the  atmospheric 

x 

photochemistry  for  only  a relatively  short  time. 


A. 3 CHLORINE  RADICALS 


C1X  = Cl  + CIO  + CIO  N0o  + HC1 . HC1  appears  to  be  the  most 
abundant  species  [see  Sundararaman  (1976),  p.  12;  Hudson  (1977), 
p.  97],  for  which 

50  km 

A / [HC1  ] dz  = 1.07  x 103*4  molecules  (3) 

11  km 

Three  measurements  of  CIO  have  been  reported  by  Anderson  et 
al . (1977);  they  differ  by  about  a factor  of  10. 

50  km 

Ad  / [CIO]  dz  = 

11  km 


j 0 . 0 9 x 103  molecules 
I ( largest  value  ) 
j 0.02  x 10^4  molecules 
' (smallest  value  ) 


The  other  species,  CIO  HO^  and  Cl  have  not  been  measured  in 
the  stratosphere,  and  thus  we  are  forced  to  rely  on  model  calcu- 
lations . 


50  km 

Ad  J [CIO  N02]  dz 

tropopause 


14.8  x 103^  molecules 
(Crutzen,  Fig.  la) 

< 6 x 10^  molecules 

(Murcray  et  al . , 1973, 
upper-bound  estimates) 


A 


d 


50  km 


[Cl] 


tropopause 


dz 


10.0  x 10  molecules 
( Crutzen,  Fig.  la ) 

1 x lO^O  molecules 
(Anderson  et  al . 
smallest  value) 

5 x 10^0  molecules 
largest  value) 


(6) 


Thus,  we  conclude  that  the  stratospheric  burden  of  C1X  lies 

34 

in  the  range  of  (1.0-1. 2)  x 10  molecules,  with  a large  uncer- 

# 

tainty . 

* 

Fig.  11a  of^Crutzen  et  a[ . reports  the  burdens  of  HC1  and  CIO 
as  37  x 10J‘-  and  3 x 1032  molecules,  respectively,  so  that  their 
mean  burden  of  C1X  would  be  0.4  x lO^4  molecules;  Fig.  4 of 
Wofsy  (1978)  gives  0.7  x 1034  molecules. 
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A. 4 SOME  OTHER  QUANTITIES 


The  stratospheric  burden  of  water  is  listed  as  (1.6-2. 2)  x 
lO1^  kg  [cf.  Oliver,  1977,  pp.  4-9],  or  (5-7)  x 10^  molecules. 
Regarding  stratospheric  aerosols,  it  is  impossible  to  make  a 
meaningful  statement,  since  it  is  considered  that  occasional 
large  volcanic  eruptions  are  the  principal  source  of  this  mate- 
rial. Thus,  in  1972-1973,  4 to  5 years  after  the  last  signifi- 
cant explosive  eruption  (that  of  Fernandina  in  1968),  the  global 

O 

loading  was  of  order  10  kg,  but  the  eruption  of  Fuego  in  1974 
probably  put  (4-6)  x 10^  kg  of  aerosols  into  the  stratosphere. 
See  Appendix  D for  a more  detailed  discussion. 
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APPENDIX  B 


STRATOSPHERIC  INJECTION 
OF  NITROGEN  AND  OTHER 
BY  ATMOSPHERIC  NUCLEAR 


OF  OXIDES 

materials 

EXPLOSIONS 


STRATOSPHERIC  INJECTION  OF  OXIDES 
OF  NITROGEN  AND  OTHER  MATERIALS 
BY  ATMOSPHERIC  NUCLEAR  EXPLOSIONS 

B.l  INTRODUCTION 

Foley  and  Ruderman  (1972)  pointed  out  that  the  theory  of 
the  depletion  of  stratospheric  ozone  by  oxides  of  nitrogen  by 
the  reaction  chain 

03  + NO  N02  + 02  (1) 

N02  + 0 - NO  + 02  (2) 

0 + 03  - 2 02  (3) 

lends  itself  to  verification  by  studying  global  ozone  changes 
related  to  the  Soviet  atmospheric  nuclear  test  series  of  1961 
through  1962. 

Explicitly,  an  atmospheric  nuclear  explosion  produces 
nitric  oxide,  mainly  by  heating  air  above  2000°  K;  during  the 
cooling  to  ambient  temperature,  some  5000  to  10,000  tons  of  NO 
per  Mt  yield  remains,  and  rises  into  the  stratosphere  for  yields 
in  the  Mt  range  (see  Section  B.2). 

For  orientation  on  the  scale  of  the  bomb-induced  perturba- 
tion, the  largest  single  test  series  (in  terms  of  total  yield) 
was  the  Soviet  series  of  Fall  1962,  which  accounted  for  some  180 
Mt  of  total  yield.  In  terms  of  the  current  estimates  of  NO  pro- 
duction per  Mt  (see  Section  B.2,  especially  Table  B-l),  this 

34 

would  correspond  to  an  injection  of  (0.7  to  2.7)  x 10  NO  mole- 
cules in  the  stratosphere  over  a period  of  several  months 
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(5  August  to  25  December  1962).  This  injection  in  the  northern 
polar  cap  of  the  stratosphere  (injection  at  75°  N latitude)  may 
be  compared  with  the  global  loading  of  stratospheric  HO  (=  NO 

•3/j  y x 

+ HNOO,  estimates  of  which  range  from  ('4-15)  x 10  molecules 
(see  Appendix  A).  We  see  that  the  fractional  enhancement  of 
HO  or  NO  in  the  northern  polar  cap  region  of  the  stratosphere 

^ y 

was  very  large.  In  terms  of  the  models  of  stratospheric  ozone 
of  the  1973  to  1974  time  period,  which  emphasize  control  of  its 
concentration  by  NO  and  also  by  HO  (=  OH  + H0?),  a detectable 
reduction  in  stratospheric  ozone  might  have  been  expected.  The 
earlier  (1973  to  197*0  models  predicted  a maximum  hemispheric 
ozone  reduction  of  8 to  10  percent,  while  the  later  ones  predict 
*4  to  6 percent,  cf  Chang  et  al.  (1977).  The  data  analysis  for 
total  ozone  gives  95  percent  confidence  interval  bounds  of  a 3 
percent  variation  in  total  ozone  (J.  K.  Angell,  private  communi- 
cation), while  H.  S.  Johnston  (private  communication)  finds  that 
when  a 10  percent  pulse  injection  is  applied  to  the  ozone  data 
with  the  Angell-Korshover  (1973)  smoothing,  its  effect  disappears 
in  the  noise.  Thus,  one  can  best  summarize  the  situation  by 
remarking  that  while  there  is  no  real  disagreement  between  pre- 
dictions and  observations,  existing  atmospheric  data  on  nuclear 
explosions  do  not  provide  a statistically  significant  demonstra- 
tion of  the  catalytic  destruction  of  ozone  by  oxides  of  nitrogen. 

To  date,  however,  only  one-dimensional  model  calculations 
have  been  made;  two-dimensional  models,  with  lat itudinally  and 
seasonally  varying  dynamics,  may  indicate  space-time  regions 
where  the  effects  should  be  most  apparent,  so  that  a more  criti- 
cal comparison  with  the  data  could  be  undertaken,  bearing  in 
mind  that  after  1955  a ground-based  network  of  some  *45  to  90 
Dobson  spectrophotometers  and  M-83  filter  ozonometers  has  been 
operational,  supplemented  from  1970  to  1977  by  the  backscattered 
ultraviolet  (BUV)  instrument  on  the  Nimbus-IV  satellite. 
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The  present  discussion  describes  the  physics  of  HO  produc- 
tion by  atmospheric  nuclear  explosions  (Section  B.2),  the  alti- 
tude of  injection  as  a function  of  yield  (Section  B.3),  some 
other  chemical  factors  (Section  B.4)  and  a listing  of  all  atmos- 
pheric thermonuclear  explosions  as  a function  of  date,  location, 
and  yield  (Section  B.5).  The  result  of  Sections  B.2  to  B.4 
provides  a listing  of  bomb-produced  N 0 ^ injections  as  a function 
of  position  and  time  (Section  B.6,  especially  Table  B-12),  which 
is  intended  to  provide  input  data  for  two-dimensional  computer 
codes  of  stratospheric  photochemistry  and  dynamics. 

B.2  N0x  PRODUCTION  IN  A NUCLEAR  FIREBALL 

A nuclear  weapon  emits  a very  large  amount  of  energy  in  a 
very  short  time  (^10~“  sec);  1 Mt  = 10^  Kcal  = 4.2  x lO’*"' 
joule.  Most  of  the  close-in  action  takes  place  very  rapidly, 
thus  in  a time  order  of  1 sec,  a "fireball”  is  formed,  which  is 
a mass  of  heated  air  (^10^  tons/Mt  yield)  heated  to  a tempera- 
ture of  order  6000°  K.  At  temperatures  of  this  order,  the  oxygen 
in  air  is  fully  dissociated  and  the  nitrogen  is  partially  dis- 
sociated. The  total  mass  of  the  nuclear  weapon  is  evidently 
very  small  compared  to  the  mass  of  the  heated  air,  as  an  energy 
of  1 Mt  corresponds  to  the  complete  fission  of  56  kg  of  fission- 
able material,  while  for  fusion  the  mass  is  even  lower,  by  a 
factor  of  almost  three. 

"The  color  of  the  (radioactive)  cloud  is  initially  red  or 
reddish  brown,  due  to  the  presence  of  various  colored  compounds 
(oxides  of  nitrogen,  especially  NO.-,)  at  the  surface  of  the  fire- 
ball. These  result  from  the  chemical  interaction  of  nitrogen, 
oxygen,  and  water  vapor  in  the  air  at  the  existing  high  tempera- 
tures and  under  the  influence  of  the  nuclear  radiations." 
(Glasstone  and  Dolan,  1977,  #2.08,  p.  29). 
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Nitric  oxide  (NO)  is  produced  in  the  fireball  principally 
by  the  rapid  heating  of  air  and  subsequent  "freezing"  of  the 
chemical  equilibrium  as  the  fireball  cools,  either  by  radiation, 
expansion,  or  mixing  with  cold  air.  The  precise  quantity  of 
NO  formed  depends  on  the  details  of  the  chemistry  and  on  the 
mixing  processes:  thus  as  the  freezing  temperature  (assuming 
equilibrium  in  dry  air  at  0.1  atmosphere  pressure)  changes  from 
2000  to  1800°  K,  the  mole  fraction  of  NO  changes  from  0.76  to 
0.^5  percent,  corresponding  to  a change  of  NO  production  per  Mt 
of  yield  (or  10()  tons  of  heated  air)  of  (1.6  to  0.9)  x 10^ 
molecules,  or  (8  to  4.5)  kt . (See,  e.g.,  Gilmore,  1974.) 

A secondary  mechanism  for  HO  production  comes  from  the 
ionization  of  air  by  gamma  and  X-rays,  neutrons,  and  beta  par- 
ticles. This  yields  roughly  one  ion  pair  per  35  eV  of  energy, 
and  1.3  NO  molecule  per  ion  pair  (e.g..  Section  E-3  of  Appendix 
E),  leading  to  10  x 10J^  NO  molecules  or  50  kt  of  NO  per  Mt  of 
that  fraction  of  the  bomb's  energy  that  comes  out  as  ionizing 
radiation.  However,  only  a small  fraction  of  the  bomb's  energy 
is  deposited  as  ionizing  radiation,*  so  that  this  is  not  the 
primary  mechanism  for  NO  production. 

It  is  clear  that  the  estimates  of  NO  production  are  subject 

to  large  errors.  A comparison  of  the  best  current  estimates  is 

listed  in  Table  B-l,  and  the  conclusion  from  the  work  of  Gilmore 

(1974)  is  that  NO  production  per  Mt  of  yield  is  (0.4  to  1.5)  x 
32 

10  “ molecules,  or  2 to  8 kt . 

This  discussion  corresponds  to  the  production  of  NO  in  the 
first  few  seconds  after  the  explosion,  and  before  the  "fireball" 
has  begun  its  expansion  and  buoyant  rise  to  its  stabilization 

•* 

For  a fission  weapon,  about  15  percent  of  the  energy  comes  out 
as  (prompt  and  delayed)  ionizing  radiation,  while  for  a thermo- 
nuclear weapon,  significantly  less  ionization  is  produced, 
since  this  comes  largely  from  the  fission  fragments.  See 
Glasstone  and  Dolan  (1977),  pp.  6-8. 


B-6 


altitude,  which  typically  takes  s me  1 min  ‘or  a ' -yield 
plosion.  Next , we  shall  discuss  h p ■<  ■ cloud 

stabilization,  which  is  physically  distinct  from  the  Lnit  la 
phase  of  NO  production. 

TABLE  B-l . CALCULATED  NO  PRODUCTION  BY 
LOW-ALTITUDE  NUCLEAR  EXPLOSIONS3 *b 


NO  Produced  per  Mt 

of  Yield  Calculation  Includes 


Molecules 

kt 

Shock 

Late  Fireball 

Reference 

0.5  x 1032 

3 

X 

Zeldovich  arid  Raizer  ( 1 967) 

0.3-2 

1.5-10 

X 

Foley  and  Ruderman  (1972) 

0.17-1 

1-5 

X 

Johnston  et  al  . (1973) 

0.3-1 .5 

1.5-8 

X 

Foley  and  Ruderman  (1973) 

1-2.5 

5-12 

X 

Simons  and  Caledonia  (1973) 

1 

5 

X 

Goldsmith  et  al  . (1973) 

0.4-1. 5 

2-8 

X 

X 

Gilmore  (1974) 

0.6-1 .1 

3-5.5 

X 

X 

COMESA  (1975)  (Section  6.3, 

P.  543ff) 


aThe  following  are  the  principal  mechanisms  for  NO  production: 

Shock  heating  and  subsequent  quick ' cool ing  of  air:  3.5  kt  NO/Mt  yield. 

Fireball  heating  of^ptir  and  slower  cooling  by  radiation,  expansion, 
and  mixing  in  of  cold  air:  -1,5  to  +4  kt  NO/Mt  yield. 

Bombardment  of  air  by  ionizing  radiation  (betas,  gammas,  neutrons): 

< 0.5  kt  NO/Mt  yield  (this  is  generally  smaller  than  the  other  mecha- 
nisms). (See  Gilmore,  1974.) 

^It  has  been  suggested  that  not  all  the  NO  produced  actually  rises  to  the 
stabilization  height  of  the  fireball,  as  a consequence  of  detrainment  dur- 
ing the  cloud  rise.  This  is  discussed  in  Section  B.3,3. 

B . 3 CLOUD  RISE  AND  STABILIZATION  AS  A FUNCTION  OF  YIELD 

B.3.1  Air  Versus  Surface  Bursts 

If  the  height  of  burst,  H,  for  a weapon  of  yield  Y (Mt)  is 
greater  than 

= 870  Y13 ’ ^ (meters),  (4) 
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then  the  fireball  does  not  ''touch  the  ground"  and  we  speak  of  an 

air  burst ; if  the  height  of  detonation  is  very  much  less  than  H, , 

-1- 

the  fireball  interacts  strongly  with  the  ground  and  we  call  it  a 
surface  burst. 

Many,  but  not  all,  thermonuclear  explosions  have  been  con- 
ducted in  an  air-burst  mode,  partly  to  minimize  local  radioactive 
fallout.  For  the  present  application  we  shall  assume  the  follow- 
ing : 

1.  The  cloud  rise  height  as  a function  of  yield  is  inde- 
pendent of  whether  we  refer  to  a surface  burst  or  an 
air  burst , 

2.  The  NO^  production  per  Mt  of  yield  for  a surface  burst 
is  half  the  value/range  estimated  in  Section  B.2  (see 
Table  B-l).  This  choice  is  made  in  analogy  with  the 
estimate  of  Machta  et  al . (1963)  that  a surface  burst 
produces  half  the  C-14  per  Mt  of  an  air  burst . 

B . 3 . 2 Cloud  Rise 

Figure  B-l  is  a sketch  of  the  air  circulation  inside  the 
rising  fireball,  and  Fig.  B-2  shows  the  cloud  rise  as  a function 
of  time,  to  indicate  the  phenomenology  and  scale  of  the  rise. 
Figure  B-3  shows  the  apparent  surface  temperature  of  a fireball 
as  a function  of  time;  from  the  discussion  of  Section  B.2  we 
know  that  most  NO  formation  is  complete  after  the  temperature 

drops  below  2000°  K so  that  NO  production  is  complete  before 

x * 

any  significant  cloud  rise  shown  in  Fig.  B-2. 


Tf 

The  results  of  Fig.  B-3  apply  to  a 20-kt  explosion  rather  than 
to  a 1-Mt  explosion  for  which  the  time  scale  would  be  stretch- 
ed out  by  a factor  of  several,  say  2 to  10. 
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TIME  AFTER  EXPLOSION,  minutes 


Source: 

FIGURE 


Glasstone  and  Dolan  (1977),  Fig.  2.12  and  Table  2.12 

-2.  Height  of  Cloud  Top  and  Rise  Rate  as  a Function 
of  Time  After  a 1 -Mt  Low  Air  Burst 


Source:  Glasstone  and  Dolan,  1977,  Fig.  2.123,  p.  69 

FIGURE  B-3.  Variation  of  Apparent  Fireball  Surface  Temperature 
with  Time  in  a 20-kt  Explosion 

Note:  After  a 1-Mt  explosion,  the  general  behavior  is  the  same, 
but  the  time  scale  is  somewhat  longer. 

For  this  application,  the  primary  interest  is  in  the  sta- 
bilization altitude  of  the  nuclear  cloud,  since  this  indicates 

the  altitude  distribution  at  which  the  NO  produced  by  the  ex- 

X 

plosion  is  injected  in  the  atmosphere.  In  Fig.  B-A  we  show  the 

* 

stabilization  altitude  as  a function  of  yield,  as  debris; 

1 

The  fact  that  radioactive  tracers  are  carried  on  aerosols 
which  are  subject  to  sedimentation  (unlike  gases)  appears 
entirely  immaterial,  as  sedimentation  velocities  are  of  order 
100  m/month,  which  is  10-!i  times  the  cloud  rise  velocity. 
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TOTAL  Yl El Dj megatons 


FIGURE  B-4.  Cloud  Rise  Height  as  a Function  of  Yield  for  Different  Lati- 
tudes. U.S.  tests  were  conducted  at  equatorial  latitude  (2  to  17°  N), 
Chinese  at  mid-latitude  (40°  N),  and  Soviet  at  polar  latitude  (75°  N). 
Peterson's  (1970)  curves  are  based  on  U.S.  tests;  all  U.S.  bombs  had 
yields  below  15  Mt.  Thus,  Peterson's  equatorial  curves  for  yields  above 
15  Mt  as  well  as  his  polar  curves  are  purely  estimates.  Estimates  of 
Seitz  et  al.  (1968)  for  individual  U.S.  (solid)  and  Soviet  (dashed) 
explosions  are  shown  as  vertical  lines.  Dot-dash  curves  have  been  drawn 
by  eye  through  Seitz'  Soviet  lines  to  provide  an  alternative  to  Peterson's 
polar  case.  Note  the  significant  difference  above  Mt-yields.  Some  mid- 
latitude estimates  of  Telegadas  (1977)  are  also  shown,  based  or,  Chinese 
explosions.  Within  a given  cloud  the  altitude  distribution  of  N0X  shown 
in  Table  B-2  may  be  assumed.  Representative  ozone  profiles  for  equatorial 
and  polar  latitudes  are  also  shown. 


TABLE  B-2 . INITIAL  VERTICAL  DISTRIBUTION  OF  RADIOACTIVITY 
OR  NOx  ASSUMED  WITHIN  "MUSHROOM  CAP"  (see  Fig.  B-5) 


Layer,  Fraction  of 
"Mushroom  Cap" 
(from  base  to  top) 


1 l 

7 ’ 7 

2 3 

7 ' 7 

3 4 
7 ' 7 

4 5 

7 ’ 7 

5 6 

7 ■ 7 

6 7 

7 " 7 

Source:  Peterson,  1970 


Percentage  of  Activity 
Within  Layer 

1 

14 

25 

25 

1 5 

15 

5 


we  may  assume  that  the  radioactivity  is  a tracer  of  cloud  motion, 

and  that  the  NO  follows  the  cloud,  with  the  exception  of  some 
x 

possible  detrainment,  which  will  be  discussed  below.  Peterson’s 
analysis  hypothesizes  the  distribution  of  radioactivity  within 
the  stabilized  cloud  as  indicated  in  Table  B-2.  It  refers  to 
air  bursts  and  shows  some  dependence  on  latitude;  there  is  also 
some  dependence  on  meteorological  conditions,  especially  for 
kt-yield  devices,  where  stabilization  altitudes  for  different 
shots  of  the  same  yield  may  vary  by  as  much  as  a factor  of  two. 

Note  that  Peterson's  estimates  are  based  on  U.  S.  data  which 

* 

cannot  be  relied  upon  above  about  15  Mt  and  which  were  always 


The  largest  announced  U.  S.  device  tested  was  Bravo  (28  Peb  1954), 
a 15-Mt  surface  burst  fired  at  Bikini  (11°  N latitude).  By  con- 
trast, the  largest  Soviet  device  was  a 58-Mt  air  burst  (fired  at 
Novaya  Zemlya,  75°  N latitude,  on  30  Oct  1961). 
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:arried  out  (for  Large-yi<  ; i levi  jes  at  equatorial  latitudes; 
the  estimates  for  polar  latitudes  are  based  on  a correlation 
with  atmospheric  stability,  and  are  clearly  not  as  reliable  as 
the  est  i mates  for  equal  irial  latitudes.  Here  we  may  refer  to 
tiic-  data  of  Telegadas  (197-1 , 1976)  on  Chinese  bombs  (3  Mt , 40°  n 
latitude)  which  come  from  airplane  and  balloon  samples  taken 
mainly  over  North  America,  a large  fraction  of  the  globe  away  in 
latitude  so  that  the  assumption  of  zonal  average  flow  is  crucial. 

Here  I use  assignments  of  Seitz  et  al.  of  both  weapon  yield 
and  cloud  rise  height  for  the  Soviet  bombs  as  shown  in  Fig.  B-4  , 

and  thus  draw  the  alternative  curves  for  polar  cloud  height 

* 

shown  dot-dash,  which  give  significantly  lower  cloud  heights 
than  do  Peterson’s  polar  assignments.  It  is  clear  that  this 
alternative  assignment  is  not  necessarily  better  than  Peterson’s, 
only  different;  and  the  difference  gives  merely  a suggestion  of 
the  uncertainty  of  the  cloud  height  assignment. 

B . 3 . 3 Detrai nment? 

The  rising  and  expanding  nuclear  cloud  is  extremely  turbu- 
lent, see  e.g..  Fig.  B-l,  and  entrains  a large  quantity  of  air 
during  its  rise.  One  may  ask  whether  any  of  the  hot  air,  NO  , 

X 

and  debris  is  "detrained"  from  the  cloud  and  left  behind  while 
the  rest  of  the  cloud  rises  past  it.  Thus,  Simons  and  Caledonia 
(1973)  suggested  that  only  17  tc  30  percent  of  the  NO  produced 
actually  rises  to  the  cloud  stabilization  height,  while  COMESA 
(1975),  p.  555,  estimates  that  about  80  percent  of  the  NO  pro- 
duced rises  to  stabilization  height. 

During  the  U.  S.  "Dominic"  nuclear  test  series  of  1962,  the 
stem  left  behind  when  a nuclear  fireball  rises  was  tested  for 
radioactivity  by  flying  an  aircraft  through  the  (visible)  stem, 

* 

These  curves  have  been  drawn  in  by  eye,  i.e.,  there  is  no 
"objective"  basis  for  them. 
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and  it  was  established  in  "Project  Stemwinder"  that  less  than  1 
percent  of  the  total  radioactivity  is  found  in  the  stem  1 hour 
after  detonation,  see  Perber  (1965).  This  suggests  that  no  sig- 
nificant detrainment  occurs,  and  thus  we  shall  assume  that  all 

the  NO  produced  rises  to  the  stabilization  altitude  as  shown  in 
x 

Fig.  B-4  and  Table  B-2. 

B . 4 OTHER  CHEMICAL  FACTORS 

In  addition  to  oxides  of  nitrogen,  atmospheric  nuclear- 
explosions  may  also  deposit  water  and  ozone  in  the  stratosphere. 
These  effects  are  estimated  here. 

B . 4 . 1 Water 

The  stratosphere  is  very  dry  (3  to  4.5  ppmv  of  Ho0),  whereas 

the  troposphere,  especially  the  tropical  troposphere,  can  be  very 
3 4 

wet  (10  to  10  ppmv  of  Ho0).  Thus,  large  nuclear  explosions 
can  carry  rather  large  amounts  of  water  vapor  into  the  strato- 
sphere. I will  show  here  that  these  quantities  of  water  appear 
to  be  insufficient  to  make  any  significant  change  in  strato- 
spheric water  vapor  content,  except  locally,  at  very  early  times 
after  the  events  (approximately  1 to  10  days).  No  attempt  is 
made  to  give  detailed  estimates  of  the  water-vapor  addition  for 
each  event,  but  some  guidelines  are  provided  for  a more  detailed 
investigation . 

Prom  Eq . (4)  we  see  that  a 1-Mt  "low  air  burst"  will  be 
detonated  at  about  1-km  altitude.  Table  B-3  shows  representa- 
tive temperature  and  moisture  contents  at  this  altitude  (assum- 
ing 50  percent  relative  humidity)  for  tropical  summer  and  polar 
fall  conditions  which  correspond,  respectively,  to  U.S.  and 
U.S.S.R.  test  conditions.  Data  are  also  shown  for  4-km  altitude, 
since  the  largest  Soviet  bomb  was  detonated  at  about  that  al- 
titude. The  final  column  in  Table  B-3  gives  plausible  estimates 
for  the  number  of  water  molecules  injected  into  the  stratosphere 
per  megaton  of  yield,  assuming  that  the  total  mass  of  the  fireball 
is  10  t/Mt  as  before,  and  that  initial  test  conditions  are 
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representative  of  the  total  entrained  air  mass.  Note  that  trop- 
ical explosions  carry  much  larger  quantities  of  water  into  the 

stratosphere  than  do  polar  explosions.  Note  also  that  the  total 

12 

water  content  of  the  stratosphere  is  roughly  2 x 10  kg,  7 x 
37 

10  molecules  (e.g..  Appendix  A).  We  may  approximate  the  water 

injected  by  the  1961  to  1962  series  of  atmospheric  nuclear  tests 

qp 

of  300  Mt  by  assuming  an  injection  of  1 x 10  H^O  molecules/Mt , 
which  is  approximately  1/2000  of  the  ambient  stratospheric  water 
loading.  For  time-dependent  estimates,  turn  to  Johnston  (1977), 
who  lists  cloud  volumes  as  a function  of  time  after  a 2-Mt  test 
(from  Christie,  1976).  These  volumes  are  shown  in  Table  B-4 
which  also  contains  estimates  of  the  fractional  augmentation  of 
local  water  content,  assuming  that  the  bomb-entrained  water  vapor 
is  uniformly  mixed  in  the  stabilized  (but  growing)  bomb  cloud; 
thus  eight  days  after  the  test  there  is  only  a 0.7  percent  en- 
hancement in  stratospheric  water  vapor. 

From  these  estimates,  it  can  be  concluded  that  detailed 

accounting  of  bomb-injected  water  vapor  is  not  warranted;  but  if 

such  calculations  are  considered  desirable,  the  molecular  ratios 

shown  in  Table  B-5  can  be  used.  These  estimates  are  based  on 

the  water  injection  of  Table  B-4  and  on  an  injection  of  1 x 10 

NO  -molecules/Mt. 
x 
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TABLE  B-3 . POTENTIAL  WATER  TRANSPORT  FROM  THE  TROPOSPHERE 

(JI/S  refers  to  Johnston  I si  and --Summer  (15°  N.  annual  mean) 
[NZ/F  refers  to  Novaya  Zemlya--Fall  (75°  N,  H (Jan  + July)] 


Place/ 

Time 

Height 

T,  °K 

P , mbar 

e mbar* 
w’ 

c ★★ 

w 

H2O  mols/Mt  at 

50%  Saturation 

JI/S 

1 km 

296 

904.2 

28.09 

0.0321 

32 

3.33  x 10 

NZ/F 

1 km 

264 

889.5 

3.097 

0.00350 

0.364  x 1032+ 

JI/S 

4 km 

277 

633 

8.13 

0.0130 

1.35  x 1032 

NZ/F 

4 km 

251 

598 

0.850 

0.0014 

0.15  x 1032+ 

Data  from  U.S.  Standard  Atmosphere,  1966  Supplements,  Table  5.1 

Smithsonian  Meteorological  Tables  (6th  Ed.,  1958)  p.  347f 

★ 

ew  = saturation  vapor  pressure  over  water 

irk 

s = water  molar  (or  volume)  mixing  ratio  in  air 
w 

tThis  is  the  saturation  vapor  pressure  over  ice  rather  than  over  water. 

NOTE:  If  a 1-Mt  air  burst,  detonated  at  1 -km  altitude,  picks  up  all  its 
water  at  the  height  of  detonation,  then  in  the  tropics  it  can  carry 
(at  50  percent  relative  humidity)  3,33  x 103^  H^O  mols/Mt  upward, 
say  into  the  stratosphere;  in  the  polar  region  it  can  carry  0.36  x 
103^  H2O  mols/Mt  upward. 

To  avoid  interaction  of  the  fireball  with  the  surface,  i.e.  to 
have  an  air  burst,  a higher  yield  burst  has  to  be  detonated  at 
higher  altitude,  e.g.,  the  58-Mt  Soviet  bomb  was  detonated  near 
4 km;  thus  it  would  bring  up,  on  the  same  scheme,  only  0.15  x 
10^2  mol s H?C/Mt.  Note,  however,  that  if  there  happened  to  be 
a rain  cloud  at  the  place  and  time  of  detonation,  so  that  the 
fireball  brings  up  the  cloud  water  as  well  as  the  atmospheric 
water  vapor,  the  total  water  transport  might  be  several  times 
higher  for,  say,  a 4-km  detonation. 


TABLE  B - 4 . LOCAL  WATER  AUGMENTATION  IN  THE 
STRATOSPHERE  FOLLOWING  A 2-Mt  TROPICAL  TEST 


Time,  days 

0 

2 

4 

8 

Volume  (102^  cm3)3 

0.58 

1 .6 

4.3 

10 

Air  mass  (103^  molecules)^ 

1 .2 

3.2 

8.6 

20 

Background  water  ( 1 0 3 3 molecules)0 

5 . 4 

14.4 

38 . 7 

90 

on  J 

Perturbation  water  (ICr  molecules) 

0.67 

0.67 

0.67 

0.67 

Percent  water  augmentation 

dJohnston  (1977),  Table  2a 
^Assuming  2 x 1013  mol/cm3  (19  km) 
cAssuming  4.5  ppmv 

1 2 

4 . 7 

1 . 7 

0.7 

^Assuming  a 2-Mt  burst  and  3.33  x 10 

33  mol/Mt 

TABLE  B - 5 . SUGGESTED  STRATOSPHERIC  H,0 
INJECTION  BY  NUCLEAR  TESTS  L 


Ratio  of  H^O 

Type  of  Test  Mo  1 ec u 1 es / NOx  Molecule 


U.S.  (1  Mt)  3 

l) . S . (large)  1 

So v i e t ( 1 M t ) 0.4 

Soviet(large)  0.2 


B . 4 . 2 Ozone  and  the  N 0 / N 0 ^ Ratio 

Johnston  (1977)  has  estimated  ozone  production  during  nu- 
clear tests  as  being  in  the  range  (-0.1  to  +3.6)  x 10^2  mole- 
cules/Mt  yield.  This  ozone  production  comes  partly  from  chemi- 
cal processes  involving  the  injected  NO^  production,  which  is 
(0.4  to  1.5)  x 10^L  molecules/Mt  yield  [see  Gilmore  (197^)]. 
Johnston  also  points  out  that  at  specified  altitudes  rather  long 
times  (from  fractions  of  a year  to  almost  two  years)  may  be  re- 
quired for  bomb-produced  MO  to  destroy  bomb-produced  ozone. 
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This  ozone  injection  is  not 


a significant  contributor 


to  the  long-term  effect  of  nuclear  tests,  since  one  NO  mole- 
cule presumably  destroys  a very  large  number  of  ozone  molecules 
(of  order  10  , depending  on  the  chemical  model  in  current  fashion) 
by  the  catalytic  cycle  of  Eqs.  (1  to  3)--see  Section  B.l.  Thus, 
the  possible  injection  of  perhaps  one  ozone  molecule  along  with 
each  NO  molecule  is  not  a significant  counter  to  process  (1  to  3). 

Further,  there  is  a question  as  to  whether  the  bomb-injected 
NO  is  injected  in  the  stratosphere  as  NO  or  as  N0„ . While  this 
is  significant  on  a short-term  basis  (cf.  Johnston,  1977),  for 
the  long-term  application  of  interest  here  it  is  not.  (H.  S. 
Johnston  and  F.  M.  Luther,  private  communications). 

B . 5 LISTING  OF  HIGH-YIELD  ATMOSPHERIC  NUCLEAR  EXPLOSIONS,  1952 

to  1976 

The  basic  data  source  is  U.S.  D0F.--see  Glasstone  (1964), 
Appendix  B,  updated  by  Cannon  (1977),  checked  against  the  Swedish 
compilation  of  Zander  and  Ar-askog  (1973  ),  a more  recent  compila- 
tion of  Carter  and  Moghissi  (1977),  and  an  informal  compilation 
prepared  by  K.  Telegadas  (NOAA-ARL). 

Here  I provide  a tentative  assignment  of  effective  yields 
for  all  listed  high-yield  atmospheric  explosions.  The  procedure 
used  to  assign  yields  is  the  following: 

1.  If  an  "official”  yield  figure  is  quoted,  that  value  is 
used . 

2.  For  the  older  (1955  to  1958)  Soviet  series,  Telegadas' 
total  estimates  divided  by  the  number  of  H-bombs  in  the 
series  are  used. 

3.  For  the  1961  to  1962  series,  the  values  of  Seitz  et  al . 
(1968)  are  used  as  base,  but  scaled  to  agree  approxi- 
mately with  total  announced  yields  (Federal  Radiation 
Council,  1963).  Thus,  for  the  1961  Soviet  series, 
where  Seitz  et  al.  list  a total  of  97  Mt , I scale  up 

24  percent  to  give  a total  yield  equal  to  the  accepted 
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figure  of  120  Mt ; for  the  1962  U.S.  Pacific  series, 
Seitz'  total  is  26  Mt , which  is  scaled  up  42  percent 
to  give  the  accepted  total  yield  of  37  Mt ; for  the  1962 
Soviet  series,  Seitz'  values  add  up  to  the  accepted 
figure  of  180  Mt . 

Table  B-6  gives  the  locations  of  the  various  test  sites  and 
Table  B-7  lists  the  yields  and  dates  of  thermonuclear  explosions 
in  the  atmosphere:  only  those  bombs  are  counted  here  that  are 
likely  to  deposit  some  NO  in  the  stratosphere,  which  means  near- 

X 

surface  devices  of  Mt  yield.  "Intermediate”  and  lower  yield 
devices  are  unlikely  to  affect  the  stratosphere  because  their 
clouds  barely  rise  above  the  tropopause,  while  high-altitude 
explosions,  which  are  detonated  at  altitudes  of  order  tens  of 
kilometers  (see,  e.g.,  Hoerlin,  1976)  produce  an  effect  on  the 
mesosphere  and  thermosphere,  but  not  on  the  stratosphere. 


TABLE  B-6.  TEST  SITES 


U .S.A. 


Eniwetok  (EN) 

11° 

N, 

162° 

Bikini  ( B K ) 

11° 

N, 

165° 

Christmas  Island  (XI) 

2° 

N, 

157° 

Johnston  Island  ( J I ) 

17° 

N, 

169° 

U.S.S.R. 

Novaya  Zemlya/Arctic  ( NZ ) 

75° 

N, 

55° 

Siberia 

52° 

N, 

O 

00 

Semi  pa  1 a t i ns  k ( SP  ) 

cn 

O 

o 

N, 

80° 

U . K. 

Chris tmas  Island  (XI) 

2° 

N, 

157° 

FRANCE 

Centre  d ' Exper i men ta t i on  du 
Pacifique  (CEP),  Tuamotu 

Archi pel  ago 

21° 

s, 

137° 

CHINA 

Lop  Nor  (LN) 

O 

o 

N, 

O 

o 

O'! 
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TABLE  B-7 . ATMOSPHERIC  THERMONUCLEAR  EXPLOSIONS 
(See  Notes  at  End  of  Table) 

Yield,  Mt 

Date  Country  Location  Name  Type  DOC Z & A C & M TE.LEGAUAS Other Suggested 


31  Oct  52 

USA 

EN : Ivy 

Mike 

Surface 

10.4 

10.4 

10.4 

10.4S 

12  Aug  53 

USSR 

? 

Atmos . 

H 

H 

H 

1 1 

28  Feb  54 

USA 

BK:Castle 

Bravo 

Surface 

15 

15 

15 

1 SS 

4 Aug  55 

USSR 

? 

Atmos. 

A 

A 

A 

0.75 

24  Sep  55 

USSR 

? 

Atmos. 

A 

A 

A 1 

1 0.75 

3 

0.75 

10  Nov  55 

USSR 

? 

Atmos. 

A 

A 

A 1 

23  Nov  55 

USSR 

? 

Atmos . 

Mt  range 

Mt  range 

>1 

0.75 

20  May  56 

USA 

BK: Redwing 

Cherokee 

Air  Drop 

Several 

Several 

j 

1 2 

27  May  56 

USA 

BK: Redwing 

Zuni 

Surface 

3.5 

3.5 

3,3 

i ’7  3.5S 

20  Jul  56 

USA 

BK: Redwing 

Tewa 

Barge 

5 

'-5 

5.01 

1 *S 

30  Aug  56 

USSR 

? 

Atmos . 

Large 

Large 

0.2-1  - 

2 Sep  56 

USSR 

? 

Atmos. 

A 

A 

A 1 

f 

1.5 

6 

1.5 

10  Sep  56 

USSR 

? 

Atmos . 

A 

A 

A \ 

1 

17  Nov  56 

USSR 

? 

Atmos. 

Large 

Large 

0.2-1  > 

' 1.5 

10  Apr  57  USSR  ? Large  Large  0.2-1  1 3 

6 

16  Apr  57  USSR  ? Large  Large  0.2-1  ] 3 

15  May  57  UK  XI:Grapple  Air  Drop  Mt  ranae  Mt  ranqe  >1  I 

31  May  57  UK  XIrGrapple  Air  Drop  Mt  range  Mt  range  >1/3  1 


19  Jun  57 

UK 

XI :Grapple 

Air  Drop 

range 

range 

] 

1 

22  Aug  57 

USSR 

Siberia 

Substantial 

Substantial 

>1 

2 

24  Sep  57 

USSR 

Arctic 

Mt  range 

Mt  range 

>1  6 

2 

6 Oct  57 

USSR 

? 

Substantial 

Substantial 

H 

2 

8 Nov  57 

UK 

XI : Grapple 

Air  Drop 

Mt  range 

Mt  range 

>1 

1 

23  Feb  58 

USSR 

Arctic 

Mt  range 

Mt  range 

>1  ) 

3 

27  Feb  58 

USSR 

Arctic 

Mt  range 

Mt  range 

» | 10 

3 

28  Apr  58 

UK 

XI tGrapple 

Air  Drop 

Mt  range 

Mt  range 

>i  i 

1 

12  May  58 

USA 

EN:Hardtack 

Koa 

Surface 

1.37 

1.37 

1.37  E 

1.4S 

28  Jun  58 

USA 

EN:Hardtack 

Oak 

Barge 

8.9 

8.9 

8.9  E 

8.9S 

2 Sep  58 

UK 

XI rGrapple 

Air  Drop 

Mt  range 

Mt  range 

.1 

0 

I 

11  Sep  58 

UK 

XI :Grapple 

Air  Drop 

Mt  range 

Mt  range 

c 

>1 

1 

10  Oct  58 

USSR 

Arctic 

Atmos . 

Rel . large 

Mt  range 

>1  \ 

1 

12  Oct  58 

USSR 

Arctic 

Atmos. 

Large 

Mt  range 

>1  1 

3 

15  Oct  58 

USSR 

Arctic 

Atmos. 

Large 

Mt  range 

>1  / 

3 

18  Oct  58 

USSR 

Arctic 

Atmos. 

Large 

Mt  range 

>]  \ 26 

3 

20  Oct  58 

USSR 

Arctic 

Atmos . 

Large 

Mt  range 

>1  ( 

3 

22  Oct  58 

USSR 

Arctic 

Atmos. 

Large 

Mt  range 

>1  \ 

3 

24  Oct  58 

USSR 

Arctic 

Atmos . 

Large 

Mt  range 

> 1 

3 

25  Oct  58 

USSR 

Arctic 

Atmos . 

Rel.  Large 

Rel.  Large 

0.2-1  / 

1 
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TABLE 

B-7 . 

ATMOSPHERIC  THERMONUCLEAR  EXPLOSIONS 

(Cont. 

P.  2) 

Yield, 

Mt 

Udte 

Country 

Location 

Name  Type 

DOf 

Z & A 

f.  X M 

TLL16ADAS 

Other 

• • : 

Order  of 

10  Sep  61 

USSR 

NZ 

Atmos. 

Several  Mt 

Several  Mt 

2. 4 (Seitz) 

3 

12  Sep  61 

USSR 

NZ 

Atmos. 

Several  Mt 

Severa I Mt 

' 

?.4(Se’ tz  1 

3 

14  Sep  61 

USSR 

NZ 

Atmos . 

Several  Mt 

Several  Mt 

2. 4 (Seitz 

3 

16  Sep  61 

USSR 

NZ 

Atmos. 

vi 

Mt  Range 

C.v(Seitz) 

10  Sep  61 

USSR 

NZ 

Atmos. 

%1 

Mt  Range 

0.5(Seit 2 

V 

20  Sep  61 

USSR 

NZ 

Atmos. 

•vl 

Mt  Range 

J.6(  Seitz 

1 

22  Sep  61 

USSR 

NZ 

Atmos. 

M 

Mt  Range 

1 

. 5 ( Se  1 1 z l 

2 Oct  o! 

USSR 

NZ 

Atmos. 

1 

Order  of 

Mt  Ranee 

\ 120 

0.5( Seitz  1 

4 Oct  61 

USSR 

NZ 

Atmos. 

Several  Mt 

Several  Mt 

. ] 

2.4(Seitz . 

3 

6 Oct  61 

USSR 

N/ 

Atmos. 

Several 

Several 

?.4lSeitz 

20  Oct  61 

USSR 

NZ 

Atmos. 

Severa 1 

Severa 1 

2. 4 i Seitz ; 

23  Oct  61 

USSR 

NZ 

Atmos . 

'.25 

25 

25 

2 KSei  tz ) 

25 

30  Oct  61 

USSR 

NZ 

Air  Durst 

58 

58 

58 

55(Seitz 

58 

31  Oct  bl 

USSR 

NZ 

Atmos . 

Sever-'.  1 

Several 

1 

2 . 4 ( Se  i t z ) 

3 

m Nov  61 

USSR 

NZ 

Atmos. 

Several 

Several 

1 

2. 4 (Seitz) 

2 May  62 

USA 

XI 

Arkansas  Air  burst 

Low  Mt 

1-10 

>1 

2 (Seitz; 

10  Jun  62 

USA 

XI 

Yeso  Air  burst 

Low  Mt 

1-10 

2 (Seitz) 

27  Jun  62 

USA 

XI 

Bighorn  Air  burst 

Mt  Range 

1 

1 

32 

6 (Seitz 

U 

30  June  62 

USA 

XI 

Bluestone  Air  burst 

Low  Mt 

1-10 

2 'Seitz) 

3 

11  Jul  62 

USA 

XI 

Pamlico  Air  burst 

Low  Mt 

1-10 

1 

2 (Seitz) 

3 

5 Aug  62 

USSR 

NZ 

30 

•30 

30 

0 

30  Seitz 

30 

Order  of 

20  Aug  62 

USSR 

NZ 

Several  Mt 

Several 

D 

7 (Seitz.! 

7 

22  Aug  62 

USSR 

NZ 

Low  Mt 
Order  of 

wt  Range 

M 

0 

3 (Seitz) 

3 

25  Aug  62 

USSR 

NZ 

Several  Mt 

Several 

1 

D 

7 (Seitz) 

7 

2/  Aug  62 

USSR 

NZ 

Severa 1 

Several 

•1 

D 

7 (Seitz, 

7 

8 Sep  62 

USSR 

NZ 

Mt  Range 

Mt  Range 

•1 

D 

1 (Seitz) 

1 

15  Sep  62 

USSR 

NZ 

Several  Mt 

Several  Mt 

>1 

0 

• ' 2 

7 

16  Sep  62 

USSR 

NZ 

Several  Mt 

Several  Mt 

•1 

0 

7 (Seitz) 

7 

Id  Sep  62 

USSR 

NZ 

A few  Mt 

A few  Mt 

1 

D 

3 (Seitz) 

3 

19  Sep  62 

USSR 

NZ 

Mult i -Mt 

Multi-Mt 

1 

D 

24  (Seitz) 

24 

21  Sep  62 

USSR 

NZ 

A few  Mt 

A few  vt 

•1 

D 

3 

25  Sep  62 

USSR 

NZ 

Atnos. 

Multi-Mt 

Multi-Mt 

•1 

D 

27  (Seitz) 

27 

27  Sep  62 

USSR 

NZ 

Atmos. 

<30 

30 

< 3U 

24  (Seitz) 

24 

18  Oct  02 

USA 

J.I. 

Chama  Air  burst 

Lot  Mt 

1-10 

>1 

C 

* (Seitz) 

3 

22  Oct  62 

USSR 

NZ 

Atmos. 

Several  Mt 

Several  Mt 

•1 

(Seitz) 

30  Oct  62 

USA 

J.I. 

Housatonic  Air 

Mt-range 

-I 

>1 

C 

8 (Seitz) 

lo 

25  Dec  62 

USSR 

NZ 

Atmos. 

B 

B' 

B' 

3 (Seitz) 

:c,  3 

17  Jun  67 

China 

Lop  Nor 

Air  drop 

•'-3(H) 

3(H) 

v3(H) 

3 

3 
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TABLE  B-7.  ATMOSPHERIC  THERMONUCLEAR  EXPLOSIONS  (Cont.  p.  3) 


Date 

Country 

Location 

Name 

Type 

DOE 

7 A A 

field. 

C A M T l : E .ADA' 

Other 

ugqestf  I 

_ 4 Aui;  ’ •* 

France 

CEP 

Canopus 

Bal loon 

Low  Mt 

2.5 

I(H  ( 

2.5 

: ep  bb 

* ranee 

CEP 

Procyuf- 

Bal loon 

1 >W  Mt 

1 

( 

1 

i.7  Dec  bo 

Ch  i na 

LCq  Nor 

Air  crop 

-3 

•■3(h)  3 

3 

i-i  Sep  69 

China 

Lop  'ior 

Air  drop 

3 

3 3 

3 

JO  Mdy  7 , 

F rar.  « 

Ora  non 

Ea  l loon 

I m t * r <* 

0.1-1 

1 

1 

3 Jut  70 

r ranee 

CEP 

Licorne 

Bal  loon 

If  termed. 

1 

i * 

1 

1-  7d 

">« 

, '• 

Air  drop 

%3  3 

3 

u ’I 

frar.  e 

Rhea 

High 

1 

0.2-1 

i 

» •>*■"  ’> 

China 

Atmos . 

2-3 

3 3 

2.5 

*.OV  7t 

China 

• 

At-  os. 

4 

- 

• 

- ^ ' • 1 j •!  - 1 1 t i ;ufie  ' *•  s are  exc I uded . 

*tes  .iff  sometimes  date  of  AEC  announcement  rather 

t *i  «»  ,'iut  late. 

» often  not  specified  m the  earlier  events. 

is  generally  not  specified  'i.e.,  air  drot,,  balloon, 
- r f a : r‘ , e t : . T h j s one  ay  • s s . ^ that  the  effective 
• t r Or responds  to  the  full  yield  rather 
tna»*  only  50'  as  tor  a Surface-  nteractive  bur  t 
-u fc-f  I- -.a lions  ’isted  on  Table  B-4. 

Meld 

C ■ m ^ a ' L»lo?  stone  '1964),  s u p p ! ero  ented  w i 1 1 . 1 a n no  n 1977  , 
l i A •’ears  Zander  A ra  sleep  { T 9 7 ^ ' 

. . A M 'I’eans  Carter  *.  Mo<;hissi  (197F', 

Teiegadas  eans  prvate  communication,  April  1977 
; e i 1 4 mears  Setiz  et  al.  ( 1 9 1 a > 

■ dash.  ( - 1 means  yield  not  specified 

A blink  ;pa  means  detonation  not  listed  -typically,  event 
after  publication  of  reference 

ear?  y Id  :•  o t specified:  "ay  not  : p thernionucl ea r “ 

H means  ' tne**monucl  ea  r ; yield  rot  specified 

B neans  AE  announcement:  a number  of  a tmc  . l h e r i c tests  held 
23-25  Ore.  1 962  , largest  about  20  Mt,  others  low  to  a few  Mt" 

B irei-s  .0  Mr;  others  listed  without  yield 
l m-.ans  Teiegadas  quotes  4 Mt  for  Chama  ♦ Housatonic 

means  " surf ace- i nteract i ve  shot:  rise  height  correspon  ‘ o 
yield  quoted,  *«P  injection  should  correspond  to  half  the 
effective  yield1 

D -eans  Teiegadas  quotes  total  yifld  of  180  Mt  for  Soviet 
1962  Series,  in  agreement  with  Federal  Radiation  Council'' 

E means  Teieyddas  quotes  a total  yield  of  Mt  for  US  tests 
from  ZB  Aoril  to  26  July  1958" 

r:*-  t1-  t '..lo,-,!1  Pesearch  Council  M 96  3:  li  ts  total  yields  of  1961 
tests  a.  3 7 V t ; USSR,  180  Mt 


i - ■ 


For  general  orientation,  Table  B-8  provides  the  total  mega- 


tonnage of  the  various  atmospheric  test  series,  Table  B-9  gives 

# 

the  groupings  used  by  Seitz  et  al.  (1968)  of  the  high-yield 
stratospheric  injections  for  the  1961  to  1962  series,  and  Table 
B-10  gives  an  aggregated  list  of  ''intermediate-range"  atmos- 
pheric nuclear  explosions. 

The  following  comments  should  be  made: 

a.  These  data  are  all  taker,  from  published  information  in 
the  open  literature. 

b.  The  data  are  clearly  incomplete  and  inconsistent — see 
the  procedure  of  items  1 to  3 above  for  "assigning" 
yields  . 

c.  There  have  been  a large  number  of  atmospheric  nuclear 
explosions  with  energy  significantly  below  1 Mt , which 
tend  to  stabilize  near  the  local  tropopause  according 
to  current  cloud  rise  models.  They  will  not  reduce  the 
ozone  level  very  much;  in  fact,  as  a result  of  the 
"smog"  reactions  they  may  even  enhance  the  ozone  level. 
Table  B-ll  is  an  attempt:  to  put  these  various  factors  in 
context  by  correlating  the  mean  size  of  explosions  with 
their  cloud  altitude  and  thus  with  their  effect  on  ozone; 
it  is  clear  that  while  there  have  been  very  many  "small" 
bombs,  the  overwhelming  megatonnage  and  thus  NC>x  was 
carried  by  the  "relatively  few"  H-bombs,  and  therefore  it 
is  not  appropriate  to  concern  oneself  unduly  with  the 
production  of  ozone  from  the  "small"  bombs  through  the 
tropospheric  smog  reactions. 

d.  The  overall  division  into  "small,"  "intermediate,"  and 
"large"  bombs  is  somewhat  ambiguous.  I have  tried  to 
count  every  announced  event  in  one  and  only  one  cate- 
gory, except  for  sub-kt  events,  which  are  ignored. 

_ 

Seitz  et  al.  (1968)  studied  the  global  transport  of  radioactive 
debris  (Sr-90)  from  the  various  high-yield  tests  in  the  1961  to 
1962  series.  Thus,  their  assignment  of  injection  heights  and 
yields  is  of  some  relevance  for  the  present  injection  of  oxides 
of  nitrogen. 
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GENERAL  ORIENTATION:  ESTIMATED  TOTAL  YIELD 


1976  China  (Mid- Lati tude) 


TABLE  B - 9 . INJECTION  GROUPINGS  OF  1 96  1 -1  962 
TEST  SERIES  USED  BY  SEITZ  et  al  . ( 1 968) 


Injection 

Dates 

Location 

1 

10,  12,  14,  16,  18,  20, 

22  September  1961 

U.S.S.R. 

(Novaya  Zemlya) 

2 

2 , 4 , 6 October  1 961 

U.S.S.R. 

(Novaya  Zemlya) 

3 

20  , 23  , 30  , 31  October  1961; 

4 November  1961 

U.S.S.R. 

(Novaya  Zemlya) 

4 

2 May  1962 

U .S .A. 

(Chris  tma  s Island) 

5 

10  June  1962 

U.S.A. 

(Christmas  Island) 

6 

27,  30  June  1962 

U.S.A. 

(Chris  tma  s Island) 

7 

11  July  1962 

U.S.A. 

(Christmas  Island) 

8 

5 August  1962 

U.S.S.R. 

(Novaya  Zemlya) 

9 

20,  22,  25,  27  August  1962 

U .S.S.R. 

(Novaya  Zemlya) 

10 

8 September  1962 

U.S.S.R. 

(Novaya  Zemlya) 

1 1 

15,  16,  18,  19,  21 , 25, 

27  September  1962 

U.S.S.R. 

(Novaya  Zemlya) 

12 

18,  30  October  1962 

U.S.A. 

(Johnston  Island; 

22  October  1962 

U.S.S.R. 

(Novaya  Zemlya) 

1 3 

24  December  1962  (two  tests) 

U.S.S.R. 

(Novaya  Zemlya) 
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Location 

— 

Type 

Number/Yield  in  Lt 

IRDA 

Z & A 

C4M 

Nov  52 

USA 

EN-.lvy 

Air  Drop 

!/(S00) 

l/SOO 

l/(  -200) 

Jun  54 

USA 

BK:Castle 

Surface 

1/100 

1/100 

1/110 

Aug-Nov  56 

USSR 

? 

- 

- 

3 

Apr 

USSR 

Various 

- 

- 

- 

3 

Mar  58 

USSR 

NZ 

- 

3/ ( < 1 000 ) 

o 

o 

o 

ro 

10 

Sep-Oct  61 

USSR 

SP 

- 

6 

6 

6 

Sep-Nov  61 

USSR 

NZ 

- 

7 

7 

7 

Apr-Jul  62 

USA 

XI : Dominic 

Air  Drop 

1 7 

1 7 

1 7 

Aug-Sep  62 

USSR 

NZ 

- 

2 

2 

2 

Oct  62 

USA 

JI  :Dominic 

Air  Drop 

2 

2 

2 

Oct -Dec  62 

USSR 

NZ 

- 

10 

10 

10 

May  66 

China 

Lop  Nor 

Air  Drop 

l/( lower  end  of  I 
Range) 

1/(200-500) 

1/(20-60) 

Sep  66 

France 

CEP 

Bal loon/Barge 

2(smal 1 ) 

2/(120.150) 

low  yield  only 

Oct  66 

France 

CEP 

Barge 

1(200-300) 

1/(200-300) 

1/(200-300) 

Dec  66 

China 

Lop  Nor 

Tower 

1 / ( few  100  kt) 

1/300 

l/(>200) 

Jul-Aug  68 

France 

CEP 

Balloon 

2/(500, <1000) 

2/(500, <1000) 

2 

May-Aug  70 

France 

CEP 

Balloon 

3 

3 ( <20 ) 

5 

Jul  71 

France 

CEP 

Balloon 

1 

1/(400-500) 

1/(20-200) 

Mar  72 

China 

Lop  Nor 

1/(20-200) 

1/(20-200) 

1/(20-200) 

Jul  74 

China 

Lop  Nor 

1 

1/(200-1000) 

Sep  76 

China 

Lop  Nor 

- 

1/(20-200) 

- 

NOTES:  Notation  is  same  as  for  Table  B-5 
Intermediate  range  is  20  to  1000  kt 

This  is  a greater  level  of  aggregation  than  in  Table  B-5 
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TABLE  B - 1 1 . AGGREGATED  EFFECTS  OF  BOMB - PRODU C E D 
NOx  ON  STRATOSPHERIC  OZONE 


Yield,  Rang** 

Total 

wona tonnage 
(1 050-197?) 

Total  fio,  of 
Bombs* 
(1950-1976) 

A1 1 i tude  of  In  ’ect  ion 

of  c tratcspher ic  r?  nr* 

Low:  1 -20  kt 

-1C 

1°0 

tropo pause 

Increases  ozr.re  by  smog  reactions 
but  is  lost  rapidly  by  ra incut. 

I ntermeri iate : 
20-1000  h t 

-10 

f 4 

iust  above  tropopause 

Rather  little. 

High:  ^.1 

-600 

87 

in  stratosphere 

Reduces  ozone  tv  rata lytic  cycle- 
lono  1 i fe  t ire. 

Note: 

The  breakdown  by 
abovr-  10  km)  are 

countries 

omitted. 

is  shown 

below.  High-altitude 

explosions  (henht  of  burst 

Country 

A-Hombs  (1-. 

IP  kt) 

(20-l°'J*  h-Bon»bs  ! 1 ,f\  > 

124 

21  14 

1-146  Mt) 

U.S.S.R. 

36 

28  66 

'-If 0 Mt) 

U.K. 

11 

7 

' -7  ••t) 

Trance 

16 

9 c 

' -7  "t) 

China 

6 

6 ( 

: -if  i't  i 

Note  that  high-altitude  explosions  are  omitted. 
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B . 6 SUGGESTED  STRATOSPHERIC  NOx  INJECTIONS  DUE  TO  PAST 
ATMOSPHERIC  NUCLEAR  EXPLOSIONS 

Table  B-12  lists  my  best  estimates  of  stratospheric  NO 

x 

injections  due  to  atmospheric  nuclear  explosions.  The  estimates 
were  derived  by  combining  the  various  input  data  and  models  dis- 
cussed previously.  For  orientation,  Fig.  B-5  shows  the  NO^ 
injection  into  the  stratosphere  by  months  from  1952  through  1976 

in  terms  of  megatons  of  effective  yield  using  the  conversion 

B2 

factor  (0.5  to  1)  x 10  NO  molecules/Mt . For  the  ambient 

X g 

stratospheric  NO^  loading  we  use  the  value  of  5 x 10  tons 
(as  HNCg)  from  A.  J.  Grobecker  et  al . (1974),  p.  10.  By  far 
the  most  significant  NOx  injections  were  the  Arctic  injections 
of  1961  to  1962,  which  exceed  10  percent  of  the  global  strato- 
spheric NOx  loading  for  several  months  and  thus  might  be  expected 
to  show  observable  effects  on  total  ozone,  at  least  at  high 
latitudes . 
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TABLE  B-12.  STRATOSPHERIC  INJECTION  OF  NO, 


BY  NUCLEAR  EXPLOSIONS 


Country  Location 


Effective 
Yield,  Mt 


NO*  Injection 

1032  mo  1 s Altitude  karuj*-, 


TABLE  B-12.  STRATOSPHERIC  INJECTION  OF  NO 
BY  NUCLEAR  EXPLOSIONS  (Cont.  p.  2) 


Date 

Country 

Locat ion 

Effective 
Yield,  Mt 

N0X  Injection 

1 0^  mo  1 b 

Altitude  Range,  km 

4 Nov  61 

USSR 

N7 

3 

3 

10  - 20.7 

2 Way  62 

USA 

XI 

3 

3 

14.7  - 23.5 

10  Jun  62 

USA 

XI 

3 

3 

14.7  - 23.5 

27  Jun  62 

USA 

XI 

11 

11 

18.1  - 30.9 

30  Jun  62 

USA 

XI 

3 

3 

14.7  - 23.5 

11  Jul  62 

USA 

XI 

3 

3 

14.7  - 23.5 

5 Aug  62 

USSR 

NZ 

30 

30 

14.4  - 29.6 

20  Aug  62 

USSR 

NZ 

7 

7 

11.3  - 23.7 

22  Aug  62 

USSR 

NZ 

3 

3 

10  - 20.7 

25  Aug  62 

USSR 

NZ 

7 

7 

11.3  - 23.7 

27  Auq  62 

USSR 

N7 

7 

7 

11.3  - 23.7 

8 Sep  62 

USSR 

NZ 

1 

1 

7.9  - 16.0 

15  Sep  62 

USSR 

NZ 

7 

7 

11.3  - 23.7 

16  Sep  62 

USSR 

NZ 

7 

7 

11.3  - 23.7 

18  Sep  62 

USSR 

NZ 

3 

3 

10  - 20.7 

19  Sep  62 

USSR 

NZ 

24 

24 

13.9  - 29.1 

21  Sep  62 

USSR 

NZ 

3 

3 

10  - 20.7 

25  Sep  62 

USSR 

NZ 

27 

27 

14.1  - 29.2 

27  Sep  62 

USSR 

N7 

24 

24 

13.9  - 29.1 

18  Oct  62 

USA 

JI 

3 

3 

14.7  - 23.5 

22  Oct  62 

USSR 

NZ 

7 

7 

11.3  - 23.7 

30  Oct  62 

USA 

Jl 

10 

10 

17.9  - 30.0 

23  Dec  62 

USSR 

NZ 

20 

20 

13.6  - 28.0 

25  Dec  62 

USSR 

NZ 

3 

3 

10  - 20.7 

17  Jun  67 

CHINA 

U 

3 

3 

13.1  - 22.7 

24  Auq  68 

FRANCE 

CEP 

2.5 

2.5 

14.4  - 22.9 

8 Sep  68 

FRANCE 

CFP 

1 

1 

12.5  - 19.5 

27  Dec  68 

CHINA 

LN 

3 

3 

13.1  - 22.7 

29  Sep  69 

CHINA 

LN 

3 

3 

13.1  - 22. 7 

30  Way  70 

FRANCE 

CEP 

1 

1 

12.5  - 19.5 

30  Jun  70 

FRANCE 

CEP 

1 

1 

12.5  - 19.5 

14  Oct  70 

CHINA 

LN 

3 

3 

13.1  - 22.7 

14  Aug  71 

FRANCE 

CEP 

1 

1 

12.5  - 19. j 

26  Jul  73 

CHINA 

LN 

2.5 

2.5 

12.5  - 21.6 

17  Nov  76 

CHINA 

LN 

4 

4 

13A  - 24.7 

NOTES:  For  locations  see  Table  2.2  where  identified,  otherwise: 

?M  (USSR,  1953-1957)  take  50-5?°N,  78-8fl°E 
Arc  (USSR,  1958)  take  HZ:  ?S°N,  55°E 

For  effective  yields,  S denotes  surface  burst:  heioht  corresponds  to 
actual  yield,  NO^  injection  is  half  that  of  equivalent  air  burst. 

1? 

For  NO-  injection,  see  Table  B-l.  I suggest  1x10  NO  mols/Mt  as  a 
plausible  choice.  x 

For  altitude  range,  see  Fig.  B*5.  For  equatorial  bursts  I suggest 
Peterson's  (1970)  curves,  for  polar  bursts,  Seitz  (1968)  curves  as 
shown  there.  Table  B-2  from  Peterson  (1970)'  is  the  recommended 
NO  injection  within  the  altitude  range  quoted  here. 
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ANTHROPOGENIC  RELEASE  OF  HALOCARBONS 
IN  THE  ATMOSPHERE,  1955  TO  1975 

C.l  SPECIES  CONSIDERED 

Table  C-l  lists  a number  of  chlorof luoro-  and  chloro- 
carbons  which  are  of  possible  significance  as  a source  of  free 
chlorine  atoms,  or,  more  generally,  of  C1X  (=C1  + CIO  + HC1)  in 
the  stratosphere,  and  which  may  thus  serve  to  reduce  the  concen- 
tration of  stratospheric  ozone.  Bromine-containing  molecules 
are  discussed  in  Section  C.4  below.  The  standard  notation  is 
used  for  chlorof luorocarbons : thus,  in  molecule  FC-xyz,  the 
hundreds  digit  "x"  is  the  number  of  carbon  atoms  in  the  molecule 
minus  one,  the  tens  digit  "y"  is  the  number  of  hydrogen  atoms  in 
the  molecule  plus  one,  and  the  units  digit  "z"  is  the  number  of 
fluorine  atoms  in  the  molecule.  The  undesignated  atoms  are 
chlorine,  and  if  x=0,  the  first  digit  is  dropped.  Thus,  FC-11 
is  C F Clg , etc  . 

In  Table  C-l,  I quote  the  mass  of  material  that  contains 

34 

10  chlorine  atoms,  since,  presumably,  it  is  these  atoms  that 
deplete  the  ozone.  For  orientation,  if  the  mixing  ratio  of  the 
C1X  species  is  q ppbv  (=  q x 10~"  molecules  per  "air  molecule"). 

this  corresponds  to  a total  stratospheric  loading  of  q x 2.47  x 

34  * 

10  molecules. 

Table  C-l  also  lists  estimated  tropospheric  injection 
rates  and  mean  tropospheric  lifetimes  for  the  species  quoted  to 

* 

This  assumes  a global  11-km  tropopause  and  the  number  densities 
<f  J.3.  Standard  Atmosphere,  1976. 
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provide  a crude  estimate  of  their  relative  impact  on  the  strato- 
sphere. Evidently,  a species  with  a snail  atmospheric  injection 
rate  and  a relatively  short  lifetime  in  the  troposphere  is  less 
likely  to  contribute  to  stratospheric  chlorine  than  do  FC-11  and 
FC-12 , with  relatively  large  injection  rates  and  no  established 
photochemical  sinks  in  the  troposphere. 

The  specific  tropospheric  lifetimes  of  various  materials 
are  due  to  reaction  with  OH  radicals.  There  is  a current  con- 
troversy over  whether  the  tropospheric  concentration  of  these 

6 *2 

radicals  is  (1  - 3)  x 10  /cm  , which  is  the  range  of  measured 

5 ' 3 

values,  typically  at  midday,  as  against  *4  x 10  /cm  on  a di- 
urnal average,  which  gives  the  range  of  5 to  10  in  effective 
tropospheric  lifetimes  shown  in  Table  C-l  (see,  e.g.,  Singh, 
1977). 

C . 2 ANTHROPOGENIC  EMISSION  RATES  OF  CHL0R0FLU0R0-  AND  CHL0R0- 

CARBONS,  1955-1975 

The  production  rates  of  halocarbons  derived  from  different 
sources  show  some  variation,  and  thus  rather  detailed  explana- 
tions are  presented  of  how  the  data  given  here  are  derived  from 
earlier  sources. 

The  emission  rates  present  an  even  more  complicated  prob- 
lem, which  has  been  discussed  by  McCarthy  et  al.  (1977)  for  the 
fluorocarbons.  Again,  detailed  descriptions  are  included  of  how 
the  present  estimates  were  derived  from  "primary"  sources. 

For  fluorocarbons  FC-11,  FC-12,  and  FC-22,  the  latest  esti- 
mates for  the  global  production  and  release  rates  are  given  in 
Table  C-2.  Explicitly,  I have  used  the  primary  data  as  com- 
piled by  A.  Grant  (1977),  who  used  as  a data  base  the  records 
of  some  20  companies  (plus  subsidiaries)  that  reported  produc- 
tion and/or  release  data  [see  Table  1 of  McCarthy  et  al . (1977)] 
The  data  of  Grant  (1977)  for  FC-11  and  FC-12  have  been  analyzed 
by  E.  I.  duPont  de  Nemours  (1977)  to  provide  worldwide  produc- 
tion and  release  data. 


TABLE  C-2.  GLOBAL  PRODUCTION  AND  RELEASE  RATES  FOR 


FC-11  , -12, 

FC-11 

and  -22 

THROUGH 

FC- 

1976  (in 

12 

1 0 6 k g ) 

FC 

-22 

Year 

Production 

Rel ease 

Production 

Rel ease 

Production 

Release 

1931 

0.5 

1932 

0.1 

1933 

0.3 

0.1 

1934 

0.7 

0.2 

1935 

1 .0 

0.2 

1936 

0.1 

1.7 

0.4 

1937 

0.1 

3.1 

0.6 

1 938 

0.  1 

2.8 

0.9 

1939 

0.1 

0.1 

3.9 

1.8 

1 940 

0.2 

0.1 

4.5 

1 . 7 

1941 

0.3 

0.  1 

6.3 

2.3 

1942 

0.3 

0.1 

5.9 

2.9 

1 943 

0.  4 

0.2 

8.2 

3.5 

1944 

0.4 

0.2 

16.  7 

4.7 

1945 

0.4 

0.3 

20.  1 

6.1 

0.1 

1946 

0.  7 

0.6 

16.6 

11.9 

1 947 

1 . 3 

1.2 

20.1 

19.0 

0.1 

1948 

3.0 

2.3 

24.8 

22.2 

0.2 

1 949 

4.5 

3.7 

26.1 

24.2 

0.3 

0.1 

1950 

6.6 

5.4 

34.9 

27.4 

0.8 

0.1 

1951 

9.  1 

7.5 

37.0 

30.8 

1 . 0 

0.2 

1952 

13.7 

10.9 

38.5 

32.6 

1.6 

0.  3 

19  3 3 

17.6 

14.9 

48.2 

37.0 

2.2 

0.4 

1954 

21.4 

18.7 

51  . 3 

42.2 

2.9 

0.6 

1955 

27.0 

23.2 

60.3 

47.5 

3.5 

0.8 

1956 

33.5 

29.  1 

71.9 

55.4 

6.1 

1.2 

1957 

35.3 

32.8 

77.9 

63.1 

6.4 

1 . 1 

1 958 

31.3 

31.2 

77.8 

66.4 

7.7 

1.4 

1959 

37.9 

32.3 

92.5 

73.9 

11.4 

2.7 

1 960 

52.  5 

42.2 

105.0 

88.1 

12.3 

3.3 

1 961 

63.9 

54.1 

114.7 

98.6 

12.3 

3.6 

1962 

82.1 

67.6 

135.0 

113.1 

15.6 

4.7 

1 963 

98.0 

82.6 

154.0 

132.5 

17.8 

5.7 

1964 

116.6 

97.9 

178.4 

153.8 

22.6 

6.9 

1 965 

129.1 

111.7 

199.2 

173.6 

25.4 

8.2 

1 966 

148.1 

125.1 

226.1 

192.9 

31  .9 

9.9 

1967 

167.8 

141.9 

253.5 

217.6 

37.8 

12.5 

1968 

192.1 

161.6 

279.3 

248.9 

46.1 

16.5 

1969 

227.3 

187.0 

310.0 

271  . 1 

56.6 

20.6 

1970 

249.2 

212.3 

334.8 

296.1 

59.3 

23.5 

1971 

275.5 

233.2 

356.3 

317.7 

65.2 

26.5 

1 972 

315.8 

260.3 

391  .1 

343.5 

70.8 

31  . 3 

1973 

359.5 

295.6 

434.7 

378.0 

77.4 

36.1 

1974 

379.7 

324.1 

455.0 

408.4 

88.7 

42.3 

1975 

326  .5 

314.9 

396.1 

396.3 

74.4 

47.6 

1 976 

355.2 

308.2 

427.9 

380.1 

94.1 

60.2 

Note:  For  FC-11  and  -12  the  data  from  Table  5 of  DuPont  (1977).  For  FC-22,  the  production 
figures  are  from  Grant  (1977)  corrected  for  global  production  as  indicated  in  the 
text  (Section  2).  The  FC-22  release  figures  are  obtained  from  the  production  figures 
by  multiplying  by  the  release/production  ratio  from  Table  7 of  McCarthy  et  al . (1977) 
on  a year-by-year  basis;  but  for  1976  the  1975  ratio,  47.1/73.6  * 0.640,  was  used. 
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For  FC-22 , the  production  data  are  from  Grant  (1977)  for 
the  20  primary  companies,  corrected  by  a factor 

, , 11.5  + 9.9 

1858.9  + 11.5  + 9.9 

on  an  annual  basis.  The  terms  here  are  the  following: 

1858.9  = cumulative  production  of  FC-22  by  the  20  re- 
porting companies  from  the  start  of  produc- 
tion through  1976 

11.5  = cumulative  production  of  F>C-22  estimated 
for  Eastern  Bloc  countries 

9.9  = cumulative  production  of  FC-22  in  India  and 
Argentina 

[These  numbers  come  from  Grant  (1977),  Schedules  4, 

11 , and  12 . ] 

The  release  data  for  FC-22  are  obtained  from  the  production 
data  obtained  in  this  way,  by  multiplying  them  by  the  ratio  of 
release/production  on  an  annual  basis  from  McCarthy  et  al . 
(1977),  Table  7;  except  that  for  1976  the  correction  factor 
used  is  from  McCarthy  et  al.  (1977)  for  1975,  which  is  47.1/73*6 
0.640. 

The  overall  release  data  for  FC-11,  FC-12 , and  FC-22  are 
shown  in  Fig.  C-l  as  functions  of  time. 

For  the  Soviet  Bloc,  the  estimates  of  Schedule  12  of  Grant 
(1977)  have  been  used  here.  Presumably,  these  estimates  should 
be  updated  by  the  new  data  of  Borisenkov  and  Kazakov  (1977) 
given  in  Table  C-3,  which  indicate  that  U.S.S.R.  production  rep- 
resents 2.5  percent  of  the  world  total  in  aggregate,  as  compared 
with  4.5  percent  for  the  whole  Eastern  Bloc  (including  China), 
as  the  earlier  estimate  of  Grant  (1977)  indicated.  However, 
these  quantities  lie  within  the  error  bounds  of  the  current 
estimates . 

For  purposes  of  comparison,  the  estimate  of  Schedule  12  of 
Grant  (1977)  for  total  Soviet  Bloc  production  from  1950  to  1976 
is  FC-11,  172  x 106  kg;  FC-12,  226  x 10£  kg.  Note  that  the 
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FIGURE  C-l.  Global  Emissions  of  Fluorocarbons  -11 


listed  production  figures  for  the  U.S.S.R.  are  lower  than  the 

# 

estimates  for  the  Eastern  Bloc,  especially  for  FC-11. 

TABLE  C-3 . FLUOROCARBON  PRODUCTION  IN  THE  U.S.S.R. 

1968-75  (BORISENKOV  AND  KAZAKOV,  1977)  in  106  kg/yr 


Compound 

1968 

1969 

1970 

1971 

1972 

1973 

1974  1975 

FC-11 

1 .4 

2.2 

2.5 

2.9 

3.6 

4.0 

6.1  7.4 

FC-12 

9.6 

12.0 

13.3 

15.9 

18.2 

19.8 

25.2  31.2 

Total 

11.0 

14.2 

15.8 

18.8 

21  .8 

23.8 

31.3  38.6 

Total  for 

1968-1975 

: FC-11  - 

30.1 

FC-12  - 

145.1 

Grand 

total 

- 175.2 

U.S.S.R.  share  in  world  production  for  1968-75  is  1.3%  for  FC-11, 
4.9%  for  FC-12 . 

In  1974-75  about  45%  of  the  total  U.S.S.R.  production  of  FC-11  and 
FC-12  was  consumed  in  the  manufacture  of  consumer  chemicals,  as 
against  50%  in  the  U.S.A.  and  80%  in  Great  Britain. 


For  carbon  tetrachloride,  different  estimates  for  produc- 
tion and  release  rates  are  shown  in  Table  C-4.  The  values  of 
Singh  et  al . (1976)  come  from  Fig.  C-2  which  shows  that  most 
U.S.  (and  in  fact  global)  emissions  of  C Cl^  occurred  in  the 
19^03;  in  fact,  the  integrated  burden  of  atmospheric  injections 
of  carbon  tetrachloride  from  the  start  of  production  to  1955  is 
1.33  x 1CK  kg,  or  2.07  x 10^  Cl  atoms,  which  corresponds  to  16 
years  of  injection  at  the  1973  injection  rate  of  8l  x 10  kg. 

Thus , if  the  atmospheric  lifetime  of  carbon  tetrachloride  is 

* 

Dr.  F.  A.  Bower  (duPont)  comments  that  (1)  the  data  quoted  by 
Borisenkov  and  Kazakov  apply  to  the  U.S  S.R.  only,  while  the 
Grant  estimates  apply  to  the  whole  Eastern  Bloc,  including  the 
U.S.S.R.,  East  Germany,  Poland,  Czechoslovakia,  Rumania,  Bul- 
garia, and  the  People's  Republic  of  China,  (2)  these  data  will 
presumably  give  rise  to  slight  changes  in  future  estimates  of 
A.  Grant  and  Co.,  and  that,  (3)  the  ratio  of  FC-11  to  FC-12  as 
estimated  by  Grant  does  seem  irregular  when  compared  with 
Borisenkov  and  Kazakov.  Item  (3)  needs  clarification,  even 
though  its  effect  on  total  worldwide  emissions  is  not  likely 
to  be  significant. 
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TABLE  C-4 . GLOBAL  PRODUCTION  AND  RELEASE  RATES  OF 
CARBON  TETRACHLORIDE  THROUGH  1976  (in  106kg/yr) 
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Release  rates  from  U.S.,  W.  Europe,  and  Japan  have  been  read  from  Fig.  1 of  Singh  et  al . 
(1976),  added  and  multiplied  by  1.1  to  account  for  qlobal  production  (see  Fig.  2 of  Singh 
et  al.,  1976).  The  figures  may  be  interpolated. 


TOTAL  U.S.  DISPERSIVE  USAGE  (CURVE  4) 
TOTAL  U.S.  EMISSIONS  (CURVE  5) 

TOTAL  JAPANESE  EMISSIONS  (CURVE  6) 
TOTAL  EUROPEAN  EMISSIONS  (CURVE  7) 


CURVE  ! / 

i 

\/ 

CURVE  3 

, \ 1 

CURVE  4 / 

\' 

\ 

/ 

i 

i 

V . -'••'’’■CURVE  6 

CURVE  5 

> 

CURVE  7 
CURVE  2 \ 

\ l 

v 

Source:  Singh  et  al . , 1976.  Copyright  1976  by  the  American 
Association  for  the  Advancement  of  Science. 

FIGURE  C-2.  Production  and  Emission  Trends  for  CC 1 ^ 

Note:  Singh  recommends  using  1.1  times  (U.S.  + European  + 
Japanese  emissions;  for  Global  Emissions. 


longer  than  this,  the  principal  effects  arise  from  the  existing 
burden  rather  than  from  the  current  injections. 


The  other  estimates  come  from  the  sources  indicated  in  Table 
C-4 . With  the  exception  of  the  estimates  of  Little  (1975),  which 
are  much  lower  than  the  others,  they  agree  quite  well.  All  work- 
ers agree  that  the  principal  effect  comes  from  past  rather  than 
current  releases.  The  problem  is  that  the  pattern  of  use  of  car- 
bon tetrachloride  has  changed  dramatically  since  its  toxicity 
was  recognized.  In  the  1940s  and  1950s  it  was  largely  used  dis- 
persively  (as  cleaning  fluid,  in  fire  extinguishers,  etc.),  but 
now  it  is  overwhelmingly  used  in  the  production  of  fluorocarbons. 
Thus,  the  difference  between  assuming  3 percent  or  5 percent  re- 
lease of  total  production  makes  a very  large  difference  in  the 
estimate.  Further,  there  are  also  uncertainties  in  production, 
especially  outside  the  U.S.A.  [perhaps  + 20  percent;  see  also 
Altshuller  (1976)]. 

The  suggested  estimate,  listed  in  Table  C-5  and  in  Fig.  C-3, 
Is  simply  the  arithmetic  mean  of  all  the  different  values  re- 
corded in  Table  C-4. 

Production  and  release  figures  for  methyl  chloroform  are 
shown  in  Table  C-6,  listing  release  figures  from  Dow  Chemical 
Co.  (Neely  and  Plonka,  1978),  from  Little  (1975),  and  from  Mc- 
Connell and  Schiff  (1978).  The  differences  between  the  various 
estimates  are  not  as  large  as  in  the  case  of  carbon  tetrachlo- 
ride. Again,  the  recommended  estimate  for  releases  is  the  arith- 
metic mean  of  all  the  values  recorded  for  that  year  in  Table  C-6. 
This  is  shown  in  Table  C-6  and  Fig.  C-3. 

For  the  other  chlorocarbons  listed  in  Table  C-2,  the  emis- 
sion rates  are  given  in  Table  C-5,  using  the  data  of  A.  D.  Little 
(1975)  and  the  procedure  of  Annex  A;  in  Fig.  C-4  the  principal 
anthropogenic  global  emissions  in  this  category  are  shown  as  a 

function  of  time. 
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Note:  Data  for  carbon  tetrachloride  come  from  Table  C-4,  and  for  methyl  chloroform  from  Table  C-6.  In  each  case, 
we  list  here  the  arithmetic  mean  of  all  the  release  fiyures  quoted  in  the  approximate  table.  All  the 
other  figures  come  from  Little  (1975)  using  the  analysis  of  Annex  A to  this  Appendix. 


1975 


TABLE  C-6.  GLOBAL  PRODUCTION  AND  RELEASE  RATES  OF  METHYL  CHLOROFORM 

(in  10^  kg/yr) 


Year 

Production 

Dow 

Chemical9 

A.D.  Little,  Inc.*3 

McConnel 1 
and  Schiffc 

1951 

0.1 

0.1 

1952 

0.2 

0.2 

- 

1953 

1.0 

1.0 

- 

1954 

3 

3 

- 

1955 

8 

8 

_ 

1956 

14 

12 

- 

1957 

21 

20 

- 

1958 

22 

21 

- 

1959 

32 

30 

- 

1960 

38 

36 

_ 

1961 

40 

38 

- 

_ 

1962 

59 

56 

- 

- 

1963 

58 

51 

- 

- 

1964 

64 

57 

- 

- 

1965 

83 

73 

_ 

_ 

1966 

116 

109 

144 

147 

1967 

146 

131 

- 

- 

1968 

158 

145 

177 

- 

1969 

168 

148 

192 

- 

1970 

181 

155 

217 

_ 

1971 

191 

167 

222 

253 

1972 

267 

230 

261 

- 

1973 

350 

340 

324 

396 ; 407  ; 46  6 

1974 

390 

363 

349 

- 

1975 

374 

365 

272 

_ 

1976 

438 

416 

373 

413 

Cumulative: 

3222 

2975 

a.  Neely  and  Plonka  (1978),  Table  1, 

estimated  by  the  Corporate 

Product 

Department  of  Dow  Chemical  Company 

b.  Little 

(1975),  using 

the  procedure 

of  Annex  A of  this 

Appendi x 

c.  McConnell  and  Schiff 

(1978)  quote  various  estimated  production 

figures 

for  some  years;  I convert 

to  release  by  multiplying 

by  the 

ratio  release/production  for  that 

year  from  Neely  and 

Plonka 

(1977) 
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C . 3 LATITUDINAL  DISTRIBUTION  OF  EMISSIONS 


For  a two-dimensional  computer  code,  it  is  necessary  to 
estimate  the  latitude  distribution  of  anthropogenic  halocarbon 
emissions.  The  A.  D.  Little  (1975)  and  A.  Grant  (1977)  reports 
provide  the  following  data,  which  are  not  strictly  comparable; 
these  are  used  here  to  provide  an  initial  best  estimate. 


1.  Source:  A.  D.  Little,  p.  II-9 

World  fluorocarbon  production,  1973 


Region 


Percent  of 
Global  Production 


U.S.A. 

Europe 

Japan 

Other 


Source:  Bureau  of  Domestic  Commerce  Estimate 

(from  U.S.  International  Trade  Commission) 


World  fluorocarbon  production,  1973 


Region 


Percent  of 
Global  Production 


U.S.A. 

Europe 

Japan 

Communist  Bloc 
Other 


Source:  A.  Grant,  Schedule  7 


Total  production  of  fluorocarbons,  ratio 
southern/northern  hemisphere 


FC-11 


0.035 


FC-12 


0.042 


0.039 


FC-22 


0.035 


Source 


A.  Grant,  Schedules  1,  11,  12,  total  produc- 
tion, 1931  to  1976,  fraction  of  global  pro- 
duction 


Substance 


FC-11 


FC-12 


FC-22 


"India  and  Argentina" 

0.00093  ) 

> 0.0015 

0.0018  ) 

0.0054 


"Eastern  Bloc' 


0 .048 


0.043  | 

0.0062 


0.045 


J 
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5.  Source:  Table  C-7  [from  A.  Grant  (1977),  Schedule  1 
for  listing  of  primary  countries,  Schedules 
11  and  12  for  others]  lists  the  countries 
where  fluorocarbons  are  produ-ced,  together 
with  the  effective  range  of  latitudes 
assigned  to  each  country. 


TABLE  C-7.  LISTING  OF  COUNTRIES  PRODUCING  FLUOROCARBONS, 
INCLUDING  LATITUDE  RANGE  TO  BE  CONSIDERED 

NO.  OF  COUNTRIES 

OR  SUBSIDIARIES  LATITUDE  RANGE 

COUNTRY  IN  TABULATION*  CONSIDERED 


Argentina 

1 

35-41 

S 

Australia 

2 

28-39 

S 

Brazi  1 

2 

5-33 

s 

Canada 

1 

42-48 

N 

France 

2 

42-50 

N 

Great  Britain 

2 

50-55 

N 

India 

- 

8-32 

N 

Italy 

1 

39-45 

N 

Japan 

4 

33-45 

N 

Mexico 

1 

15-20 

N 

Netherl ands 

2 

51-53 

N 

Spain 

2 

35-43 

N 

U.S.A. 

6 

30-48 

N 

U.S.S.R. /COMECON 

- 

40-60 

N 

West  Germany 

2 

47-54 

N 

Grant  ( 1 977  ) , Schedul e 1 

From  the  data  quoted  in  items  1 through  4 above,  and  in- 
cluding Table  C-7,  it  is  clear  that  the  primary  production  of 
fluorocarbons  is  in  the  latitude  range  of  30  to  55°  N,  which  in- 
cludes the  U.S.A.,  Europe,  Japan,  and  most  of  U ,S .S .R ,/COMECOM. 
Bearing  in  mind  that  the  category  "other,  7 percent"  in  item  1 
includes  COMECON  at  about  4 percent — see  item  3,  this  would  leave 
2.5  percent  unaccounted  for;  and  this  may  be  compared  with  item 
3,  4 percent  in  the  Southern  Hemisphere. 


\ 
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The  difference;  between  these  items  should  riot  be  considered 
'nil' lean'.  . Taking  into  account  the  difference  between  produc- 
■ I m,  sal  i ■ a.  imi  1 ' >n,  and  emissl  -n  figures,  he  lifference  be— 
' : i . ■ %v  ed  product!  n and  product!  n rut  e (;,•■:•  annum ) , 

and  the  difference  between  the  different  chemicals  listed  in 

* 


Tabl  • 

;-l,  I r ;commend 

the  follow:; 

i ng 

two  options  for 

L i t 1 0/ 11 

. 

le  lis  t rl  but  Lon  o f 

emissions 

of 

all  halocarbons 

Option  I 

Option  II 

50  to  60°  II 

5$ 

30 

to  55°  N 96% 

30  to  50°  N 

90% 

30 

to  40°  S 4% 

30°  II  to  30°  S 

2 % 

30  to  40°  S 

n al 

J h 

where  the  emission  should  be  distributed  uniformly  within  each 
r<  Levant  latitude  band.  Option  I is  recommended  as  primary, 
with  Option  II  as  an  alternative. 


C . 4 BROMOCARBONS 

It  has  been  suggested  by  Wofsy  et  al.  (1975)  that  bromine 
injections  into  the  stratosphere  can  also  reduce  the  concentra- 
f strat  S{  heric  ozone.  For  the  present  catalog,  this  is 
unlikely  to  be  important  because  the  major  scarce  of  strato- 
spheric bromine  is  presumably  natural--cf.  Wofsy  et  al.  (1975)  — 
and  the  relevant  anthropogenic  production  was  relatively  small 
during  the  1955  to  1975  time  period  and  has  not  changed  very 
much . 

The  largest  single  use  of  bromine — representing  perhaps  60 
percent  of  global  production — is  ethylene  dibromide  (1,  2 dl- 
bromoethane ) used  as  a scavenger  in  leaded  gasoline.  The  bro- 
mine (or  equivalently  chlorine)  is  added  so  that  the  lead  (from 
3 — — — 

For  chlorccarbons , there  are  estimates  of  the  breakdown  of  1973 
Is  i ns  ini  the  jateg  ri<  s "".f.,"  "Europe"  and  "Other"— 

. Ltt2  . • "• . : . III-!  ) , but  these  do  not  differ 
irasti  . from  th  present  estimates  for  fluorocarbons . 
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tetraethyl  lead  used  as  anti-knock  agent)  will  be  scavenged  as 
PbBr„,  etc.,  rather  than  deposited  inside  the  enmlne.  Like  most 

<L 

other  bromides,  this  lead  bromide  will  be  exhausted  as  fine  par- 
ticles which  are  mostly  washed  out  by  precipitation  or  else  fall 

to  the  ground  by  sedimentation,  rather  than  being  transported 

# 

out  of  the  "planetary  boundary  layer"  to  enter  the  "free 
troposphere . " 

Approximately  2 percent  of  man-made  ethylene  dibromide,  as 
well  as  mostNjnan-made  methyl  bromide,  are  used  as  soil  and  grain 
fumigants.  Soirie  fraction  of  these  materials  enters  the  atmos- 
phere unchanged  arid  so  supplements  the  natural  sources  of  bro- 
mine from  seaweed,  etc.,  in  the  free  atmosphere.  In  the  atmos- 
phere, ethylene  dibromide  and  methyl  bromide  will  react  rather 
slowly  witn  OH,  presumably  with  lifetimes  comparable  to  those  of 
corresponding  chlorocarbons  (see  Table  C-l , where  lifetimes  of 
order  0.2  to  0 . years  are  listed).  Thus,  the  relative  contri- 
bution of  these  bromocarbor.s  to  stratospheric  bromine  will  be  com- 
parable to  the  relative  contribution  to  stratospheric  chlorine  of 
materials  like  ethylene  dichloride  end  methyl  chloride,  rather 
than  that  of  materials  like  FC.-11  and  FC-12 . 

Table  C-8  lists  industrial  production  figures  for  methyl  bro- 
mide and  ethylene  dibromide  for  the  time  period  1955  to  1976,  and 
Table  C-9  lists  my  best  estimate  of  sources  of  bromine  that  enter 
the  free  atmosphere.  It  is  clear  that  the  total  man-made  source 
of  bromine  is  not  only  relatively  small,  but  also  its  variation 
over  the  1955  to  1975  time  period  is  relatively  small  within  any 
reasonable  error  or  variation  in  the  natural  source.  At  present, 
while  a constant  ground-level  source  of  (950  + 250)  x 10^  kg 

The  lowest  1-km  portion  of  the  earth's  atmosphere  in  which  large 
day-night  temperature  variations  occur,  and  in  particular,  rapid 
transfer  of  water  occurs  between  the  atmosphere  and  land/sea. 

Most  air  pollution  processes  take  place  in  the  boundary  layer, 
and  only  a small  fraction  of  air  pollutants  ever  rise  above  this 
region  to  enter  the  "free  troposphere"  which  lies  above  and  out- 
side the  boundary  layer, 

C-  20 


« 


- mi 


Br./yr  is  appropriate  in  a model  of  atmospheric  bromine,  includ- 
ing its  potential  effects  on  ozone,  no  significant  variation  in 
this  source  appears  to  have  occurred  during  the  interval  of 
interest  here,  1955  to  1975. 


TABLE  C - 8 . WORLD  INDUSTRIAL  PRODUCTION  OF 
METHYL  BROMIDE  AND  ETHYLENE  DIBROMIDE 


Annual  Production  (in  10^  kg) 


Year 

Methyl  Bromide 

Ethylene  Dibr 

1955 

5 

_ 

1960 

5 

120* 

1965 

8 

- 

1970 

13 

190 

1972 

16 

205 

1974 

13 

190 

1976 

16 

160** 

*Ear1y  1960s 

* * 

The  decrease  Is  genuine,  associated  with  In- 
creasing U.S.  use  of  unleaded  gasoline. 


Note:  E.  Bauer's  estimates,  based  on  data  from 

the  Bureau  of  Mines;  the  U.S.  International 
Trade  Commission;  and  Kirk-Othmer  (1965), 

Vol . 3,  pp.  763-772.  Most  primary  data 
listed  are  for  U.S.  production,  and  have 
been  scaled  on  the  assumption  that  U.S.  pro- 
duction represented  75  percent  of  worldwide 
production  in  1960,  70  percent  in  1973. 

Errors  of  at  least  10  to  20  percent  are  prob- 
able, l.e.,  not  all  figures  listed  are  sig- 
nificant. "Ethylene  dibromide"  is  1,  2,  di- 
bromoethane . 


TABLE  C-9.  POTENTIAL  ATMOSPHERIC  SOURCES  OF 
BROMINE  IN  1975  (in  106  kg  Br./yr) 


Source 


Wofsy  et  al . (1975)  Present  Estimate 


Natural  Source  900  900a 

Methyl  Bromide  100  12  to  14b 

Ethylene  Dibromlde  1 00c  4 to  5b 

aF rom  Wofsy  et  al . ( 1 975) 

bDr.  Wm.  Hunt,  private  communication 

cWofsy  et  al  . ( 1 975  ) assume  that  50  percent  of  bromine  from 
ethylene  dibromide  can  enter  the  lowest  atmosphere  In  Inor- 
ganic form.  However,  Dr.  Wofsy  in  a private  communication 
stated  that  most  of  it  will  be  washed  out,  and  a negligible 
fraction  (<  1%)  will  enter  the  stratosphere. 


A 
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ANNEX  A TO  APPENDIX  C 


PRESENT  ANALYSIS  OF  DATA  FROM  LITTLE  (1975) 

TO  OBTAIN  ANNUAL  RELEASE  RATES* 

Little  (1975)  quoted  estimated  world  chlorocarbon  emissions 
for  1973  (Table  CA-1).  They  also  list  U.S.  production  from 
1960-74--see  pp . 11-13,  -14,  op.  ait.  In  general,  U.S.  produc- 
tion rates  are  about  half  worldwide  production  rates--see  Table 
CA-1;  Little,  for  1973,  estimates  annual  worldwide  ratios  of 
emission/production  to  be  equal  to  annual  U.S.  emission/produc- 
tion ratios — see  p.  III-27,  op.  ait. 

Here  I supplement  U.S.  production  figures  for  1955  to  i960 

from  Kirk-Othmer  (1965),  Vol . 5,  pp . 107-201,  and  for  1975  to 
1976  from  U.S.  International  Trade  Commission  figures  (Reports 
804  and  833,  respectively). 

I also  assume,  (1)  that  the  ratio  of  annual  world  produc- 
tion rate/U.S.  production  rate  remains  constant  from  1955  to 

1974  (at  the  ratios  given  in  Table  CA-1  for  each  chemical  spe- 

* * 

cies),  and  (2)  that  the  ratio  of  annual  world  emission/annual 
world  production  remains  constant  from  1955  to  1974  at  the  1973 
ratios  estimated  by  Little  (1975). 


This  is  not  necessarily  a satisfactory  analysis,  but  it  is 
objective  in  terms  of  the  limited  data  base. 

# to 

Galbally  (1976)  estimates  the  U.S. /Global  production  ratio 
for  CCI4  as  90  percent  in  1955  and  50  percent  in  1972;  the 
Little  (1975)  figure  for  1973  is  50  percent — see  Table  CA-1 


C-23 


TABLE  CA-1 . ESTIMATED  WORLD  CHLOROCARBON  EMISSIONS,  1973 


Regional  Percent  of 

1otald 

World  Emissions 
(106  kq) 

Chlorocarbon 

United  States 

Europe 

Other 

Carbon  tetrachloride 

50 

40 

10 

41  .7 

Chloroform 

50 

40 

10 

12.4 

Ethyl  chloride 

55 

30 

15 

14.6 

Ethylene  dichloride 

35 

40 

25 

565.2 

Methyl  chloride 

60 

35 

5 

7.9 

Methyl  chloroform 

60 

30 

10 

324.2 

Methylene  chloride 

55 

30 

15 

346.4 

Perchloroethylene 

45 

45 

10 

609.0 

Trichloroethylene 

30 

55 

15 

648.3 

Vinyl  chloride 

25 

40 

35 

351  .6 

aRounded  to  nearest  5 percent. 

Source:  A!  D.  Little  (1975),  p.  Ill -29 


For  carbon  tetrachloride  and  methyl  chloroform  there  exist 
other  estimates — see  Section  C.2,  especially  Tables  C— 4 and  C-6. 
In  Figs.  C-3  and  C-4  I compare  different  estimates  for  emis- 
sion rates  of  these  two  species,  and  conclude  that  the  various 
estimates  agree  quite  well  for  methyl  chloroform,  but  not  for 
carbon  tetrachloride.  The  disagreement  for  carbon  tetra- 
chloride is  unlikely  to  be  very  significant,  since  it  is  the 
burden  of  pre-1955  injections  rather  than  the  annual  injections 
thereafter  that  are  probably  most  important  for  the  present 
atmospheric  loading  of  this  material,  but  the  deviations  be- 
tween these  different  estimates  does  raise  a warning  flap  about 
possible  errors  in  the  present  estimates  of  atmospheric  emis- 
sion rates. 
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APPENDIX  D 


STRATOSPHERIC  INJECTIONS  DUE  TO 
VOLCANIC  ERUPTIONS,  1955  TO  1976 

D. 1 DATA  AND  MODEL 

A listing  of  major  volcanic  eruptions  from  1955  to  1976  is 
given  in  Table  D-l.  The  list  is  not  complete — Dr . S.  Self 
(Dartmouth  College)  is  preparing  a much  more  comprehens ive  list- 
ing— see,  e.g.,  Barbera  and  Self  (1978) — but  it  probably  in- 
cludes all  eruptions  that  sent  major  amounts  of  material  int 
4 he  stratosphere.  Table  D-S  lists  the  suggested  stratospheric 
injection  due  to  the  major  eruptions  of  Mt . Agung  in  1963-  We 
show  the  estimates  of  Cadle  et  ai.  (1976,  1977)  based,  ultimately* 
on  some  unpublished  estimates  of  Lamb  (see  Mitchell,  1970),  ad- 
justed to  agree  with  the  mass  of  fine  dust  estimated  by  Mossop 
(19611),  and  also  those  of  Self  (private  communication)  which  are 
based  solely  on  volcanic  parameters,  in  particular  on  the  report 
of  Zen  and  Hadlkusumo  (1964).  In  view  of  the  very  limited  data 
base  and  the  different  methods  used,  the  relatively  good  agree- 
ment of  these  two  data  sets  must  be  considered  rather  satisfactory. 

Regarding  the  height  of  injection,  we  may  assume  that  for  a 
given  eruption  the  material  is  distributed  uniformly  between  the 
tropopause  and  an  observed  maximum  height.  Values  for  these 
heights  are  listed  in  Table  D-3  for  the  largest  volcanic  erup- 
tions; once  again,  we  may  note  the  very  limited  data  base. 

Reference  should  be  made  to  Rig.  4 of  Lamb  (1970),  to  Friend 
(1972),  and  especially  to  Wilson  et  al . (1978). 
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TABLE  D-l.  MAJOR  EXPLOSIVE  VOLCANIC  ERUPTIONS,  1955-76 


Stratospheric 


Date 

Volcano 

Position 

Seven  ty  " 

♦ b 

I nput 

Jun  1955 

Puyehue  (and  Rancoj,  Chile 

40-41°  S,  72°  W 

2 

3 

(0.3  km  tephraj 

No 

30  Mar  1956 

Bezymianny,  Kamchatka 

56°  N , 160°  E 

1 

(1  kmJ  tephra, 

Yes 

21-24  May  1960 

Puyehue  (etc),  Chile 

39-45°  S,  72-73°  W 

3 

Yes 

1962 

Tokati 

43°  N,  143°  E 

Probable , 
but  smal 1 

17  Mar  and  16  May  1963 

Agung , Bali 

8.5°  S,  115.5°  E 

1- 

2 

Yes 

12  Nov  1964 

Sheveluch,  Kamchatka 

56.4°  N,  161°  E 

1 

a 

(1  km  tephra; 

?? 

28  Sep  1965 

Taal,  Luzon 

14°  N,  121°  E 

3 

Yes 

1966 

Redoubt,  Alaska 

60.5°  N,  1 53°  W 

Probabl e 

1966 

Kelud 

8°  S,  112°  E 

Probable, 
but  smal 1 

12  Aug  1966 

Awu,  Celebes 

3.5°  N,  125.5°  E 

2 

( ■ 3 kmJ  dust) 

Yes 

4-7  Dec  1967 

Deception  Island 

63°  S,  60.5°  W 

3 

No 

1967-68 

Redoubt  and  Trident 

No 

11-12  Jul  1968 

Fernandina  Is.,  Galapagos 

0.5°  S,  72°  W 

2 

( K 1 kinJ  ?) 

Yes 

5-6  May  1970 

Hekla,  Iceland 

64°  N,  19.5°  W 

Yes 

1970 

Deception  Island 

No 

1971 

Mt.  Hudson,  Chile 

46°  S,  73°  W 

Probable 

1972 

Alaid,  Kurile  Islands 

51°  N,  155°  E 

Probable 

23  dan-16  Feb  1973 

Helgafell,  Iceland 

63.5°  N,  20.3°  W 

No 

1974 

ApiSian 

8°  S,  119°  E 

Probabl e 

10-23  Oct  1974 

Fuego,  Guatemala 

14.5°  N,  91°  W 

2* 

Yes 

1975 

Plosky  Tolbuchik,  Kamchatka 

56°  N,  160°  E 

Probabl e 
but  small 

23  Jan  1976 

Augustine,  Alaska 

59.4°  N,  153.4°  W 

3 

Yes,  but 
small 

A qualitative  measure  of  the  overall  magnitude  of  the  eruption  (not  necessarily  its  stratospheric 
injection!),  see  Sapper  (1927)  and  Lamb  (1970).  1 is  biggest,  2 is  an  order  of  magnitude  smaller, 

3 is  an  order  of  magnitude  smaller  yet,  in  terms  of  either  mass  or  volume  of  material  moved.  Note 
that  no  successful  scaling  has  been  achieved  in  that  the  volume  of  material  injected,  the  enerey 
released,  and  the  amount  and  chemical  composition  of  gases  emitted  are  not  simply  related.  It  is 
likely  that  only  studies  of  individual  eruptions  will  have  success  in  predicting  their  stratospheric 
injection,  in  terms  of  the  total  mass  and  composition. 

^From  S.  Self,  private  communication,  May  1978. 
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TABLE  D - 2 STRATOSPHERIC  SOURCE  FUNCTION  OF  THE 
AGUNG  ERUPTION  CLOUD 


Mass,  kg 


Estimate  of 


Estimate  of 
Cadle  et  al . 
(1977) 

S.  Self  (private 
communication, 

May  1978) 

Suggested  h 
Value 

c 

Total  erupted  material 

5 x 1012 

1.4  x 1012 

12 

2.6  x 10 

r-  . d 

Fine  ash 

1 x 1010 

9 x 109 

9.5  x 109 

Total  gas 

- 

2.9  x 1010 

- 

Water  vapor  e 

1.2  x 1011 

2.2  x 1010 

5.1  x 1010 

Total  S as  H^SO^  f 

1.8  x 1010 

2.1  x 109 

6.1  x 109 

HC1  9 

1.2  x 109 

7.3  x 108 

9.4  x 108 

d Based  on  17  March  and  16  May  eruptions. 

^ Geometric  mean  of  the  Cadle  et  al . and  Self  estimates. 
c Most  of  this  material  does  not  enter  the  stratosphere. 

d Cadle  et  al (1977)  do  not  specify  the  particle  size;  Self  refers  to  particles 
< 2 pm  diameter. 

e 12 

The  mean  stratospheric  water  vapor  loading  is  n,  2 x 10  kg,  see  Appendix  A. 

f 9 34 

5.1  x 10  kg  H,S0.  contain  3.3  x 10  S atoms  so  that  presumably  (3-6)  x 
34  ^ ^ 

10  OH  radicals  are  consumed  in  the  formation  of  sulfuric  acid. 

n ft  *3  7 

9 9.4  x 10  kg  HC1  contain  9 x 10  molecules. 


D-5 


In  Table  D-3,  we  also  list  recommended  source  functions  for 
the  major  eruptions  between  1955  and  1975.  The  follow!  lg  com- 
ments should  be  noted: 

• Uncertainties  in  absolute  injected  amounts  by  a factor 
2 to  5 (up  or  down)  exist,  but  the  ratios  quoted  in 
Table  D-3  may  have  a smaller  uncertainty. 

• The  S/Cl  ratio  differs  drastically  between  different 
types  of  eruptions.  For  Icelandic  and  Hawaiian  erup- 
tions (effusive,  not  explosive)  the  ratio  is  of  order 
10  to  20:1,  while  for  Central  American  volcanoes,  which 
are  more  explosive,  the  ratio  is  of  order  1:1.  [See 
Anderson  (1975),  Stoiber  and  Rose  (1973)  ] • 

• The  average  gas  composition  used  by  S.  Self  is  75  percent 
H^O,  15  percent  C0o,  5 percent  S0,2 , 2.5  percent  HC1, 

2.5  percent  others,  but  evidently  large  variations  from 
this  mean  composition  occur. 

• In  particular,  for  Bezymianny  and  Sheveluch  there  is 
very  little  data  available. 

D . 2 POSSIBLE  INFLUENCES  ON  STRATOSPHERIC  OZONE 

The  stratospheric  injections  due  to  the  Agung  eruption,  as 
listed  in  Table  D-3,  should  be  compared  with  the  'normal"  strato- 
spheric loading  of  the  relevant  materials.  Thia  is  done  now  for 
each  component  in  turn.  (See  also  Stolarski  and  Butler,  1977, 
Turn  et  al.,  1978.) 

D . 2 . 1 Aerosols/Sulfuric  Acid 

In  the  lower  stratosphere  there  is  normally  a region  of  fine 
particles,  the  density  and  composition  of  which  is  quite  variable. 
Typical  loadings  in  nonvolcanically  perturbed  times  (e.g.,  1972 
to  1973,  prior  to  the  eruption  of  Agung)  are  of  order  of  one  par- 

> O, 

tide  of  radius  'v  0.15  um  per  cm",  with  a vertical  extent  of  5 km 
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centered  about  18  to  20  km,  the  so-called  "Junge  layer"  [cf. 

q p 

Cadle  et  al.,  1 97 5 J - This  corresponds  to  5 x 10  particles/m  - 

Q 

column,  or  a global  loading  of  order  10°  kg  (if  each  particle 
has  radius  0.15  um  and  a density  of  2 to  3 gm/cm  ) . 

Stratospheric  aerosols  are  generally  considered  to  be  made 
up  of  a mixture  of  impure  sulfuric  acid  (including  some  ammonium 
sulfate)  and  of  silicates,  i.e.,  fine  volcanic  ash.  The  rela- 
tive amount  of  these  materials  found  by  different  observers  at 
different  times  and  places  varies  significantly.  Table  D-3  lists 
a current  estimate  of  the  stratospheric  loading  due  to  the  Agung 
eruption,  with  comparable  amounts  of  sulfuric  acid  and  silicates, 
and  a total  mass  injection  of  order  several  times  1010  kg.  This 
is  consistent  with  the  current  concept  that  explosive  volcanic 
eruptions  are  the  major  source  of  the  Junge  aerosol  layer. 

Regarding  the  effect  of  this  injection  on  the  stratosphere, 
it  is  generally  considered  that  stratospheric  aerosols  of  vol- 
canic origin  affect  the  climate  [cf.,  e.g.  Oliver  (1976)],  but 
it  is  possible  that  they  may  also  affect  the  stratospheric  photo- 
chemistry; thus,  Cadle,  Crutzen,  and  Ehhalt  (1975)  discuss  hetero- 
geneous reactions  that  may  take  place  on  the  stratospheric  aero- 
sols. Also,  if  the  sulfur  enters  the  stratosphere  as  SO^  which 
is  there  oxidized  to  sulfuric  acid  by  reaction  with  OH  radicals, 
locally  and  for  a time  of  order  a month  this  oxidation  of  S0„ 
may  also  provide  a significant  loss  mechanism  for  the  OH  radi- 
cals, and  could  thus  affect  the  local  ozone  concentration. 

D.2.2  Hydrochloric  Acid 

The  total  stratospheric  loading  of  C1X  = Cl  + CIO  + HC1  + 

34  8 

CIO  NO^  is  approximately  1 x 10  molecules  worldwide  (6  x 10 

kg  as  HC1) — see  Appendix  A.  Thus  we  see  that  the  stratospheric 
injection  of  HC1  due  to  the  present  model  of  Agung  is  almost  as 
large  as  the  global  loading  of  C1X.  This  might  be  expected  to 
show  an  effect,  either  by  reducing  the  ozone  through  the  cata- 
lytic C1X-0-.  cycle  or  else  by  trapping  H0y  and  N0y  as  H0C1  and 


D.2.3  Water 


The  normal  strain. :nhori  • Lou  3Jn  • < f water  is  juot-..-d  as 
12 

[1.6  to  2.2)  x 10  kg — see  Appendix  - - t hal  th<  eruption  : ' 

Agung  presumably  led  to  a 5 to  7 uerc-.-nl  enhancement  of  global 
stratospheric  water  vapor. 

We  see  that  at  the  ii  j e :tion  ev  1 1 by  id! 

al.  (197<  , L977) » a maj < r v : inic  rupt  i i,  sucl  a:  that  f 
Agung,  would  be  ex:  :ted  t produc  si  fi  ■ • • ffe  :ts  n 1 he 

stratosphere,  although  the  specifi  • ff<  :ts  ar  sul  Jecl  1 mucl 
uncertainty.  Angel!  and  K rsh  v r (]  o'  fii  . ai  tncreas:  in 
total  ozone — as  well  as  a change  in  the  quaslbiennlal  oscilla- 
tion — occuring  shortly  after  the  eruption  of  Agung.  A number 
of  workers  have  presented  annual  mean  volcanic  injections  [e.g.. 
Cicerone  (1975),  Ryan  and  MukherJ ee  (1975)],  which  is  appr  j riat: 
in  asking  for  sources  and  sinks  of  JlX,  aer  sols,  el  .• . , 1 it  vi- 
dently  not  for  the  present  purpose. 
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TIME-DEPENDENT  PRODUCTION  OF  NO  AND  OH  IN  THE 
STRATOSPHERE  RESULTING  FROM  COSMIC  RAY  IONIZATION 
IN  THE  1955-1975  TIME  FRAME 

E.l  INTRODUCTION  AND  SUMMARY 

Ruderman  and  Chamberlain  (1975)  and  Nicolet  (1975),  enlarg- 
ing on  an  earlier  calculation  of  Warneck  (1972),  pointed  out  that 
galactic  cosmic  rays  (GCR)  produce  ionization  in  the  stratosphere 
which  leads  to  increased  NO  production  as  a result  of  ion  chem- 
istry, with  subsequent  effects  on  stratospheric  ozone.  The  GCR 
intensity  is  modulated  by  the  solar  wind  and  thus  varies  with  the 
11-year  sunspot  cycle,  leading  to  an  11-year  oscillation  in  NO 
(=  NO  + NO^)  production  with  its  maximum  at  sunspot  minimum. 
Crutzen  et  al . (1975) — see  also  Heath  et  al . ( 1977 ) --pointed  out 
that  some  large  solar  flares  produce  bursts  of  high-energy  pro- 
tons, the  solar  cosmic  rays;  these  "Solar  Proton  Events"  (SPE) 
also  produce  ionization  in  the  stratosphere/mesosphere  that  would 
be  expected  to  affect  NO  and  hence  stratospheric  ozone.  SPEs 
occur  mainly  near  (but  not  at)  sunspot  maximum.  Both  GCRs  and 
SPEs  are  modulated  by  the  geomagnetic  field  and  thus  their 
effects  are  maximized  at  high  geomagnetic  latitudes,  where  charg- 
ed particles  of  relatively  low  energy,  momentum,  or  "rigidity" 
can  penetrate  most  easily  into  the  earth's  atmosphere. 

There  are  also  other  sources  of  atmospheric  ionization  which 
may  exhibit  some  time-dependence,  such  as  X-rays,  EUV  (in  partic- 
ular Lyman-alpha,  121.6  nm)  , auroral  electrons,  etc.,  which  all 
lead  to  ionization  at  altitudes  above  60  to  80  km.  At  these 
altitudes,  there  is  very  little  ozone,  and  the  lifetime  of  NO  is 
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so  short  that  all  the  NO  produced  as  a result  of  ionization  will 
be  destroyed  photochemically  before  it  can  move  down  to  altitudes 
where  there  is  a significant  amount  of  ozone  (mainly  below  35  km). 

Two  other  effects  which  have  been  considered  are  the  follow- 


a.  relativistic  electron  precipitation  (REP)  events  (Thorne, 
1977,  1978),  which  can  be  shown  to  be  small  in  the  altitude  range 
of  interest  here,  below  A5  km--see  Figs.  E-7  and  EB-2 . 

b.  high  energy  OCRs  (above  0.1  to  0.3  GeV)  produce  chlo- 

rine isotopes  36,  38,  39  in  the  atmosphere  through  nuclear  re- 
actions with  atmospheric  argon  (see  Lai  and  Suess,  1968).  While 
these  Chlorine  isotopes  do  undergo  chemical  reactions  (Rowland, 
1978  ) arid  thus  can  affect  the  ozone,  the  total  number  of  atoms 
produced  is  so  small  (£,  10  times  the  number  of  HOY  or  NC) 

molecules  produced  by  GCRs)  that  its  effect  on  ozone  will  be 
negligible . 

Here  we  summarize  all  currently  known  time-dependent  sources 
of  ionization  in  the  stratosphere  which,  by  producing  either 
oxides  of  nitrogen  (M0^  = NO  + NO.-,)  or  oxides  of  hydrogen  (H0y  = 

OH  + H 0P),  would  be  expected  to  have  produced  time-varying  re- 
ductions in  stratospheric  ozone  during  the  1955  to  1975  time 
period . 

The  conclusions  of  this  discussion  are  the  following: 

• One  ion  pair  produces,  on  the  average,  1.3  NO  molecules, 
and  2.0  OH  radicals  (see  Section  E.3  and  Fig.  E-5). 

• The  largest  and  principal  time-varying  effects  of  (solar 
and  galactic)  cosmic  ray  ionization  on  ozone  occurs  in  the 
polar  cap  regions  of  high  geomagnetic  latitude  A >,  60°  to 
65°,  which  are  shown  in  Fig.  E-2 . Also  see  Table  E-l  and 
Fig.  E-3  for  an  effective  representation  in  terms  of  geo- 
graphic latitude  9 for  use  in  a two-dimensional  photo- 
chemical/dynamic model  of  ozone. 


• At  high  latitudes,  most  of  the  ozone  lies  below  25  to  30 
krr  ( see  Fig . E-l ) . 

• Because  of  the  altitude  variation  of  photochemical  life- 
times for  NO  and  OH,  and  also  the  characteristic  times  for 
vertical  diffusion,  NO  produced  below  approximately  45  km 
can  affect  the  bulk  of  the  ozone,  but  because  of  its  rela- 
tively short  lifetime,  the  OH  radical  is  most  unlikely  tc 
diffuse  far  enough  (say,  >,  1 km)  to  affect  ozone  other 
than  locally  (see  Annex  A to  this  Appendix). 

• The  ionization  profile  due  to  GCRs  may  be  expressed  as  a 
function  of  both  time  (relative  to  the  year  YM  of  GCR 
maximum,  i.e.,  solar  sunspot  minimum)  and  of  latitude 
(where  we  distinguish  between  the  polar  cap  region,  de- 
fined in  terms  of  the  curves  C of  Fig.  E-2  (or  A > 60°) 

o 

and  the  rest  of  the  earth,  i.e.,  curve  "1-C  " (or  A < 60°) 

c 

as  follows : 

IONIZATION  PROFILE  DUE  TO  GCRs  IN  TERMS 
OF  EQUATIONS  OF  SECTION  E.4 

High  Geomagnetic  Latitude,  Low  Geomagnetic  Latitude, 


z > 20  km,  Eq  (9)  z > 31  km,  Eq  (12) 

15  km  < z < 20  km,  Eq  (10) 

10  km  < z < 15  km,  Eq  (11)  18  km  < z < 31  km,  Eq  (13) 


The  ionization  profile  due  to  SPEs  may  be  taken  as  uniform 
with  latitude  over  the  polar  caps  defined  in  terms  of  the 
inner  curves  Ci  of  Fig.  E-2,  or  equivalently,  for  geomag- 
netic latitudes  A >,  65°.  The  altitude  profile  of  ioniza- 
tion due  to  SPEs  is  shown  in  Section  E.5,  Fig.  E-7 . 


The  overall  magnitude  of  the  effect  of  GCRs  and  SPEs  may 
be  represented  as  follows: 


E-5 


. 


j 


At  solar  sunspot  minimum,  the  total  number  of  ion 
pairs  produced  by  GCRs  in  the  stratosphere  over 
both  polar  caps  (A  _>  00°),  as  well  as  the  total 
number  of  ion  pairs  produced  in  the  stratosphere 
over  lower  latitudes  (A  £ 60°)  each  exceed  the 
appropriate  number  of  ion  pairs  produced  at  sun- 
spot maximum  by  approximately  6 x 10^"  ion  pairs 
per  year  (see  Section  E.4). 

A large  SPE  produces  ^ 20  x 10^2  ion  pairs  (see 
Table  E-3) . 

During  the  1955  to  1975  time  frame,  there  have 
been  just  three  major  series  of  SPEs  that  should 
be  considered  from  the  standpoint  of  ozone  deple- 
tion, namely  those  of  May-July  1959,  November 
I960,  and  August  1972  (see  Table  E-3  and  Fig.  E-7 ) • 

The  Soviet  atmospheric  nuclear  test  series  of  Fall 
1961  (the  largest  ever,  corresponding  to  180  Mt 
total  yield)  injected  (70  to  270)  x 10^2  NO  mole- 
cules into  the  stratosphere. 

The  global  stratospheric  loading  of  NOy  (=  NO  + 

N02  + HNO3)  is  of  order  (4  to  15)  x 10^4  molecules. 

— The  global  stratospheric  loading  of  H0x  radicals 

(=  OH  + HC,)  is  of  order  (7  to  16)  x 1032  molecules, 

* 

but  with  very  large  uncertainties. 

E . 2 PENETRATION  OF  IONIZING  RADIATION  INTO  THE  EARTH'S 
ATMOSPHERE 

We  have  to  consider  X-rays  and  gamma-rays,  electrons  and 
heavy  particles,  principally  protons.  All  electromagnetic  ion- 
izing radiation  of  wavelengths  longer  than  0.3  nm  is  absorbed 
above  80  km,  and  the  overall  intensity  of  X-rays  of  wavelength 


See  Appendix  A for  a discussion  of  this  point. 
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less  than  0.3  nm  is  so  low  that  this  again  may  be  ignored  in 
the  present  context.  Regarding  particulate  radiation.  Fig.  E-l 
shows  the  effective  penetration  depth  of  electrons  and  protons 
for  vertical  incidence.  The  following  points  should  be  noted: 

• For  oblique  incidence,  the  actual  penetration  depth  is 
slightly  less  than  the  value  shown.  Thus,  the  penetra- 
tion depth  is  reduced  by  about  5 km  for  angles  of  inci- 
dence greater  than  60° , and  by  10  km  for  angles  of  inci- 
dence greater  than  80°  (see  Reid,  1974). 

• For  charged  particles,  the  energy  loss  due  to  ionization 
(average,  35  eV/ion  pair  in  air)  maximizes  during  the 
last  portion  (say,  1 to  2 scale  heights)  of  the  path  [cf. 
Fermi  (1950,  p.  33)]. 

• For  protons  at  energies  in  excess  of  100  MeV,  nuclear 
collisions  are  a significant  source  of  energy  loss 
through  the  production  of  "stars"  and  "showers."  The 
physics  here  is  significantly  different  [cf.  Fermi 
(1950,  p.  215)]. 

The  ionization  in  the  atmosphere  leads  to  NO  and  OH  species, 
which  in  turn  affect  the  ozone.  For  a discussion  of  NO  and  GH 
production,  see  Section  E.3.  The  next  question  is  whether  the 
NO  and  OH  molecules  produced  by  ionization  at  a given,  relatively 
high  altitude  z^  can  get  down  to  where  the  ozone  is.  See  Fig. 
E-l,  which  shows  that  at  low  and  mid-latitudes  this  is  mainly 
below  40  km,  while  at  high  latitudes  it  is  mainly  below  30  km. 
High  (geomagnetic)  latitudes  are  particularly  significant  be- 
cause the  earth's  magnetic  field  deflects  charged  particles 
away  from  lower  geomagnetic  latitudes,  but  not  from  the  polar 
cap  regions  in  which  the  flux  of  charged  particles  is  approxi- 
mately equal  to  the  flux  in  free  space. 

Figure  E-2  shows  the  polar  cap  regions  that  are  particularly 
affected  by  cosmic  rays,  both  solar  and  galactic.  The  inner 
closed  curve  corresponds  roughly  to  geomagnetic  latitude  A > 
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Source:  Reid,  1974 

FIGURE  E-2.  The  Northern  and  Southern  "Polar  Caps."  SPE  effects  are 
important  inside  curves  , while  GCR  effects  are  most 
important  inside  C . 


1 *3  p 

65°  (total  area  4.8  x 10  J m in  both  hemispheres);  this  is  the 

# 

region  in  which  solar  protons  produce  the  biggest  effect . The 
outer  curve  C , which  corresponds  approximately  to  geomagnetic 
latitude  A > 60°  (total  area  6.9  x 10J'^  ) is  that  over  which 
the  galactic  cosmic  ray  flux  has  its  plateau  value,  and  should 
thus  be  used  for  GCR  effects. 

There  is  a question  of  how  these  curves  should  be  modeled 
in  a two-dimensional  photochemical/dynamic  model  which  uses  geo- 
graphic latitude  as  one  parameter.  Of  course,  there  is  no  very 
simple  correlation  between  geographic  and  geomagnetic  latitude. 

The  following  procedure  is  suggested.  Let  the  polar  caps  be 
divided  into  annular  zones  of  geographic  latitude  having  5-deg 
width,  i.e.,  going  from  60  to  65°,  65  to  70°,  etc.  For  each 

annular  zone,  we  ask,  what  fraction  f.  (j  = i,o)  of  the  area  of 

J 

the  zone  lies  inside  curve  C.  of  Fig.  E-2.  Table  E-l  and  also 

J 

Fig.  E-3  show  the  area  A of  each  zone,  the  f.  (K  = north , south)  , 

v-J 

and  also  the  effective  area,  weighted  by  f.  , that  should  be  used 

•J 

here.  Note  the  asymmetry  between  Northern  and  Southern  Hemis- 
pheres; the  effects  extend  from  40  to  90°,  but  maximize  in  the 
60  to  70°  latitude  range. 

Note  in  Table  E-l  that  E f (N+S)^  = 5.8  x 101-5  , whereas 

1 

I O p 

the  area  of  the  polar  caps  defined  as  A > 65°  is  4.8  x 10  : m , 
while  E fo^w+°^A  = 7.9  x 10^  m^  as  against  6.9  x 10^  m^  for 
A > 60°.  This  difference  arises  because  of  the  zonal  asymmetry 
of  the  actual  polar  caps  shown  in  Fig.  E-2. 

A tentative  discussion  of  the  penetration  of  NC  and  HO^ 
from  40  to  60  km  down  to  where  the  bulk  of  atmospheric  ozone  is 
~n 

Dr.  D.  J.  Hofmann  (private  communication)  has  pointed  out  that 
during  an  SPE  the  geomagnetic  field  is  distorted  so  that  the 
protons  can  penetrate  to  lower  latitudes  than  normally.  This 
is  an  effect  of  several  degrees  [see,  e.g.,  Zmuda  and  Potemra 
(1973)1  which  is  unlikely  to  be  of  primary  importance  in  the 
present  context. 

* # 

A 5-deg  interval  is  used  here  simply  because  it  appears  roughly 
consistent  with  the  overall  data,  that  is,  using  finer  intervals 
is  not  justified. 
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Located  is  given  in  Annex  EA;  this  suggests  that  HO  radicals 
are  destroyed  so  rapidly  by  photochemistry  that  they  cannot  be 
transported  through  significant  vertical  distances  ("significant" 
here  means  more  than  a few  kilometers).  NO  molecules  are  als 
destroyed  very  rapidly  above  about  45  km,  but  at  lower  altitudes 
their  photochemical  lifetime  is  so  long  that  they  can  be  trans- 
ported downward  to  where  they  can  affect  the  bulk  >f  the  z >ne. 

Thus,  we  conclude  that  ionization  produced  abov<  ki  Li 
only  of  local  significance  in  affecting  atmospheric  osorv  , but 
ionization  produced  at  lower  altitudes  will  affect  th-  . u, 
an  overall  basis.  Accordingly: 

• Primary  solar  electromagnetic  radiation  is  not  si  mil L- 
cant  because  it  oroduces  ionization  r rinc ir  allv  -.r  ev-- 
80-km  altitude. 


• 

Solar  protons 

of  energies 

100  MeV  or  greater  are 

found 

in  re  Lati  v<  . y 

large  solar 

proton  events  (formerly 

’ kn<  wn 

ns  PGA  from  t 

heir  effects 

on  ionospheric  radif  p: 

• : Iga- 

tion) . Phi  Large  SPEs  ar«  relative! y rar<  and  juil 
complex  ; :•  oes. , which  ->  •••  i ioousse-l  ir  Se  :ti  E.5. 

• fala  ■ Lc  (as  iistinct  from  solar)  c smi : rays  ( >CRs)  have 

9 

energy  in  the  1 to  10  GeV  range  (1  GeV  = 10'  eV) , and 
show  a modulation  with  solar  activity,  having  a maximum 
at  solar  minimum.  These  processes  are  discussed  in 
Section  E.4. 

E.3  PRODUCTION  OF  NO  AND  OH  BY  IONIZATION  IN  THE  STRATOSPHERE 

When  protons  or  other  ionizing  particles  travel  through  the 
middle  atmosphere,  they  lose  energy  and  produce  ions  at  the  ca- 
nonical rate  of  35  eV/ion  pair.  As  a result  of  charge  transfer 
and  ion-molecule  reactions,  the  ions  actually  observed  are  not 

those  formed  initially,  which  are  N+,  N+,  ot,  and  0+ , but  rather 

^ d + 

those  of  lowest  ionization  potential,  which  are  NO  in  dry  air  and 

-f  + 

and  also  , Hj-;h.  , etc.  (hydrated  protons)  in  the  actual  atmo- 

sphere, which  contains  i t 10  ppmv  of  water  vapor. 

E-13 


The  principal  ions  in  the  lower  mesosphere  are  H-V  and 
H.O^  [cl'.,  e.g.,  Heicklen  (19765  p.  127)  and  Ferguson  (1975) 
for  a review  of  stratospheric  ion  chemistry].  There  are  no  de- 
tailed results  on  stratospheric  ions,  although  observations  from 
a balloon  at  27  to  32  km  indicate  mass  numbers  73+,  91+ , 109+, 

which  correspond  to  H+(H  0)  , for  n = 4,  5,  6 — see  Olson  et  al. 

2 n 

(1977). 

The  NO  and  OH  molecules  are  produced  by  a variety  of  reac- 
tions. Some  currently  accepted  schemes  are  the  following: 

NO • The  source  of  odd  nitrogen  is  impact  dissociation  or 
dissociative  ionization  of  atmospheric  N-,  by  fast  (10  to 
100  eV)  secondary  electrons  produced  from  the  incident 
protons 


/N2 

+ 

2e" 

(a) 

2 + e\ 

— -N 

+ 

N + 2e~ 

(b) 

Xx2N 

+ 

e~ 

(c) 

followed  by  the  reactions 

+ e -*■  2N  (some  atoms  are  excited) 
N + 0^  (a) 

NO+  + 0 (b) 

NO  + e -*•  N + 0 (some  atoms  are  excited) 


N + 02  - NO  + 0 
N + 0^  NO  + 0. , 


(1) 


(2) 

(3) 

(*0 

(5) 

(6) 


which  produce  NO.  There  is  also  the  destruction  reaction 
for  odd  nitrogen 


N + NO  -*■  N,,  + 0 


(7) 
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Different  versions  of  this  scheme  give  the  following  number 
of  N0x  molecules/ion  pair: 

1.4  (factor  2)  Dalgarno  (1967) 

1.2'-  1.3  Frederick  (1976) 

1.5  Crutzen,  Isaksen,  and  Reid  (1975) 


OH.  After  the  initial  ion  formation,  charge  transfer  and 

ion-molecule  reactions  lead  to  the  formation  of  ot,  then 

oj,  and  then  either  to  (H^0+  + OH),  or  else  to  H^0+.  OH, 

which  reacts  with  another  H„0  to  make  H,-0~.  Thus,  one  OH 

2 5 2 

radical  is  produced  in  the  formation  of  the  hydrated  hydro- 
nium  ion.  The  ion  is  neutralized  by  dissociative  recombina- 
tion which  yields  yet  another  OH  radical, 


H30+  + e"  - OH  + H2.  (8) 

This  is  the  scheme  of  Swider  and  Keneshea  (1973),  who  find 
two  OH  radicals/ion  pair;  Frederick  (1976)  uses  a different 
scheme  but,  again,  predicts  two  OH  radicals/ion  pair  in  the 
stratosphere/mesosphere  below  60  km. 

In  the  present  context,  following  M.G.  Heaps  (1978b)  as 
described  by  Fig.  E-4 , we  suggest  that  one  should  assume  that 
1.3  NO  molecules  and  2.0  OH  radicals  are  produced  per  ion  pair 
at  the  altitudes  of  interest  here,  i.e.,  below  45  km. 


E . 4 STRATOSPHERIC  IONIZATION  DUE  TO  GALACTIC  COSMIC  RAYS  (GCRs): 
VARIATION  WITH  SOLAR  CYCLE* 

_ O _ 1 

The  ionization  rate  Q (ion  pairs  cm  Js  ) due  to  GCRs  shows 
considerable  variation  due  to  latitude,  as  well  as  solar  sunspot 
cycle,  as  is  illustrated  in  Fig.  E-5 . 

The  latitude  dependence  can  be  simplified  such  that  for  the 
polar  cap  regions  (corresponding  to  those  geographic  latitudes 
for  which  the  magnetic  latitude  A > 60° — see  also  Fig.  E-2  and 

# 

This  section  has  been  extensively  revised  by  M.G.  Heaps. 
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ALTITUDE, km 


60 


Source:  M.G.  Heaps,  1978  b 

FIGURE  E-4.  NO  and  OH  Production  Per  Ion  Pair. 
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related  discussion)  the  Ion  pair  production  rate  may  be  repre- 
sented as  (Heaps,  1978a) 


Q 


C + D cos 


( YM-year ) 


P 


(9) 


(Heaps,  1978a)  where  C = 0.38,  p 
at  the  altitude  of  interest,  and 
1976  (year  of  OCR  maximum,  i.e., 


= 0.13,  p = pressure  (mbar) 

YM  = 1954,  1965,  and  nominally 
solar  sunspot  minimum) 


year  = current  year  of  interest 

P = period:  9 years  going  down  from  YM 
13  years  going  up  from  YM 


Over  the  polar  cap  regions  the  value  of  Q depends  linearly 
on  pressure  as  shown  in  Eq  (9)  down  to  about  20  km  (^  50  mbar) 
(see  Fig.  1 in  Heaps,  1978a),  but  then  increases  at  a much 
slower  rate  down  to  approximately  15  km  and  actually  begins  to 
decrease  below  15  km. 


Below  20  km  these  approximations  may  be  used  over  the  polar 

caps 


Q = 


Q = 


2.6  + 0.8  cos 


(YM-year ) 


p°-5  (i5_20  km, 


70  + 21  cos 


(YM-year ) 


n 

(J 


p 


-0.2 


(10-15  km, 


50-110  mbar) ( 10 ) 


110-250  mbar ) ( 11 ) 


At  middle  and  low  latitudes  corresponding  to  A < 60°  (or, 
lying  outside  curves  Cq  of  Fig.  E-2)  the  ionization  rate  above 
31  km  may  be  expressed  as  (Heaps,  1978a) 


Q = (A  + B sin  A)p  (above  31  km,  p < 10  mb) 


(12) 


(Heaps,  1978a)  where 
minimum  and  B = 0.51 
sinusoidal  variation 


A = 0.0^6,  B = 0.75  for  a solar  sunspot 
for  a solar  sunspot  maximum.  A quasi- 
with  solar  sunspot  cycle  similar  to  that 
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1 


expressed  in  Eq  (9)  may  be  inferred,  but  there  is  no  extensive 
lower  latitude  data  to  confirm  this. 

At  lower  latitudes  the  ionization  rate  due  to  GCRs  actually 
increases  faster  than  the  ambient  pressure  in  the  31-18  km  range 
such  t-hat 


Q 


(A  + B siri  A)ppn  (18-31  km,  10-70  mbar) 


(13) 


where  p = 10  mbar,  n = 0.8  + 0.6  cosA  and  y = 1 - n.  A and  B 

o 

are  the  same  as  in  Eq  (11). 

In  the  15-18  km  range  the  ionization  rate  may  again  be 
approximated  as  proportional  to  the  pressure  and  Eq  (11)  may  be 
used  again  if  it  is  multiplied  by  the  factor  "*"  ■ 

The  ionization  rates  given  in  Eqs  (9)  to  (13)  can  be  inte- 
grated over  the  volume  of  the  stratosphere  to  yield  a total  ion 
pair  production  for  solar  sunspot  maximum  (i.e.,  cosmic  ray 
minimum),  which  should  serve  as  a baseline  and  integrated  again 
for  solar  sunspot  minumum  to  obtain  the  difference  in  production 
over  the  solar  cycle. 

Over  the  polar  cap  regions  a nominal  tropopause  level  shall 
be  taken  as  10  km,  which  shall  serve  as  a lower  boundary  on  the 
spatial  integration.  Using  the  values  of  Q from  Eqs  (9)  to  (11) 
one  has 


/ TTP 

0 0 0 z 


Q(r  dr  sinadotd<j> )dt  = / / QdVdt 

t C 


// 


(14) 


i 


for  total  yearly  production,  which  yields 

11  QdVdt  = 6.5  x 1032  ion  pairs/yr 


t C 
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for  the  total  Ion  pair  production  in  the  stratosphere  over  each 
polar  cap  for  a solar  sunspot  minimum  (i.e.,  OCR  maximum),  and 


3-4 


x 10 


32 


ion  pairs/, yr  , 


for  a solar  sunspot  maximum. 

The  middle  and  low  latitude  regions  also  contribute  appre- 
ciably to  the  total  ion  pair  production  because  of  the  much 
larger  surface  area  involved,  even  though  the  ionization  rate  is 
reduced.  For  these  lower  latitudes  a nominal  tropopause  level 
of  15  km  is  chosen  to  serve  as  the  lower  bound  on  the  spatial 
integration.  Using  the  values  from  Eqs  (12)  and  (13)  over  the 
appropriate  altitude  range,  one  has 


rl  yr  tt/2  2ir  z ? r r 

I III  Q(r  dr  sinadad'})  )dt  = / f 


QdVdt 


(15) 


0 tt/6  0 


t ( 1-C  ) 
o 


for  the  total  yearly  production  over  middle  and  low  latitudes. 
Her e , Jp  _r,  ^ JV  means  the  integration  goes  over  a hemisphere 
but  excludes  the  polar  cap. 


QdVdt 


U-c0) 


13  x 1032 


ion  pairs/yr, 


the  total  ion  pair  production  in  the  stratosphere  over  a hemis- 
phere (0  to  60°)  for  a solar  sunspot  minimum  (i.e.,  GCR  maximum). 

QdVdt  = 10  x lC3<r'  ion  pairs/yr 

for  the  solar  sunspot  maximum. 
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Thus,  one  can  see  that  the  total  stratospheric  ion  pair 
production  over  both  polar  caps  during  a GCR  minimum  is  on  the 
order  of  7 x 103^  ion  pairs/yr,  with  a total  global  stratospheric 
production  on  the  order  of  27  x lO3^  ion  pairs/yr.  These  values 
represent  the  "baseline"  values.  During  a solar  sunspot  minimum 
(GCR  maximum)  an  additional  6 x 10''2  ion  pairs/yr  are  produced 
over  the  polar  caps,  with  an  additional  total  global  production 
of  12  x 10  32  ion  pairs/yr. 


When  one  adopts  the  odd-nitrogen/ion  pair  production  ratio 
of  1.3,  the  baseline  values  for  odd-nitrogen  production  become 


8.8  x 10 32  NO/yr 


stratospheric  production  over  both 
the  polar  caps  during  GCR  minimum 


8.1  x 1032  NO/yr 


additional  polar  cap  production  dur- 
ing GCR  maximum 


3.5  x 1033  NO/yr 


total  global  stratospheric  production 
during  GCR  minimum 


1.6  x 10 


33 


NO/yr 


additional  global  production  during 
GCR  maximum 


Considering  a total  stratospheric  loading  of  (4  to  10)  x 1033 
NO  molecules,  the  additional  production  over  the  polar  caps 
during  GCR  maxima  should  be  a slowly  varying  perturbation  of  < 2 
percent/yr,  with  a total  global  perturbation  < 4 percent/yr. 


With  an  odd-hydrogen/ion  pair  production  ratio  of  2.0  below 
30  km,  the  baseline  and  additional  values  of  odd-hydrogen  become 
3 7 

1.4  x 10  J OH/yr  stratospheric  production  over  the 

polar  caps  during  GCR  minimum 


1.2  x 103j  OH/yr 


additional  polar  cap  production  dur- 
ing GCR  maximum 


5.4  x 1033  OH/yr 


total  global  stratospheric  production 
during  GCR  minimum 


2.4  x 1033  OH/yr 


additional  global  product  ion  during 
GCP.  maximum 
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The  annual  production  of  OH  due  to  GCRs  if.  still  smaller 

than  the  production  rate  from  the  action  of  sunlight  and 

* 

m H,  . Therefore,  the  global  ■ ffect  n 1 ading  is  3t  ill 
Likely  to  be  a perturbation  of  a few  percent,  although  the  { - 
tential  effect  in  the  winter  polar  stratosphere  may  be  much 
larger . The  H0y  loading  over  the  winter  polar  cap  may  even 
double  due  to  the  enhanced  ionization  during  a GCR  maximum. 

E . 5 STRATOSPHERIC  IONIZATION  DUE  TO  SOLAR  PROTON  EVENTS  (SPEs) 

In  contrast  to  galactic  cosmic  rays,  which  are  a very 
steady  source  of  ionization,  modulated  principally  by  the  11- 
year  solar  cycle,  solar  proton  events  are  highly  variable  impul- 
sive events  resulting  from  some  (but  not  all)  solar  flares.  For 
orientation.  Table  E-2  lists  major  SFEs  during  the  time  frame 
1956  to  1972.  Here  "major’’  refers  to  a Shea-Smart  index  (de- 
fined in  footnote  of  Table  E-2)  greater  than  X20. 

In  fact,  the  Shea-Smart  index  has  significant  limitations 
for  the  present  application  because  different  solar  flares  pro- 
duce significantly  different  numbers  and  energy  distributions 
of  high  energy  protons — see  Annex  E to  this  Appendix,  item  1. 

Given  the  total  spectrum  for  a given  SPE  (see  Annex  B to 
this  Appendix,  item  1),  it  is  possible  to  calculate  the  time- 

integrated  Ionization  due  to  the  event.  Figure  E-6  shows  this 

# # # 

for  a number  of  major  events.  Next,  we  ask  which  SPEs  should 


Note  that  while  the  estimated  H0X  loading  of  the  stratosphere 
is  only  (7  - 16)  x 1032  n Lecules  (see,  e.g.,  append'.?:  A) 
the  lifetime  of  H0X  is  so  short  (less  than  30  min  above  40  kr.-- 
see  Table  EA.-l)  that  the  effect  on  these  annual  production 
rates  is  relatively  small. 

#* 

No  major  events  were  recorded  In  1955,  and  there  were  no 
major  events  during  1973  to  1975. 

«## 

It  also  shows  the  ionization  due  to  relativistic  electron  pre- 
cipitation (REP)  events  (Thorne,  1977,  1978)  which  occur  a; - 
proximately  3 to  4 percent  of  the  time.  Their  contribut i 
evidently  net  very  important. 
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TABLE  E - 2 . SIGNIFICANT  SOLAR  PROTON  EVENTS,  1 956  -1  973 


DATE 

SHEA  & Smapt 
RATING4  ■' 

PEAK  FLUX  RATING 
( > 30  Mev)b,<: 

OLE  (X) 

PCA  ( d B ) 

FLARE 

CAUSING1 

2 3 

Fet* 

56 

X 34 

3 & 

>4550  (c  )( 
19000  (d)f 

13  (c.b.d) 

3 

31 

Aug 

66 

X 21 

- 

3 (c ,d) 

4.9  (c.b.d) 

3 

1 3 

Nov 

56 

X30 

Nl 

- 

5.4  (c.d) 

2 & 

22 

Jun 

57 

X 30 

34 

3 (c  ,d) 

4.9  (c.b.d) 

? 

3 

Jul 

57 

X30 

- 

- 

6-9.2  (c.b.d) 

3 + 

29 

Aug 

57 

X30 

24 

- 

8.2-9  (c.d) 

? 

31 

Aug 

57 

X30 

NL 

- 

4.9-5  (c.d) 

3 

2 

Sep 

57 

X30 

NL 

- 

7.2-9  (c.d) 

24 

21 

Sep 

57 

X 30 

NL 

- 

5.1  (c.d) 

3 

20 

Oct 

57 

X 30 

- 

- 

5-7.8  (c.b.d) 

3 + 

10 

Feb 

58 

X30 

- 

- 

3.2-12  (c.b.d) 

2* 

25 

Mar 

58 

X30 

34 

- 

10-12  (c.d) 

? 

7 

Jul 

58 

X40 

34 

- 

>15-23.7  (c.b.d) 

3 + 

16 

Aug 

58 

X40 

2 & 

- 

1 2 . 1 -> 1 5 (c.b.d) 

3 + 

22 

Aug 

58 

X30 

24 

- 

>10-10.6  (c.b.d) 

3 

26 

Aug 

58 

X40 

34 

- 

16.6  (c.b.d) 

3 

22 

Sep 

58 

X30 

NL 

- 

4-5  (c.d) 

2 

10 

May 

59 

X40 

34 

- 

>15-22  (c.b.d) 

3 + 

10 

Jul 

59 

X40 

34 

- 

>15-20  (c.b.d) 

3 + 

14 

Jul 

59 

X40 

4& 

- 

>15-23.7  (c.b.d) 

3 + 

16 

Jul 

59 

X42 

4& 

>10  (c)t 

1 5 (d)( 

>15-21 .2  (c.b.d) 

3 + 

30 

Mar 

60 

X30 

U 

- 

3-7  (c.b.d) 

2 

29 

Apr 

60 

(1)30 

2 & 

* 

11.2-14  (c.d) 

2 + 

4 

May 

60 

(1)24 

2 & 

1290(c)! 

♦ 175(d)! 

3.4-5  (c.b.d) 

3 

6 

May 

60 

(1)40 

1& 

- 

8 . 7 - > 1 5 (c.b.d) 

3 + 

3 

Sep 

60 

221 

2 & 

4.5  (c.d) 

2.5-4  (c.b.d) 

2 + 

12 

Nov 

60 

444 

" 5 " & 

>135  (c)l 
1100  (d)f 

>14->22  (c.b.d) 

3 + 

1 5 

Nov 

60 

443 

3& 

>88  (c)  ! 
1150  ( d ) f 

>22  (c.b.d) 

3 

20 

Nov 

60 

(3)32 

34 

1?  15!) 

3-8  (c.b.d) 

.2 

E-2  3 


TABLE  E - 2 . (cont) 


DATE 

SHEA  & SMART 
RATING3 

PEAK  FLUX  RATING 
( > 30  Mev)b’2 

GLE  (S) 

PCA  ( dB ) 

FLARE  , 
CAUSING 

12  Ju  1 

61 

(2)40 

26 

- 

17-20  (c,b,d) 

3 

16  Jul 

61 

333 

36 

24  (c) 

8.7-10  (c,b,d) 

3 + 

20  Jul 

61 

(1)22 

28 

7 (c) 

5 (c,b) 

3 

26  Sep 

63 

(1)30 

06 

- 

3-4.6  (c.b.d) 

3 

7 Jul 

66 

121 

U 

2.5  (c) 

2-2.5  (c,b) 

28 

2 Sep 

66 

230 

26 

- 

12-14  (c,b,d) 

38 

28  Jan 

67 

1 33 

16 

) 2 1 (c)l 
(20  ( d ) I 

7 (c.b.d) 

7 

23  May 

67 

330 

1 6 

- 

11  (c.b.d) 

38 

9 Jun 

68 

230 

1 6 

- 

6.5  (c.b.d) 

3B 

29  Sep 

68 

121 

1 6 

(1.4  (c)l 

1 2 (d)  1 

1-7  (c.b) 

2B 

31  Oct 

68 

230 

1 6 

- 

6.5  (c.b.d) 

3B 

1 Nov 

68 

230 

NL 

- 

4-5.9  (c.d) 

2N 

18  Nov 

68 

233 

26 

(14  (c)l 

113  ( d ) f 

12.5  (c.b.d) 

18 

5 Dec 

68 

230 

1 6 

- 

4.7  (c.b.d) 

7 

25  Feb 

69 

123 

1 6 

716  (c)  ( 
115  (d)  ) 

2.1  (c.b) 

2B 

30  Ma  r 

69 

112 

16 

(8.8  (cH 
»5  (d)  ( 

1.3  (c.b) 

7 

11  Apr 

69 

330 

26 

- 

>12  (c.b.d) 

N 

2 Nov 

69 

330 

26 

- 

14.5  (c.b.d) 

N 

No  Shea  and  Smart8  or  Svestka  and  Simonc  ratino  after  December  1969. 


DATE 


SHEA  6 SMART 
SYNTHESIZED 
RATING 


PEAK  FLUX  RATING 
( >30  Mev)b'2 


GLE  (S) 


PCA  ( d B ) 


TABLE  E-2. 


( cont ) 


NOTE  : 

Shea  and  Smart's  ratings  are  listed  below.  Only  events  greater  than  (X20)  are  given 
in  this  table. 


Solar  Proton  Event  Classification  System  (Smart  and  Shea,  1971) 


First  digit 

Second  digit 

Third  digit 

0 i g i t 

>10  MeV  proton 
flu*  (cm^  s ster)’1 

Daylight  polar 
cap  absorption 
at  30  MHz 

Sea  level 
neutron  monitor 
Increase 

-2 

10‘2-00_1 

- 

- 

-1 

o 

o 

o 

- 

- 

0 

10°-<101 

No  increase 

No  Increase 

I 

T 01 - «7  02 

0.5  dB 

<31 

2 

102-003 

1 . 5 - <4 . 6 dB 

3- <1  OX 

3 

103-  O04 

4.6-05  dB 

lO-OOOX 

4 

>104 

> 15  dB 

> 100X 

X • measurements  not  available 
()  « the  digit  Is  uncertain  or  implied 


2.  The  Sauer  and  Stonehocker  peak  flux  rating  (Ref.  b,  Table  1,  p.  5)  used  here  is 

3 4 2 

analogous  to  Shea  and  Smart's  first  digit,  i.e.,  3 for  a flux  10  -10  protons/cm  -sec-ster. 

We  write  3S  instead  of  3 to  indicate  that  Sauer  and  S tonehocker ' s flux  refers  to  proton 
energy  >30  MeV,  as  against  Shea  and  Smart's  index  which  refers  to  energies  >10  MeV. 

NL  - means  "event  not  listed" 

4 dash  (-)  means  "no  enhancement  quoted" 

References : 

a.  Shea  and  Smart  (1977)  - data  from  1955  through  1969 

b.  Sauer  and  Stonehocker  (1977)  - data  from  1956  through  1972 

c.  Svestka  and  Simon  (1975)  - data  from  1955  through  1969 

d.  Pomerantz  and  Ouggal  (1974)  - data  from  1942  through  1973 


be  considered.  Some  details  are  reviewed  in  Annex  B ‘o  m 


Appendix,  where  item  2 presents  some  general  factors,  wr.ile-  i'  :r 
3 explains  one  method  of  calculating  the  total  ionization  due  to 
protons  of  energy  above  100  MeV.  Table  E-3  summarizes  the  re- 
sults and  indicates  those  events  which  ought  to  be  considered  as 
contributing  at  least  several  times  ic'^/  ion  pairs  below  45  k-r 
or  so . 

The  following  comments  should  be  made: 

(a)  The  estimates  for  total  ionization  due  to  Reid,  and 
those  made  here  for  the  events  between  1956  and  I960, 
use  a slightly  different  treatment  of  the  same  data. 
Reid  provides  the  various  altitude  profiles  of  ioniza- 
tion shown  in  Fig.  F.-6 , while  in  item  3 of  Annex  B so 
this  Appendix  we  simply  compute  the  total  energy  car- 
ried by  protons  of  100  MeV  or  more  (which  penetrate 
down  to  perhaps  32-km  altitude)  and  estimate  the  total 
ionization  from  the  standard  factor  of  35  eV/lon  pair 
in  air.  The  agreement  between  the  November  I960 
events  is  excellent  (17  vs  13.7,  for  the  number  of 
ion  pairs  in  units  of  10^);  for  the  February  1956 
event  the  two  estimates  differ  by  a factor  of  2.7, 

but  it  is  clear  that  this  was  a much  smaller  event 
overall . 

(b)  We  suggest  that  the  ionization  profile  for  the  May- 
July  1959  events  should  be  scaled  from  that  for  the 
November  I960  event  by  multiplying  by  a scale  factor 
proportional  to  the  respective  ionization,  which  is 

(4.8  + 1.7  + 5.7  + 6-1)  _ 18.3  _ , 

(8.2  + 5.5)  13.7  ~ J ' 

(c)  The  events  of  February  1956,  May  I960,  and  September 

# 

1966  are  so  small  that  they  can  probably  be  ignored; 

_ 

Using  the  scaling  of  (b)  above,  the  May  I960  event  is  approxi- 
mately one-eighth  times  as  large  as  the  February  1956  event . 
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TABLE  E -3 . SOME  LARGE  SOLAR  PROTON  EVENTS  AND  THEIR  IONIZATION 


^ee  Appendix  D,  Item  3 

^Three  figures  are  not  significant 

^G.C.  Reid,  private  communication,  Dec.  1977 


I have  no  calculations  for  the  November  1969  event, 
which  was  certainly  smaller  than  the  other  events 
listed . 

(d)  In  conclusion,  the  events  of  summer  1959,  November 
I960,  and  August  1972  should  probably  be  taken  into 
account,  although  their  effects  on  NO^ — and  probably 
also  on  ozone-~is  certainly  less  than  that  of  the 
largest  individual  nuclear  test  series  (see,  e.g.. 
Appendix.  B)  . 
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ANNEX  A TO  APPENDIX  E 
DOWNWARD  TRANSPORT  OF  NOx  AND  HOx 
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ANNEX  A TO  APPENDIX  E 


DOWNWARD  TRANSPORT  OF  NOx  AND  HOx 

The  question  discussed  here  is  the  following.  NO  and  HO 
are  produced  by  ionization  throughout  the  stratosphere  and  meso- 
sphere, and  they  will  affect  the  ambient  ozone  at  the  altitude 
of  production.  However,  much  ionization  occurs  at  relatively 
high  altitudes  (above  60  km)  while  most  of  the  o zone  lies  at 
relatively  low  altitudes  (below  25  km),  and  so  whether  the  N0x 
and  H0x  will  affect  the  bulk  of  the  ozone  depends  on  the  rela- 
tive rate  of  transport  as  compared  to  photochemical  destruction. 
These  factors  are  compared  here  in  a very  crude  way.  The  con- 
clusion is  that  below  about  45  km  the  photochemical  lifetime  of 
NO  is  sufficiently  long  that  this  material  can  be  transported 
down  to  where  the  ozone  is,  but  at  higher  altitudes  NO  is  de- 

X 

stroyed  so  very  rapidly  that  transport  is  simply  too  slow.  Over- 
all, the  lifetime  of  H0x  is  so  short  that  vertical  transport  of 
this  species  is  never  significant. 

Throughout  the  stratosphere,  and  up  to  at  least  60  km 
altitude,  NO  is  produced  principally  by  the  reaction 

N20  + 0(1D)  2 NO  . ( EA-1 ) 

There  are  a number  of  photochemical  reactions  which  shuffle  ’’odd 
nitrogen”  = N + NO  + NO,., , in  particular 


I should  like  to  thank  Dr.  J.  Zinn  (LASL)  for  providing  the 
information  on  chemical  production  and  destruction  mechanisms 
and  lifetimes  used  here. 
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NO  + hv  N + 0 


N + 0j  ->  NO  + 0 (EA-2) 

N + 0^  ->  NO  + 0;) 

Net:  N + 0^  + hv  -+  NO  + 2 0 


NO 

+ 

03  ->  NOp 

+ 

°2 

NO., 

+ 

hv  ->  NO 

+ 

0 

NOp 

+ 

0 -*■  NO 

+ 

°2 

(EA-'j ) 


Net : NO  + NO,,  + 20^  + hv  2 NO  + 2 0 + 2 0^ 

and  combining  these 

N + 0^  + hv  -*■  NO  + 20 

NO  + NOp  + 2 03  + hv  -*■  2 NO  + 20  + 2 O0  (EA— ^ ) 

N + NO  -+  N0  + 0 


Net:  2 N + NO,,  + 2 0^  + 2 hv  NO  + 3 0 + 2 Op  + Np 

so  that  we  get 


3 odd  nitrogen  1 odd  nitrogen  (EA-5) 

principally  by  the  reaction 

N + NO  -*•  N + 0 (EA-6) 


HOy  is  made  and  destroyed  principally  by  the  reactions 

0(1D)  + HpO  -*•  2 OH  (EA-7) 

OH  + HOp  - HpO  + Op  (EA-8) 


E-34 


The  characteristic  lifetimes  for  both  N0x  and  H0x  are  given  in 
Table  EA-1 . 

TABLE  EA-1.  PHOTOCHEMICAL  LIFETIMES  OF  N0X  AND  H0X* 


Altitude, 

Lifetime  of  N0x , 

Lifetime  of  HO 

km 

sec 

sec 

"high" 

"shorter" 

- 

52 

3 x 105  (=  3.5  days) 

1600  (30  min) 

40 

7 x 107  (2  yr) 

1000  (20  min) 

★ 

Source:  J.  Zinn 


Transport  is  crudely  represented  by  a one -dimensional 
’’eddy  diffusivity"  profile;  in  fact,  the  magnitude  of  the 
transport  coefficient  in  the  40-  to  60-km  altitude  range  is  not 
well  established,  but  the  general  character  of  transport  is  that 
vertical  transport  is  relatively  fast  in  the  40-  to  60-km  range 
but  much  slower  below  25  km  where  the  ozone  is.  Figure  EA-1 
shows  some  currently  used  eddy  diffusivity  profiles. 

Note  that  K is  small  in  the  altitude  range  from  the  tropo- 
pause  to  about  30  km,  but  large  at  both  lower  and  higher  alti- 
tudes. That  is,  transport  is  slow  between  the  tropopause  and 
30  km,  but  fast  at  both  lower  and  higher  altitudes.  To  make  a 
quantitative  estimate  of  the  overall  rate  of  vertical  transport 
is  beyond  the  scope  of  this  study,  but  a qualitative  estimate 
can  be  made  as  described  below. 

From  dimensional  considerations,  transport  through  a dis- 
tance Az  near  altitude  z takes  a time  At 

At  m (Az)2/aK  (EA-9) 
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where  K = mean  eddy  diffusivity  at  altitude  z,  and  a 'v  2 is  a 

numerical  constant.  Thus,  using  the  Wofsy  1975  K-profile  of 
Pig.  EA-1  and  5 -km  steps  for  Az , one  obtains  the  results  of 
Table  EA-2,  which  show  that  transport  in  the  mesosphere  is  very 
fast  ('v  1 week),  slower  in  the  upper  stratosphere  (months), 
very  low  in  the  lower  stratosphere  (years),  and  fast  in  the 
troposphere  ('v  1 week). 


FIGURE  EA-1.  Three  Recent  Eddy  Diffusivity 

Profiles  (From  Bauer  and  Gardner, 


TABLE  EA-2.  REPRESENTATIVE  DIFFUSION  TIME  STEPS 
THROUGH  5- KM  HEIGHT  INTERVALS  [FROM  EQ.  ( EA-9 ) ] 


Altitude 

2 

z , km 

K,  cm  /sec 

t,  days 

60 

2.5  x 105 

6 

55 

2 x 105 

7 

50 

1 

X 

,05 

14 

stratopause 

45 

6 

X 

104 

24 

40 

3 

X 

1 o4 

48 

35 

2 

X 

io4 

70 

30 

1 

X 

.o4 

140 

25 

7 

X 

,o3 

210 

20 

5 

X 

io3 

290 

1 5 

3 

X 

io3 

20 

(through  1 km) 

14 

3 

X 

io4 

30 

(through  4 km) 

10 

1 

X 

io5 

14 

tropopause 

5 

0 

1 

X 

io5 

14 

We  thus  see  that  the  photochemical  lifetime  of  HO  is  so 

X 

short  that  vertical  transport  of  this  species  grouping  is  cer- 
tainly not  significant  in  the  present  context.  For  NOx , things 
are  more  complex;  evidently,  above  about  45  km  the  photochemical 
lifetime  is  so  short  that  downward  transport  is  not  significant, 
but  at  lower  altitudes  vertical  transport  has  to  be  taken  into 
account . 


E-3  J 


r 


ANNEX  B TO  APPENDIX  E 


SELECTED  TOPICS  RELATED  TO  THE 
STRATOSPHERIC  IONIZATION  DUE  TO  SPEs 


ANNEX  B TO  APPENDIX  E 

SELECTED  TOPICS  RELATED  TO  THE 
STRATOSPHERIC  IONIZATION  DUE  TO  SPEs 

1.  THE  INTEGRAL  ENERGY/RIG  ID ITY  SPECTRUM  FOR  DIFFERENT  EVENTS 

Table  EB-1  shows  the  energy  spectra  of  a number  of  large 
SPEs  over  the  1956  to  1961  time  frame.  Note  the  large  differ- 
ence between  the  ratios  of  integrated  intensities  in  excess  of 
30  MeV  vs  100  MeV  for  the  different  events  (last  two  columns). 
Note  also  with  reference  to  Fig.  E-l  that  only  particles  with 
energies  in  excess  of  180  MeV  can  penetrate  at  vertical  inci- 
dence down  to  25  km,  which  is  near  the  upper  limit  of  the  ozone 
layer  at  high  latitudes.  Also,  from  Tables  E-2  and  E-3  we  see 
that  relatively  few  of  the  events  contain  significant  numbers 
of  1 to  10  GeV  protons  which  give  rise  to  a "Ground  Level  En- 
hancement" (GLE). 

It  is  customary  to  characterize  high-energy  particles  in 
terms  of  their  magnetic  rigidity  P,  which  is  a measure  of  their 
ability  to  penetrate  into  the  atmosphere  through  the  geomagnetic 
field  which  acts  as  a magnetic  analyzer,  removing  charged  part- 
icles of  low  rigidity  from  the  flux  incident  on  the  top  of  the 
atmosphere.  Magnetic  rigidity  is  defined  as 

P = pc/q  (EB-1) 

where  P = momentum,  q = charge  of  the  particle,  and  c = speed 
of  light.  P is  customarily  measured  in  volts  (actually  in  GV), 
corresponding  to  the  fact  that  pc  is  measured  in  eV,  and  q in 
units  of  the  electron  charge  e ; it  is  related  to  the  kinetic 
energy  K of  the  particle  through  the  relativistic  relation 
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TABLE  EB-1.  SOME  MAJOR  SOLAR  COSMIC  RAY 
OBSERVATIONS  DURING  1956-1961 
(After  McDonald,  1963) 


Solar  Harr  Data 

Cuamic  Kay  Data 

Maximum 

Onset  or  Kise 

Peak  Flux 

Integrated 

Helio- 

Optical 

Ru  tiio 

Time 

Decay  Time 

• no. 

rn3~'J 

Intensity 

• b) 

' h ) 

sec 

->> 

1 n<>. 

cm”2 ) 

/raphie 

Position 

por 

mum 

> JO 

> 100 

> W 

> 100 

> 10 

> 100 

> 

10 

> IUO 

lljlF 

tance 

• deg » 

• [ T ' 

il'Tl 

MeV 

MeV 

MeV 

MfV 

MeV 

MeV 

Me 

V 

MeV 

1956 

F-b  23 

3^ 

.N22  W~4 

0340 

20.000 

03-11 

6-8 

34 

10 

16 

6.200 

5.000 

6 5* 

10* 

3.2  X 

10* 

1957 

lan  JO 

3-t- 

"25  W.30 

1120 

2.000- 

100- 

3 X 

10* 

1G7 

i 

3.000 

200 

\ldr  23 

3+ 

SO  £78 

1005 

10.000 

1004 

10) 

1 20) 

1.200- 

- 50 

2 < 

10* 

5 X 

10* 

1.500 

Jul  7 

3-f- 

N25  W08 

0115 

2.000 

0112  : 

32 

16-20 

1.5O0- 

2.000 

80 

i X 

10* 

7 x 

10* 

1958 

Vug  16 

34* 

S 14  WS0 

0440 

9.200 

0*42 

10 

18 

200 

2 X 

10: 

\uk22 

3 

N 18  W10 

1 148 

10-12 

34 

20 

8-12 

500 

20 

5 x 

107 

1 X 

10* 

Vug  26 

3 

N20  W54 

0027 

6.  *00 

0025 

^ 9 

' 12) 

1.100 

5.3  X 

107 

Mar  10 

34- 

m 

10.000 

18-22 

12  18 

22 

10-14 

6.000- 

8.000 

1.000 

7 X 

10* 

7.5  X 

107 

10S9 

Jul  10 

3+ 

N22  E70 

0222 

15.000 

0224 

30-40 

18-20 

40 

20 

4.000 

1.200 

8.8  x 

10* 

1.0  < 

10* 

Jul  It 

3-f 

N 16  to: 

0349 

6.300 

0332 

1620 

1248 

18 

9-12 

10.000- 

1.200 

1.1  X 

10* 

6.3  x 

107 

12.000 

juj  J6 

B 

B 

i 

12-14 

4-5 

.10 

18 

16.000- 

18.000 

1.500 

8.1  X 

10* 

1.3  x 

10* 

Vpr  1 

3 

N 12  W10 

0859 

2-3 

< l 

12 

4-6 

50 

6 

2.7  X 

10* 

1.5  x 

10s 

Vpr  5 

2 + 

NIO  W61 

02*5 

8.000 

0203 

12 

40 

2 v 

:o« 

Apr  28 

3 

S05  E3-4 

0130 

8-10 

34 

18 

300 

20 

2.5  X 

107 

7 X 

10> 

24 

8 

Mar  l 

3-*- 

NU  W90 

1020 

2-3 

< 1 

8 

4 

200 

40 

7 X 

10* 

7 X 

10* 

I960 

Sep  3 

3 

N 18  E88 

0110 

12.000 

0108 

1246 

7-9 

32 

26 

240 

o0 

4 < 

107 

7 X 

10* 

Nov  12 

34- 

N27  W02 

1329 

10.000 

1320 

12-16 

8-10 

18-24 

14-18 

12.000 

2.500 

1.4  x 

10* 

3 5 X 

10* 

Nov  15 

J-i- 

N30  W32 

0221 

11.000 

0227  j 

10-16 

3-5 

16-20 

8-12 

6.000 

2.400 

m 

10* 

1.2  x 

10* 

\..V  20 

3 

N28  W 1 1 3 

J020 

1 

1 

34 

- 1 

10-16 

4-6 

1.000 

4.000 

EE 

Q 

6 < 

10* 

ini  n 

3 

>06  E32 

1700 

2.500 

8- 10 

4 

22  26 

18 

20 

3 

2 X 

10* 

3 X 

10* 

Jul  12 

**07  F.22 

1030 

7.500 

842 

9 

16-20 

12 

120 

15 

1.0  < 

10* 

1.6  < 

10* 

*1 

Jul  18 

W- 

**06  W60 

1010 

5.000 

6-10 

2 3 

24 

12 

2.500 

600 

2.1  X 

10* 

48  X 

107 

/ if  JO 

n>7  $90 

. 1600 » 

2.500 

1.6 

1.5 

6-8 

1 

100 

70 

9 X 

10* 

12  x 

10* 

,p  n 

mi  no 

2223 

•S  -3  . 

10* 

9 7 X 

10* 

« it  J to  *0  *M  pH. 


. tlley  -1965,  pp.  17-20. 
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(EB-2) 


(K  + N c2)2  - p2c2  + M V 


where  Mq  = rest  mass  of  particle  (0.9382  GeV  for  protons,  0.5110 
MeV  for  electrons).  Thus,  the  relation  of  rigidity  and  kinetic 
energy  for  protons  is  given  in  Table  EB-2. 

TABLE  EB-2.  RELATION  BETWEEN  RIGIDITY  AND 
KINETIC  ENERGY  FOR  PROTONS* 


Kinetic  Energy, 
GeV 


p2/2M 


o °) 
0.0053 


0 

1 

1 

2 

3 

10 


0.21 
0.125 
.433 
.000 
.27 
. 205 
.0 
.0 


20.0 

pc  (->-“) 

See  Eq  (EB-2) 


Rigidity, 

GV 

pc  (-*  0) 
0.10 
0.20 
0.50 
1 .00 
1 .70 
2.0 
3.0 
3.825 

10.9 

20.9 
pc  (->■  °°) 


The  integral  energy  spectrum  of  solar  protons  is  generally 
described  by  an  exponential  expression  in  terms  of  the  rigidity 
of  the  form 


J(>  P)  = JQ  exp ( ■ 


P/V> 


(EB-3) 


where  J(>  P)  is  the  flux  of  particles  with  rigidity  greater  than 

2 


P,  expressed  in  protons/cm  -sec-ster. 
for  a number  of  large  SPEs  is  given  in  Table  EB-3. 
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A list  of  J and  P values 
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TABLE  EB-3.  INTEGRAL  EXPONENTIAL  SPECTRA 
FOR  SOLAR  PROTONS-  -EQ  (D-3) 


Event  Date 

Optical  Flare  Max. , LIT 

Time  of  Measurement,  UT 

2 

JQ  Protons/cm  -sec-ster 

P , MV 

0 

23  Feb.  56 

0342 

1930 

23  Feb. 

140 

300 

10  May  59 

2118 

0500 

12  May 

30x1 03 

65 

0520-0915 

12  May 

18xl03 

70 

1335-1750 

12  vay 

2.8x1 03 

66 

10  July  59 

0222 

0755-1625 

11  July 

7.8xl03 

90 

1225-2345 

12  July 

2.1x1 03 

90 

14  July  59 

0349 

1330-1555 

15  July 

35x1 03 

68 

16  July  59 

2145 

1200 

17  July 

3. 2x1 03 

125 

1555-1935 

1 . 1 xl 03 

113 

1400 

18  July 

800 

90 

1350-1640 

2.8xlOJ 

70 

1445 

20  July 

3C0 

66 

4 May  60 

1020 

1630 

4 May 

6.7 

280 

1700-0200 

5 May 

5.8 

235 

12  Nov.  60 

1329 

1930 

12  Nov. 

1100 

280 

2000 

1250 

240 

0200 

13  Nov. 

3400 

185 

0800 

3400 

155 

1305 

3000 

120 

1830-2000 

2800 

95 

2230 

1000 

105 

0500 

14  Nov. 

1000 

95 

15  Nov.  60 

0221 

0500 

15  Nov. 

250 

375 

1130 

320 

240 

1030-1230 

285 

175 

2130 

3300 

120 

0930 

16  Nov. 

1000 

100 

1430-1730 

300 

100 

IF  00 

2 Aug 

6 

178 

2' 00 

2 Aug 

11 

130 

4 Aug  72 

0640a 

0600 

4 Aug 

5.4xl04 

32 

0900 

4 Aug 

7 . 8x1 03 

96 

2100 

4 Aug 

3.63x1 05 

84 

0600 

5 Aug 

2.6x1 0^ 

70 

0600 

6 Aug 

7.5xl03 

57 

7 Aug  72 

1534b 

0600 

7 Aug 

3 .9x1 03 

49 

0600 

8 Aug 

9 . 1 xlO3 

65 

0600 

9 Aug 

1 ,58xl04 

40 

0000 

10  Aug 

220 

63 

Source:  Results  through  November  1960  from  Freier  and  Webber,  1963. 

*G.C.  Reid,  private  communication.  These  data  for  J and  P are  computed  and  analyzed  from  flux  data  at 
10,  30,  60  MeV.  Thus,  the  higher  energies  are  not  reliable.  Note  that  the  Aug  72  flares  had  a very 
large  flux  of  relatively  low  energy  protons. 

bMcMath  region  11976  produced  a number  of  large  solar  flares  between  1 Aug  and  11  Aug,  in  particular  8 
giving  a large  optical  (Hu)  signal,  6 with  large  10-cm  radio  noise,  and  10  with  large  1-8A  x-ray  flux-- 
see  Hakura  (1972).  We  list  the  Hu  maxima  of  the  two  class  3B  flares. 
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2.  IONIZATION  DUE  TO  SFEs:  GENERAL  COMMENTS 


liven  th  Integral  speetr  J(>  P)  of  Eq  (EB-3)  and  t he 
duration  of  the  event,  it  is  possible  to  compute  the  ionization 
profiles  shown  in  Fig.  E-7  for  several  events--see  also  Figs. 

EB-1  and  EB-2,  which  compare  ionization  rates  due  to  PPEs , GCRs . 
and  a number  of  other  sources.  The  ionization  rate  of  an  CPE 
(or  other  transient  event)  has  to  be  integrated  over  its  dura- 
tion, which  may  extend  for  several  days.  This  integration  is 
evidently  not  a trivial  task  as  the  spectrum  J(>  P)  changes  sig- 
nificantly— see  Table  EB-3  where  we  show'  how  J and  P decrease 

o o 

with  time  for  numerous  events  e.g.,  those  of  10  May  1959,  16 
July  1959,  12  November  I960. 

It  is  clear  that  the  ionization  rates  of  Figs.  EE-1  and 
EB-2,  as  well  as  the  total  ionization  of  Fig.  E-7,  are  subject 
to  significant  errors;  thus,  in  Fig.  EB-2  we  see  that  the  ioni- 
zation rate  profile  due  to  the  November  1969  event  calculated  by 
two  different  workers  differ  by  two  orders  of  magnitude  over  a 
significant  altitude  range.  I do  not  think  that  it  is  possible 
to  estimate  the  error  in  these  estimates,  but  a factor  2 to  3 up 
or  down  in  the  total  ionization  profile  such  as  Fig.  E-7  seems 
plausible . 

Overall,  it  is  not  clear  that  there  Is  any  definite  corre- 
lation between  the  Shea-Smart  index  and  the  total  ionization, 
say  below  45  km  (see  Figs.  E-7  and  EB-1).  Nevertheless,  it  does 

appear  that  only  a very  few  of  the  largest  SPEs  produce  more  than 
32 

10  ion  pairs  in  the  middle  stratosphere,  and  thus  we  shall  re- 
strict consideration  to  the  events  listed  in  Table  E-3 • 

3.  IONIZATION  DUE  TO  SPEs:  SOME  CALCULATIONS 

H.W.  Dodson  has  computed  the  solar  proton  flux  at  various 
times  following  17  solar  flares,  using  the  data  compiled  by 
Freier  and  Webber  (1963).  In  Table  EB-4  we  show  the  flux  for 
protons  of  energy  greater  than  20  MeV,  and  in  Table  D-5  we  give 
her  averaged  conversion  factors  for  energies  > 50  and  > 100  MeV. 

E-45 


t 


Source:  Zmuda  and  Potemra  (1972). 

FIGURE  EB-1.  Ion  production  rates  due  to  different  SPEs  and 

also  due  to  some  other  D-Region  sources  which  are 
defined  in  Zmuda  and  Potemra  (1972).  Note  that 
there  is  no  obvious  correlation  between  the  Shea- 
Smart  index  and  the  ionization  profile  for  these 
relatively  small  events. 
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FIGURE  EB-2.  Ion  production  rates  due  to  SPEs,  and  also  those 
due  to  GCRs  and  REPs.  Two  different  estimates  are  shown  for 
each  of  the  3 SPEs  listed  to  indicate  the  difference  in  ion 
production  rate  estimates  made  by  different  authors,  which 
presumably  gives  a measure  of  the  uncertainty  in  our  quanti- 
tative understanding  of  the  effect.  (Note  that  the  data  of 
Heaps  as  listed  do  not  extend  below  30  km).  The  sources 
used  are  indicated  on  the  figure. 
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From  Fig.  E-l  we  see  that  for  vertical  incidence,  protons 
of  energy  20  MeV  penetrate  down  to  54-km  altitude,  of  energy  50 
MeV  down  to  42  km,  and  of  energy  100  MeV  down  to  32 -km  altitude. 
Thus,  we  arbitrarily  ask  for  the  flux  above  100  MeV  for  the 
events  listed  in  Table  E-2 ; use  Tables  EB-4  and  EB-5  to  obtain 
the  flux  above  100  MeV  as  a function  of  time;  integrate  this 
flux  over  the  time  duration  of  the  event;  multiply  by  the  area 

of  the  polar  caps,  here  defined  from  curve  C.  of  Fig.  E-2  as 

17  2 ^ 

4.8  x 10  cm  ; and  obtain  the  total  number  of  ion  pairs  pro- 

O 

duced  by  each  event  by  multiplying  by  10/35  (i.e.,  average 
energy/proton  = 100  MeV,  average  energy/ion  pair  = 35  eV)  to 
obtain  the  results  shown  in  Table  E-3 . 

This  is  clearly  a very  crude  procedure,  and  one  can  argue 
that  in  view  of  the  transport  considerations  of  Annex  EA,  all 
protons  that  get  down  to  35  to  38  km  (i.e.,  of  energies  to  6 5 
to  85  MeV)  and  which  thus  deposit  most  of  their  energy  below  45 
km,  should  be  considered.  The  present  choice  makes  some  correc- 
tion for  oblique  incidence  and  has  the  advantage  that  for  the 
November  i960  events  it  is  in  very  good  agreement  with  Reid's 
calculation,  and  thus  lends  itself  to  the  scaling  procedure  sug- 
gested in  context  of  Table  E-3. 
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TABLE  EB-4 . SOLAR  PROTONS  (cm  -sec)  OF  ENERGY  GREATER  THAN 
20  MeV  AT  VARIOUS  TIMES  FOLLOWING  17  SOLAR  FLARES 


E-A9 


Source:  Table  reportedly  compiled  by  H.W.  Dodson;  provided  by  E.W.  Hones.  Tor  basic  data,  see 
Svestka  and  Simon  (197S). 


TABLE  EB-5.  CHARACTERISTIC  RIGIDITY,  P (t),  AND  RATIOS  OF  PROTON  FLUXES 
> 3 MeV,  > 50  MeV,  AND  > 100  MeV,  TO  THAT  > 20  MeV  FOR 
VARIOUS  TIMES  IN  AN  "AVERAGE"  SOLAR  PROTON  EVENT 


Time  After 
Beginning  of  Flare, 
hr 

V t 

Q(>  3 MeV) 

Q(>  20  MeV) 

Q(>  50  MeV) 

Q(>  20  MeV) 

Q(>  100  MeV) 

Q(>  20  MeV) 

3 

250 

1.9 

0.60 

0.37 

12 

160 

2.2 

0.43 

0.21 

24 

130 

2.9 

0.36 

0.13 

48 

100 

3.5 

0.29 

0.07 

72 

85 

4.0 

0.25 

0.06 

96 

80 

4.9 

0.22 

0.04 

120 

70 

5.7 

0.19 

0.032 

144 

60 

9.0 

0.13 

0.016 

Source:  Table  reportedly  compiled  by  H.W.  Dodson;  provided  by  E.W.  Hones. 
For  basic  data,  see  Svestka  and  Simon  (1975). 
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CHANGES  IN  INDUSTRIAL  AND  AGRICULTURAL 
NITROGEN  FIXATION,  1955  TO  1975 

F.l  INTRODUCTION 

It  was  suggested  by  Crutzen  (197*0  and  McElroy  (197^0  that 
the  increased  use  of  commercial  nitrogen  fertilizer,  and  also 
nitrogen  fixation  resulting  from  fossil  fuel  combustion,  may  lead 
to  a reduction  in  stratospheric  ozone . As  the  chemically  fixed 
nitrogen  (this  means  nitrogen  present  in  molecules  such  as  ammo- 
nia, amines,  amino  acids,  etc.,  or  nitrates,  as  distinct  from 
the  very  stable  and  inert  molecule  N^)  is  denitrified  (i.e,  , as 
the  compound  breaks  up  to  produce  various  end  products,  includ- 
ing in  particular  ) , some  fraction  of  the  nitrogen  appears  as 
nitrous  oxide  Nf.O,  rather  than  as  ^ . The  fraction  is  normally 
in  the  range  5 to  20  percent,  with  6 percent  as  a current  best 
estimate — see  CAST  (1976). 

The  global  cycle  of  atmospheric  fixed  nitrogen  is  not  well 
known.  Table  F-l  lists  various  current  estimates  of  sources  of 
fixed  nitrogen,  but  the  sinks  of  nitrous  oxide,  which  are  by  far 
the  largest  atmospheric  components  of  fixed  nitrogen,  are  not 

known.  There  may  be  a sink  on  the  ground  [see  Brice  et  al . 

« 

(1977)],  but  there  is  no  presently  known  sink  of  N^O  in  the  free 
troposphere . 


Brice  et  al . (1977)  find  a diurnal  variation  in  N2O  concentra- 
tion at  the  ground  in  England,  and  from  this  infer  a sink  of 
N2O.  However,  J.M.  Bremner  and  A.M.  Blackmer  (see  Denitrifica- 
tion Seminar,  1977)  do  not  observe  this  diurnal  variation  near 
Ames,  Iowa.  See  also  Freney,  Denmead  and  Simpson  (1978)  and 
Cicerone  et  al.  (1978). 
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In  the  stratosphere,  NO  reacts  with  excited  oxygen,  ( ‘ 0) 
and  is  subject  to  photodissociation. 


N20  + 0 ( ) 


.2  NO  ('v  50%) 


N2  + 02  (v,  50%) 


(la) 

(lb) 


N20  + hv  (X  < 230  nn)  N,  + 0<‘LD) 


(2) 


Reaction  (la)  is  the  principal  source  of  stratospheric  NO. 

Roughly  80  percent  of  N2>0  reaching  the  stratospheric  is  destroyed 
by  reaction  (2) — cf.  Schmeltekopf  et  al . (1977). 

It  should  be  noted  that  NO  formed  in  this  way,  by  reaction 
(la),  is  likely  to  be  formed  at  relatively  high  altitudes  (>,  3'- 
km)  because  0(1D)  is  produced  at  high  altitudes  in  the  strato- 
sphere and  above  by  photodissociation  of  ozone  (and  of  0o)  and 
is  quenched  somewhat  by  collision.  Because  of  these  relatively 
high  altitudes  of  formation,  the  NO  produced  in  this  way  will 


destroy  ozone  by  the  "well-known 

catalytic  cycle: 

NO  + 03  -» 

no2  + o2 

(3a) 

no2  + 0 ■* 

C\J 

O 

+ 

O 

•"7 

(3b) 

net  0 + 0^  -*• 

2 02 

(3c) 

There  are,  however,  important  interactions  with  the  H0x  and  C1X 
cycles,  so  that  the  net  effect  of  N20  injection  on  total  ozone 
is  uncertain — see  Section  F.3. 


The  reaction  N2O  + 0(  P)  ■>  N2  + O2  has  a temperature  dependence 
e-l  000/1 } i.e.,  is  only  significant  at  high  temperatures  such 
as  in  combustion,  but  not  in  the  atmosphere. 
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After  nitrogen  is  fixed,  then.'  is  a delay  time  oefre  the 
N2°  produced  in  denitrification  actually  affects  stratospheri 
ozone.  This  time  delay  is  the  sum  of  three  terms,  namely,  a time 
constant  t^  for  denitrification,  a time  constant  t for  transport 
and  mixing  in  the  troposphere,  and  t . , a time  constant  for  strato- 
spheric photochemistry  as  a function  of  N.^0  flux  through  the 
tropopause . 

The  denitrification  time  constant  depends  on  hcv,  the 
material  is  used,  and  may  range  from  a few  months  to  many  years 
(see  Denitrif icat ion  Seminar,  1977).  Evidently,  must  be  re- 
lated to  the  time  for  recycling  of  C0n  in  plants,  and  here  a 
representative  estimate  (Keeling,  1973>  p.  307)  is  that  two- 
thirds  of  agricultural  fixed  carbon  is  recycled  in  approximately 
one  year,  giving  an  effective  recycling  time  for  carbon  of  about 
1 . 5 years . 

Estimates  for  the  tropospheric  "turnover"  time  t,  range  from 
5 to  200  years  (cf.  Johnston,  1977),  but  a current  "best  estimate" 
is  8 to  12  years  (see,  e.g.,  Hahn  and  Junge"  1977,  and  Goldan  et 
al . , 1978).  Thus,  even  if  the  time  constant  t.,  ror  stratospheric 
photochemistry  is  short,  there  is  still  a significant  delay  be- 
tween the  time  at  which  fixed  nitrogen  is  produced  or  used  at  the 
ground,  and  the  time  when  the  N-,0  resulting  from  denitrification 
can  actually  affect  the  stratosphere. 

Anthropogenic  nitrogen  fixation  comes  from  three  principal 
sources,  namely: 

• Industrial  nitrogen  fixation,  mainly  as  ammonia  by  Haber 
synthesis  and  also  as  ammonium  salts  formed  as  a by- 
product in  coke  production. 

• Combustion,  which  leads  to  the  production  of  oxides  of 
nitrogen  (NO,  also  NO^) — essentially  NO  is  formed  in 
equilibrium  when  air  is  heated  above  2000°  K and  then 
"freezes  in,"  i.e.,  is  not  decomposed,  as  the  air  cools. 

• The  large-scale  agricultural  production  of  plants  having 
nitrogen-fixing  bacteria  in  their  root  systems, 
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principally  legumes  (cf.  Evans  and  Barber,  1977)  results 
in  a significant  amount  of  biological  nitrogen  fixation. 


In  each  of  these  instances  the  nitrogen  compound  is  either 
denitrified  directly  or  else  is  taken  up  in  plants  (as  proteins, 
amino  acids,  etc.)  or  in  industrial  products  (such  as  nylon  or 
other  plastics);  oxides  of  nitrogen  formed  by  combustion  are 
transformed  into  nitric  acid,  which  is  washed  out  of  the  atmo- 
sphere by  rainfall  and  is,  in  turn,  taken  up  in  the  biosphere. 


Here  we  review  time-dependent  changes  in  anthropogenic 
nitrogen  fixation  during  the  1955  to  1975  time  frame  (in  Section 
P.2)  in  order  to  provide  a basis  for  estimating  the  effect  of 
the  resulting  N^O  on  global  ozone.  However,  the  brief  discussion 
in  Section  F.3  indicates  that  only  a very  small  and  somewhat  un- 
certain change  in  total  ozone  could  have  occurred  over  the  1955 
to  1975  time  period. 


F . 2 CHANGES  IN  ANTHROPOGENIC  NITROGEN  FIXATION,  1955  to  1975 
F.2.1  Industrial  Nitrogen  Fixation 

Worldwide  nitrogen  production  is  shown  in  Table  P-2,  while 
Table  P-3  shows  nitrogen  fertilizer  production  and  consumption 
with  a coarse  indication  of  the  latitude  distribution  of  both 
production  and  consumption.  Overall  fixed  nitrogen  production 
has  been  rising  at  a very  rapid  rate — 12  percent  per  annum  from 
1955  to  1970,  or  10  percent  per  annum  from  1955  to  1975. 

The  following  points  should  be  noted: 

• In  view  of  the  long  turnover  time  t2  of  NgO  in  the  tropo- 
sphere, as  compared  to  the  time  it  takes  for  material  to 
be  distributed  uniformly  throughout  the  atmosphere 
(approximately  1 to  2 years),  the  latitude  variation  is 
unlikely  to  be  significant. 

• At  this  point,  I am  unable  to  recommend  a unique  value 

of  t^,  the  delay  between  the  production  of  fixed  nitrogen 
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TABLE  F-2.  WORLD  ANNUAL  INDUSTRIAL  PRODUCTION 
OF  FIXED  NITROGEN,  1955  TO  1975* 


Annual  Fixed  Nitrogen  Production 
Year  Ending  30  June  (lo"  kg  nitrogen) 


1950-1954  (avg) 

5, 

.40 

1955 

7 , 

.30 

1956 

8, 

.07 

1957 

8 

.65 

1958 

9, 

.49 

1 959 

1 1 

.80 

1960 

12, 

.87 

1961 

13 

.97 

1962 

14, 

.81 

1963 

1 7 

.11 

1964 

19 

.39 

1965 

21 

.86 

1966 

24 

.49 

1967 

28 

.69 

1968 

32. 

.14 

1969 

35, 

.89 

1970 

38, 

.78 

1971 

41  , 

.15 

1972 

43, 

.00 

1973 

46  , 

.72 

1974 

48, 

.44 

1975 

49 

.53 

★ 

This  is  essentially  all  produced  as  ammonia,  either  by  Haber 
(or  other)  synthesis,  or  as  a by-product  in  coking  plants. 
However,  South  American  guano  production  is  also  included  in 
the  statistics.  It  was  approximately  2 percent  of  global 
production  in  1959  and  tends  to  decrease  both  absolutely  and 
relatively  with  time.  Agricultural  nitrogen  fixation  (e.g., 
in  soybeans  or  alfalfa)  is  not  included,  nor  is  nitrogen  fix 
ation  in  combustion. 

Source:  U.S.  Bureau  of  Mines,  Minerals  Yearbook--rvari  ous 
years,  data  kindly  provided  by  Dr.  R.J.  Foster 


TABLE  F-3.  GLOBAL  PRODUCTION  AND  CONSUMPTION 
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and  the  eventual  appearance  of  N20.  Putting  t^  = 0 gives 
an  instantaneous  response,  which  is  certainly  not  correct, 
but  provides  an  upper  bound  to  the  ozone  changes. 

• It  is  clear  that  most  of  the  world's  supply  of  anthropogeni- 
cally  fixed  nitrogen  is  still  both  produced  and  consumed  in 
the  latitude  range  of  30  to  60°  N,  which  represents  the  de- 
veloped countries  (see  Table  F-3). 

• Note  the  discrepancy  between  global  production  figures  in 
Table  F-2  (which  lists  total  fixed  nitrogen,  estimated  by 
U.S.  Bureau  of  Mines),  and  in  Table  F-3  (which  lists  only 
fixed  nitrogen  used  in  fertilizer,  estimated  by  U.N. 
Agencies).  In  the  U.S. A,  (and  to  some  extent  in  other 
developed  countries),  25  percent  of  fixed  nitrogen  is 
used  for  applications  other  than  as  fertilizer  (explo- 
sives, resins,  fibers  plastics,  and  animal  feeds),  so 
that  the  figures  in  Table  F-2  (which  are  larger  than 
those  of  Table  F-3)  should  be  used. 

• Note  the  discrepancy  between  production  and  consumption 
figures  of  Table  F-3;  the  difference  is  due  to  stockpil- 
ing, possibly  to  diversion  to  nonfertilizer  use,  and 
possibly  also  due  to  errors  in  the  figures. 

F . 2 . 2 The  Combustion  Source  of  Fixed  Nitrogen 

The  four  estimates  for  this  quantity  listed  in  Table  F-l 

# 

range  from  20  to  40  Mt(N)/yr,  with  an  arithmetic  mean  of  29 
Mt(N)/yr.  U.S.  combustion  emissions  in  1975  are  2.2  x 10^  kg 
N02/yr,  or  6.7  Mt(N)/yr  (data  kindly  provided  by  C.O.  Mann,  U.S. 
EPA).  Assuming  that  the  U.S.  represents  25  percent  of  worldwide 
use  of  fossil  fuel,  and  presumably  25  percent  of  N0x  emissions, 
this  corresponds  to  a global  combustion  source  of  27  Mt(N)/yr,  in 
excellent  agreement  with  the  above  values.  We  shall  use  the  29 
Mt(N)/yr  value  for  1975  as  a recommended  estimate. 


1 Mt(N)  = 10^  metric  tons  of  nitrogen,  or  10°  kg 
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The  global  production  of  fossil  fuel  has  increased  by  a 
factor  2.7  between  I960  and  1971,  corresponding  to  an  annual 
growth  rate  of  5 percent  (data  from  Perry  and  Landsberg,  1977); 
data  from  SCEP,  1970,  p.  303,  between  1950  and  1967  correspond  I 

to  an  annual  growth  rate  of  4.5  percent.  Assuming  an  annual 
growth  rate  of  emissions  of  5 percent,  since  these,  presumably, 
track  the  worldwide  use  of  fossil  fuel,  gives  an  estimated  com- 
bustion source  of  10.9  Mt(N)/yr  in  1955. 

£ 

F . 2 . 3 Agricultural  Source  of  Fixed  Nitrogen 

During  the  1955  to  1975  time  frame,  the  total  U.S.  acreage 
of  soybeans  has  almost  tripled,  from  19  to  54  million  acres,  or 
7.7  to  22  million  hectares.  Assuming  a median  nitrogen  fixation 
rate  of  75  kg/hectare-year,  this  being  the  middle  of  the  range 
57  to  94  kg/hectare-year  quoted  by  Evans  and  Barber  (1977),  leads 
to  a change  in  biological  nitrogen  fixation  resulting  from  this 
change  of  U.S.  soybean  production  of  0.6  Mt(N)  to  1.7  Mt(N) — 
numbers  which  are  not  insignificant,  but  still  are  less  than  4 
percent  of  industrial  or  combustion  fixation  taken  separately. 

Nitrogen  fixation  rates  in  ke/hectare-year  are  57  to  94  for 
soybeans,  and  128  to  600  for  alfalfa,  with  other  legumes  gener- 
ally lying  between  these  limits  (see  Evans  and  Barber,  1977). 

The  U.S.  acreage  of  alfalfa  has  remained  approximately  constant, 
at  27  million,  which  would  correspond  to  an  annual  rate  of  ni- 
trogen fixation  of  (1.3  to  1.6  x 10^  kg);  the  acreage  of  other 
clovers  (timothy,  lespedeza,  etc.)  has  been  declining  slightly, 
and  both  the  acreage  and  the  nitrogen  fixation  is  lower  than 
those  for  soybeans  and  alfalfa.  I do  not  have  any  statistics 
on  non-U, 3.  acreage  of  these  crops,  but  the  U.S.  is  the  biggest 
factor  here. 

— * 

Dr.  A.S.  Newman  (USDA)  suggested  this  additional  source,  and 
Mr.  Floyd  Rolf  (USDA)  kindly  provided  the  acreage  figures. 

** 

Chinese  acreage  of  soybeans  has  not  been  changing  rapidly,  from 
8 million  hectares  in  1970  to  9.2  in  1976.  Brazil,  the  second 
largest  exporter  of  soybeans  (after  the  U.S.)  had  6.4  million 
hectares  planted  in  1976,  up  from  only  1.3  in  1970. 
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F . 3 AN  UPPER-BOUND  ESTIMATE  OF  OZONE  REDUCTION  IN  1 955  TO  1975 

DUE  TO  TIME-DEPENDENT  ANTHROPOGENIC  NITROGEN  FIXATION  IN 

THIS  TIME  PERIOD 

In  1975,  a reasonable  estimate  of  annual  nitrogen  fixation 
due  to  commercial  fixed  nitrogen  production  is  49.5  Mt(N)  (see 
Table  F-2),  with  an  additional  source  of  strength  29  Mt(N)  due 
to  combustion  (see  Section  F.2.2),  During  the  1955  to  1975  time 
frame,  the  growth  rate  of  commercial  fixed  nitrogen  production 
was  10  percent  per  annum,  while  that  of  fossil  fuel,  and  so  pre- 
sumably also  of  the  combustion  source,  was  5 percent  per  annum. 

Thus,  in  1955,  the  nitrogen  fixation  due  to  commercial  pro- 
duction was  7.4  Mt(N),  and  that  due  to  combustion  was  10.9  Mt 
(N).  During  the  20-year  time  period  1955  to  1975,  the  total 
nitrogen  fixation  due  to  commercial  production  was  approximately 
440  Mt(N)  and  that  due  to  combustion  was  370  Mt(N),  or  an  in- 
crease in  production  of  450  Mt(N)  over  the  rate  of  production  at 
the  initial  (1955)  rate. 

In  Table  F-4  we  list  best-estimate  nitrogen  fixation  rates 
during  this  time  period,  based  on  the  present  assumptions,  to  be 
inserted  in  numerical  models. 


TABLE 

F-4.  ESTIMATE  OF  ANTHROPOGENIC  NITROGEN 
RATE,  1955  TO  1975  (units:  Mt(N)/yr) 

FIXATION 

1955 

1960 

1965 

1970 

1975 

Commerci al 

Production  7.4 

11.9 

19.1 

30.7 

49.5 

Comb  us  ti on 

10.9 

14.0 

17.8 

22.7 

29.0 

Total 

18.3 

25.9 

36.9 

53.4 

78.5 

^Effect  of 

soybeans  and  other 

leguminous 

crops  (see 

Section 

F.2.3)  is  not  considered. 

Let  us  now  ask  for  the  anthropogenic  production  of  nitrous 
oxide  due  to  this  nitrogen  fixation.  Allowing  for  a denitrifi- 
cation factor  of  6 percent  (see  CAST,  1976,  p.  15),  this  being 
the  fraction  of  fixed  nitrogen  that  goes  to  make  N^O  rather  than 
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N on  denitrification,  indicates  an  anthropogenic  production  of 
49  Mt(N)  as  N^O  during  this  time  period,  of  which  27  Mt(N)  rep- 
resents the  increase  over  production  at  the  1955  rate,  which  we 
may  call  the  time-varying  part.  This  may  be  compared  with  a 
total  atmospheric  burden  of  nitrogen  in  N^O  of  1700  Mt(N)  (see 
Hahn  and  Junge,  1977).  Thus,  assuming  uniform  mixing,  the  frac- 
tion of  N^O  that  is  of  anthropogenic  and  time-varying  origin  due 
to  this  20-yr  injection  is  0.015. 

Now  the  predicted  effect  on  ozone  due  to  a doubling  of 

stratospheric  No0  is  a decrease  between  0 and  9 percent  (see 

d * 

Crutzen  and  Howard,  1977).  If  the  change  in  ozone  is  a linear 
function  of  the  change  in  I^O,  the  predicted  maximum  decrease  in 
ozone  due  to  a 1.5  percent  increase  in  NpO  is  0.14  percent. 

These  estimates  represent  an  upper  bound  insofar  as  the 
times  for  denitrification,  t^,  for  mixing  in  the  troposphere,  t^; 
and  for  formation  of  NO  by  reaction  (F.la),  t ^ , have  all  been 
ignored  and  also  that  none  of  the  N^O  formed  as  a result  of  an- 
thropogenic activity  is  lost  to  unknown  sinks.  In  Section  F.l 
we  see  that  (t^  + t^)  ^ 10  yr,  so  that  the  predicted  effects 
would  occur  with  this  delay,  and  correspondingly  averaged  out 
over  a longer  time  and  so  even  smaller. 

F . 4 CHANGES  IN  ATMOSPHERIC  N90  DURING  THE  1955  TO  1975 

TIME  FRAME 

If  the  effects  on  ozone  of  anthropogenic  nitrogen  fixation 
during  the  1955  to  1975  time  frame  are  to  be  significant,  one 
would  certainly  expect  to  detect  a change  in  atmospheric  nitrous 

*Note,  however,  that  at  least  some  models  using  current 
(January  1978)  chemistry  predict  that  an  increase  in  strato- 
spheric N2O  will  lead  to  an  increase  in  stratospheric  ozone 
(F.  Luther,  private  communication  to  R.C.  Oliver,  January 
1978).  I have  not  seen  published  references  showing  this  re- 
sult. It  is  of  interest  also  that,  using  CIAP  (1974)  chem- 
istry, the  predicted  effect  on  ozone  due  to  a doubling  of 
stratospheric  NO  was  a 19  percent  decrease  Lsee  Vupputuri 
(1974)].  2 
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oxide.  Pierotti  and  Rasmussen  (1977)  have  reviewed  the  atmo- 
spheric measurements  of  N_0  between  1949  and  1976.  While  there 
are  differences  between  different  observers  and  between  differ- 
ent methods  of  measurement,  and  while  there  has  been  no  inter- 
comparison of  calibrations  between  the  different  experiments, 
there  is  no  clear  evidence  of  any  long-term  trend  over  the 
period  as  a whole  . 

The  following  points  should  be  noted: 

• Local  measurements  taken  both  in  Germany  and  in  Massa- 
chusetts over  the  period  1966  to  1968  show  a definite 
increase  of  perhaps  20  percent  over  this  2 to  3 year 
period  [see  Schutz  et  al.  (1970),  Hahn  and  Junge  (1977)], 
but  no  comparable  effect  is  apparent  in  the  rest  of  the 
data . 

• H.  Craig  and  R.F.  Weiss  [see  report  on  GEOSECS,  Science , 
195 , 166  (1977),  also  Denitrification  Seminar  (1977)] 
compared  shipboard  samples  of  N^O  taken  at  various  times 
and  find  that  the  average  concentration  for  equivalent 
sites  in  the  Pacific  had  increased  by  4 ppb/decade,  or 
0.15  percent  per  annum. 

• Goldan  et  al.  (1978),  in  reporting  even  more  recent  mea- 
surements, indicate  that  measurements  of  N^O  will  be  made 
on  a continuing  basis;  from  their  careful  measurements 
(which,  however,  only  extend  over  a 6-month  time  frame), 
they  are  able  to  rule  out  changes  in  M^O  greater  than  4 
percent  per  annum. 

In  conclusion,  it  must  be  stressed  that  the  present  discus- 
sion is  limited  to  the  1955  to  1975  time  frame;  it  appears  that 
current  anthropogenic  injections  of  fixed  nitrogen  are  probably 
too  small  to  have  a detectable  effect  on  stratospheric  ozone, 
but  this,  clearly,  could  not  continue  to  be  the  case  indefi- 
nitely if  nitrogen  fixation  rates  continue  to  grow  exponentially. 
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